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Molecular weights have been compared by combination of the x-ray and density data for 
the crystals, calcite, diamond, lithium fluoride, sodium chloride, and potassium chloride. By 
assumption of the chemical atomic weights of carbon, chlorine, potassium, lithium, and sodium, 
values for the atomic weights of calcium and fluorine have been calculated. Thé method 
employed consists of a series of independent calculations from the results of which final weighted 
average values have been taken for the atomic weights of calcium and fluorine. This has 
yielded 40.0849 +0.0030 and 18.9967 +0.0013, respectively. These values are most probably 
the most precise values for these constants at the present time. From the results of the present 
work it becomes apparent that the determination of atomic weights from comparison of molec- 
ular weights by combination of density and x-ray data is as reliable as other standard atomic 


weight procedures. 


I. INTRODUCTION 


OLECULAR weights have been compared 
previously'~* by combination of density 
and x-ray data. For the comparison of the 
molecular weights of two crystalline substances, 
there must be known the values of the absolute 
densities, the ratio of the true x-ray grating 
spaces, and the arrangement of molecules in the 
unit cell of each substance. The comparison can 
be expressed by means of the relationship, 


M2{ (01/p2) FR,*}. (1) 


The numerical subscripts refer to the two 
crystalline substances, the molecular weights of 


* Research Associate, State Engineering Experiment 
Station, Georgia School of Technology, Atlanta, Georgia. 
'C. A. Hutchison and H. L. Johnston, J. Am. Chem. Soc 
63, 1580 (1941). 
(1943) L. Johnston and D. A. Hutchison, Phys. Rev. 62, 32 
°C, A. Hutchison, J. Chem. Phys. 10, 489 (1942). 
‘D. A. Hutchison, Phys. Rev. 66, 144 (1944). 


which are compared; M is a molecular or atomic 
weight; p is a density ; F is a factor determined by 
the crystal structures; and R, is the ratio of the 
true x-ray grating spaces of the two crystals. 
C. A. Hutchison and H. L. Johnston! have used 
Eq. (1) to compare the molecular weights of 
lithium fluoride and calcite. With the assumption 
of the atomic weights of calcium, carbon, and 
lithium, they were able to obtain an atomic 
weight for fluorine. H. L. Johnston and the 
writer? have obtained an atomic weight for 
fluorine utilizing the data on lithium fluoride and 
sodium chloride and assuming the atomic weights 
of lithium, sodium, and chlorine. C. A. Hutchison* 
has calculated the atomic weight of calcium 
utilizing the data on calcite and diamond and 
assuming the atomic weight of carbon. From this 
value for calcium and that assumed for carbon, 
he has calculated the atomic weight of fluorine. 
Recently, the writer reported values for the 
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TaBLeE I. Density and x-ray data corrected to 20°C for certain selected crystals 
used in the comparison of molecular weights. 


Substance 


Data 
designation - 


P2990 (8/ 


(degrees) 5 X106 n plane 


Reflecting 


a X105(linear) 


2.640302 
+0.00014 


A LiF 


2.71030° 
0.00005 


B Calcite 


3.51540° 
0.00006 


2.163604 
+0.00004 


1.98827¢ 
+0.00004 


Diamond 


NaCl 


KCl 


6°42'34.9" 0.30." 1.85 1 
27°51'32.1+0.3° 4 


23°26'5.8” +0.59 2.39 1 100 0.118! 


6°28'25.8" +0.49 1.32 1 


83°51/15.7"" +0.5/ 7.91 3 111 3.417 


Natural 1.04* 
Cleavage 


15°51’4.4”+0.8" 1.41 1 100 4.05" 


100 3.77" 


my Rev. 25, 618 (1925). 
. A. Bearden, Phys. Rev. 54, 698 (1938). 
. L. Johnston and D. A. Hutchison, Phys. Rev. 62, 32 (1942). 
e °D. A. Hutchison, Phys. Rev. 66, 144 (1944). 


Cc. A. Hutchison and H. L. Johnston, J. Am. Chem. Soc. 62, 3165 (1940). 
+ J. A. Bearden, Phys. Rev. 38, 2089 (1931). Bearden gives this density value as the best mean of his value and that of Defoe and Compton, 


/ Straumanis, levins and  Someree, Zeits. f. physik. Chemie B42, 143 (1939), 


9 Y. Tu, Phys. Rev. 40, 622 (1932). 

4 Larsson, Phil. Mag. [7] 3 1136 (1927). 

* M. Siegbahn, Ann. d. physik BD59, 56 (1919). 
3 Adenstadt, Ann. d. physik 26, 69 (1936). 

* Larsson, Inaugural- issertation, Uppsala (1929). 


! International Critical Tables (McGraw-Hill Book Company, Inc., New York, 1927), Vol. 3, 
™ Straumanis and levins, Zeits. f. Physik 102, 353 (1936); Straumanis, levins, and Karlsons, Doite i anorg. allgem. Chemie 238, 175 (1938). 


for calcite, lithium fluoride, and potassium 
chloride. 

The R,’s used in the present work were calcu- 
lated directly from the primary data, the x-ray 
reflection angles, as were those of C. A. Hutchi- 
son.* Since the values of the R,’s, F’s, and p’s 
have been determined with high precision for the 
five crystalline substances’ used in the present 
calculations, it is possible to compare the mo- 


lecular weights of these substances with high 


precision. 

If sufficient x-ray data were available, it would 
be possible to compose a system of x-ray-density 
atomic weights which would be independent of 
chemical determinations of atomic weights. Since 
a set of crystals cannot be found which would 
give an atomic weight system based on oxygen, 
and independent of atomic weights from other 
sources, we may refer to the results of the present 
calculations as x-ray-density-chemical atomic 
weights. 

The purpose of the present paper is the presen- 
tation of the method of calculating and the 
results obtained for the atomic weights of calcium 


» Klemm, Tilk, and Mullenheim, Zeits. f. anorg. allgem. Chemie 176, 1 (1928). 


atomic weights of calcium and fluorine which 
were calculated from the x-ray and density data 


and fluorine utilizing the x-ray and density data 
for the five crystalline substances, calcite, dia- 
mond, lithium fluoride, sodium chloride, and 
potassium chloride. 


TABLE II. Chemical atomic weights employed in 
the calculations. 


Symbol Value 
12.0104+0.0002¢ 
35.4570+0.0005° 
K 39.0960 +0.0005¢ 
LA 6.9390 +0.00024 
Na (1)22.9970+0.0005° 


« This value is a mean of all the good combustion analyses of Baxter 
and Hale and is in good agreement with that of Scott and Hurley. 
(1) Baxter and Hale, J. Am. Chem. Soc. 59, 506 (1937). 
(2) Scott and Hurley, J. Am. Chem. Soc. 59, 1905 (1937). 
+’ (1) Richards and Wells, “‘A Revision of the Atomic Weights of 
Sodium and Chlorine,’’ Carnegie Institution of Washington, Publication 
No. 28 (1905). 
(2) Noyes and Weber, J. Am. Chem. Soc. 30, 13 (1908). 
(ene 4 Zintl and Meuwsen, Zeits. f. anorg. allgem. Chemie 136, 223 
a Honigschmid, Chan, and Birckenbach, Zeits. f. anorg. allgem. 
Chemie 163, 315 (1927). 
(5) Honigschmid and Sachtleben, Zeits. f. anorg. allgem. Chemie 
178, 1 (1929). 
(6) Scott and Johnson, J. Phys. Chem. 33, 1975 (1929). 
(7) Honigschmid, J. Am. Chem. Soc. 53, 3012 (1931). 
¢ (1) Richards and Stahler, J. Am. Chem. Soc. 29, 623 (1907). 
(2) Richards and Mueller, J. Am. Chem. Soc. 29, 639 (1907). 
(3) Baxter and MacNevin, J. Am. Chem. Soc. 55, 3185 (1933). | 
(4) Honigschmid and Sachtleben, Zeits. f. anorg. allgem. Chemie 
213, 365 (1933). 
‘6 ) McAlpine and Bird, J. Am. Chem. Soc. 63, 2960 (1941). 
4(1) Richards and Willard, J. Am. Chem. Soc. 32, 4 (1910); Zeits. f. 
anorg. allgem. Chemie 66, 229 (1910). 
(2) Carnegie Institution Publication No. 25, 1 (1910). 
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WEIGHTS OF CALCIUM AND FLUORINE 


TABLE III. The atomic weight of calcium. 


Calculation 
No. 


Data for mol. 


wt. comparison 


from Table I 


F 


Ro 


Ry 


Chem. at. wt. 
on which 
calc'd. value 
is dependent 


M calcite(B ) 
Mnaci(D) 


Mealcite(B) 
Mxci(E) 


Mealcite(B) 


Mdiamond(C) 


calcite(B) 
Miir(A) 


1.09594 
+0.00001 


1.09594 
+0.00001 


0.273985 
+0.000003 


1.09594 
+0.00001 


1.076443 
+0.000020 


0.965029 
+0.000020 


3.404303 
+0.000046 


1.508187 
+0.000021 


0.4604078 
+0.0000021 


1.000000 
+0.000014 


1.000000 
+0.000014 


0.4604078 
+0.0000021 


Na, Cl, C 


Kk, Cl, C 


Most probable value (weighted average) of atomic weight of calcium—40.0849 +0.0030. 


TABLE IV. The atomic weight of fluorine. 


Calculation 
No. 


Data for mol. 
wt. comparison 
from Table I 


F 


Ro 


Ry 


Chem. at. wt. 
on which 
calc’d. value 
is dependent 


5 


Muir(A) 
Mxci(£) 


Muir(A) 
Mnaci(D) 


Mi ir(A) 


Mdiamond(C) 


Muir(A) 


1.000000 
+0.000000 


1.000000 
+0.000000 


2.000000 
+0.000000 


0.9124587 
+0.0000083 


0.639860 
+0.000023 


0.713733 
+0.000012 


1.128610 
+0.000047 


0.663048 
+0.000015 


2.1719872 
+0.0000098 


1.000000 
+0.000018 


2.1719872 
+0.0000098 


2.1719872 


K, Cl, Li 


Na, Cl, Li 


Li, C 


Na, Cl, C 


Mealcite(B) 


+0.0000098 Li, K 


Most probable value (weighted average) of atomic weight of fluorine—18.9967 +0.0013. 


II. DATA FOR THE MOLECULAR WEIGHT 
COMPARISONS AND ATOMIC 
WEIGHT CALCULATIONS 


The quantities necessary for the molecular 
weight comparisons have been corrected to 20°C 
and the values listed in Table I. The density of 
each substance at 20°C has been computed from 
the most probable reported values by means of 
the cubical coefficient of expansion, which is three 
times the linear expansion coefficient listed in 
column eight of Table I. The reflection angle, #, 
for each crystal at 20°C has been obtained from 
the most probable literature value by means of 
the expression, 


log sin =log sin 3;— (a(20—#t)/2.303), (2) 


where ¢ is the temperature at which the reflection 
angle is reported in the literature and a is the 
linear coefficient of expansion given in column 


eight of Table I. 6, the unit decrement of the 
index of refraction of a crystal, is defined by the 
expression, 6=1—jy, where u is the index of re- 
fraction. The values of the 6’s given in column 
five of Table I were calculated by means of the 
expression, 


5 (3) 


which assumes the wave-lengths used to be much 
shorter than any critical absorption wave-length 
of the crystals involved. n is the total number of 
electrons per cubic centimeter, and other symbols 
have their usual meaning. Column six of Table I 
gives the order of the reflection angles and column 
seven, the reflecting plane. 

In the calculations of the atomic weights, 
certain chemical atomic weights have been as- 
sumed. The values for these chemical atomic 
weights are listed in Table I]. 
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In all instances of these original data, reported 
values were surveyed critically, and, where neces- 
sary, recalculations from original data were made 
and arbitrary weights, based on the internal and 
external consistency of the reported data, were 
assigned to obtain most probable values. __ 


III. CALCULATIONS OF THE ATOMIC WEIGHTS 


In considering the best method of attacking 
the problem of getting the most probable values 
for the atomic weights of calcium and fluorine 
from the available x-ray and density data, two 
methods are at one’s disposal. 

One method involves a series of stepwise com- 
putations from which each atomic weight ob- 
tained is dependent on the preceding computa- 
tion. If there are given the densities, p’s, and the 
x-ray reflection angles, #’s, for the five crystalline 
substances given in Table I, any four of the 
corresponding molecular weights may be com- 
puted from the remaining one. From certain sets 
of three assumed chemical atomic weights (no 


less than three atomic weights can be chosen), the — 


remaining four atomic weights can be computed. 
This method cannot be used since, in any given 
series of stepwise calculations to obtain best 
values for the atomic weights of calcium or 
fluorine, the final probable error is so great that 
the purpose of the calculations is defeated. 

Another method consists of a series of inde- 
pendent calculations, from which a final weighted 
average value can be obtained. The latter is the 
method chosen for the calculations presented 
here. The results of each independent calculation 
are presented in Tables III and IV. 

Let us take as an example of these calculations 
the computation of the atomic weight of fluorine, 
listed as computation number 5 in Table IV, 

which utilizes the x-ray and density data for 
lithium fluoride (A) and potassium chloride (Z). 
Equation (1) becomes in this case, 


My= — Mii. (4) 
PKC1 


The R, has been computed by means of the 
expression, 


mz sin 3;(1—6;/sin? #1) 


(5) 
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where the subscripts ‘‘1’’ refer to data on lithium 
fluoride and the subscripts ‘‘2”’ refer to data on 
potassium chloride. ~ is a factor which corrects 
for the fact that x-ray reflection angles are not 
from corresponding planes in the two crystals; 
R,5 is the ratio of the wave-lengths employed in 
the determination of x-ray reflection angles; the 
other symbols have their usual meaning as given 
in the introduction. 

Other data employed in Eq. (4) are listed in 
Table I. Since both crystals have the sodium 
chloride type structure, the value of F becomes 
unity. It is seen in columns six and seven of 
Table IV that the calculated atomic weight for 
fluorine, 18.9967+0.0031, is dependent on the 
chemical atomic weights of chlorine, potassium, 
and lithium. 

The results of the series of independent calcula- 
tions for calcium are given in Table III, and those 
for fluorine in Table IV. Their probable errors 
have been obtained in the usual manner by 
application of the law for propagation of errors, 


dU? 
(—) p2+(—) (6) 
Ox oy 


where U=f(x, y, ++). 

In order to utilize the x-ray and density data 
for lithium fluoride in the determination of the 
atomic weight of calcium (see calculation number 
4 of Table III), a weighted average was taken for 
the atomic weight of fluorine from the results 
obtained in the first three computations of 
Table IV. This result, 18.9967+0.0015, is not 
dependent on an atomic weight of calcium. This 
procedure makes possible the use of all available 
data and reduces the error of the final adopted 
value for the atomic weight of. calctum by ten 
parts per million. 

In a similar manner the data on calcite have 
been used to obtain an atomic weight for fluorine 
(see calculation number 8 of Table IV). The 
weighted average, 40.0848+0.0034, obtained 
from the results of the first three calculations of 
Table III, has been used in calculation number 8. 
Thus, in this case all data have been employed 
and the final adopted value for the probable error 

5 Ry was calculated from the reflection angles for the 
CuKa; line and MoKa; line from calcite as obtained by 


J. A. Bearden and C. H. Shaw, Phys. Rev. 48, 18 (1935) 
and Larsson, Phil. Mag. ois 1136 (1927). 
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in the atomic weight of fluorine is reduced by nine 
parts per million. 

The writer has observed from results obtained 
in carrying out these computations that x-ray 
angles from “‘powder”’ method data are definitely 
less trustworthy than angles obtained by the 
Bragg reflection method even though great im- 
provements have been made on the “powder” 
method by Straumanis and co-workers.® For this 
reason ‘‘powder”’ data have been excluded from 
this work except where they were the only data 
‘ available. 

A possible source of error® in these computations 
lay in the uncertainty in the thermal coefficients 
of expansion of the crystals involved. With re- 
duction of density and x-ray data to 20°C, it was 
necessary to make correction for the uncertainties 
of the thermal coefficients. Each coefficient was 
considered separately; but in all, the correction 
amounted to multiplying the probable error 
assigned by the original investigator by a factor 
of approximately 2. 

With this description of the calculation pro- 
cedures, one may consult the summary of results 
of individual calculations in column six of 
Tables III and IV. 

The final most probable value for each 
atomic weight is a weighted average based on 
the probable errors of the results of the individual 
calculations. The final adopted value for calcium 
is 40.0849+0.0030 and for fluorine is 18.9967 
+0.0013. The probable errors are based on the 
internal consistency of the results of the indi- 
vidual calculations. For calcium the error amounts 
to 75 parts per million and for fluorine to 68 parts 
per million. These appear to be definite improve- 
ments over previously published values. The 
probable error obtained from the external con- 
sistency of individual results can be seen to be 
very small compared to that obtained from the 
internal consistency. The adopted probable error 
is the larger of the two. Since theory predicts that 
the ratio of these two values should not differ 


6 This observation has been made previously by R. T. 
Birge in connection with his work on the value of the 
on number. R. T. Birge, Am. J. Phys. 13, 63-67 


387 


from unity by more than about 28 percent, it can 
be concluded that too great allowance has been 
made for systematic errors in the original data 
and that the adopted probable error is somewhat 
greater than the true probable error. 


IV. CONCLUSIONS 


If the density and x-ray data are utilized to 
obtain other atomic weights, one obtains excellent 
checks between the chemical atomic weights and 
the calculated values. In view of these close 
checks between chemical and calculated atomic 
weights as well as the method employed in these 
calculations, considerable weight can be attached 
to the results for the atomic weights of calcium 
and fluorine obtained here. 

It will be seen that the value obtained for 
calcium, 40.0849+0.0030, is in close agreement 
with the value 40.085, obtained by Honigschmid 
and Kempter.’ It is somewhat higher than the 
recent value, 40.080+0.005, given by Birge* and 
is in disagreement with the older value of 
Richards and Honigschmid® and the more recent 
value of Kendall, Smith, and Tait.’ 

The atomic weight, 18.9967 +0.0013, obtained 
for fluorine is in close agreement with the gas- 
density value, 18.997, of Cawood and Patterson" 
using fluoroform, but in disagreement with their 
value, 18.995, using carbon tetrafluoride. It is in 
disagreement with the gas density value of Moles 
and Toral,” using silicon tetrafluoride and in 
disagreement with the mass_ spectrographic 
value, 18.999.!8 

It is apparent that the determination of atomic 
weights from the comparison of molecular weights 
by combination of density and x-ray data is as 
reliable as other standard atomic weight pro- 
cedures. 


7 Honigschmid and Kempter, Zeits. f. anorg. allgem. 
Chemie 195, 1 (1931). 

8R. T. Birge, Rev. Mod. Phys. 13, 233 (1941). 
we and Honigschmid, J. Am. Chem. Soc. 32, 1577 

10 Kendall, Smith, and Tait, Nature 131, 688 (1933). 

11 Cawood and Patterson, Trans. Roy. Soc. (London) 
A236, 77 (1936). 

12 Moles and Toral, Zeits. f. anorg. allgem. Chemie 236, 
225 (1938). . 

18 A. K. Brewer, Phys. Rev. 49, 867 (1936). 
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Force Constants of Triatomic Molecules 


Gro. GLOCKLER AND Jo-YuN TuNnG* 
Department of Chemistry and Chemical Engineering, State University of Iowa, Iowa City, Iowa 
(Received June 25, 1945) 


In the general expression of the potential energy of the 
triatomic molecule X YZ there are six force constants ¢1, ¢1’ 
(valence bond), cz (angle deformation), c3, ¢3’ (angle bond 
interaction), and cy, (bond-bond interaction), of which 
only three can be calculated from the three experimentally 
determined fundamental frequencies. For the isosceles 
triatomic molecule YX2 the number of force constants 
reduces to four (¢1, ¢2, c3, and cy). From a study of all the 
real values of these four force constants, it is seen that they 
are related to one another by ellipses and straight lines. 
The graphs of these ellipses indicate the possibility of 
choosing a certain singular point (near one end of the major 
axis) which gives the following relation between the force 


constants: 
dc4/dc2=0 or dc,/dc3=0. 


This condition and the equations usually obtained from the 
theory of small vibrations permit the calculation of a 
singular set of four force constants of the isosceles triatomic 
molecule. When these force constants are applied to isotopic 


molecules, satisfactory checks between the calculated and — 


experimentally determined frequencies are obtained. 

The ellipses showing the relation between the force 
constants indicate clearly why the simple central force field 
(c1, C2, 3) or the simple valence force field (c1, C2, ¢3), con- 
taining one interaction constant or crossterm, must yield 
imaginary values for cz and c; in certain cases. 


INTRODUCTION 


HE theory of small vibrations’? leads to 

equations connecting the fundamental fre- 
quencies of a molecule and the force constants, 
masses, distances, and angles of its structure. 
These equations can be used to calculate the 
force constants if other quantities are known. 
It happens quite frequently that fewer fre- 
quencies are available than the number of force 
constants which it is desirable to know. The 
triatomic molecule has only three fundamental 
frequencies, and, hence, only three force con- 
stants can be determined. Two cases are usually 
dealt with: the central force field and the 
valency force field. Since the latter views of 
atomic interaction agree with the ideas of chemi- 
cal valence, they are used most frequently and 
also in the following considerations. For the 
general triatomic molecule XYZ the six con- 
stants are: valence bond (¢;, ¢:’); angle bending 
(c2); angle-bond interaction (c3, c3’); and bond- 
bond interaction (c4). For the isosceles triatomic 
molecule the four force constants are ¢1, Cs, C3, 
and cy. If the constant cs is neglected (c4=0), 


*A thesis submitted in partial fulfillment of the re- 
quirements for the degree of master of science, in the De- 
rtment of Chemistry and Chemical Engineering, in the 
raduate College of the State University of Iowa. 
1K. W. F. Kohlrausch, Der Smekal-Raman Effekt 
(Verlagsbuchhandlung, Julius S ringer, Berlin, 1931). 
141 “The Raman Effect,” Rev. Mod. Phys. 15, 


the case can be called the simple valency force 
field and these three constants can be 
determined from the three fundamental fre- 
quencies. However it happens quite often that 
no consistent set (ci, C2, ¢3) can be found, and cz 
and c; turn out to be imaginary. Any new method 
which would furnish additional information per- 
mitting the evaluation of the four force constants 
of the isosceles triatomic molecule, even though 
only three frequencies are available would be of 
some value and should be investigated. Such an 
empirical method is considered in this study. 
The process is as follows: All of the real values 
of the four force constants are calculated by 
allowing c, to have all possible values which do 
not lead to imaginary numbers for cz and ¢;. 
These real values are shown in Figs. 1-6 and 8, 
as ellipses and straight lines. It is seen that an 
infinite number of sets (¢1, C2, €3, cs) can be found 
that satisfy the initial equations of the theory of 
small vibrations. However, it is felt that one of 
these sets might be distinctive, and the essence 
of this study is to find this special case. It is 
expected that the bond force constant ¢; in a 
given molecule should be of the same order of 
magnitude as the force constant of the corre- 
sponding diatomic molecule. For example the 
bond constant of OH in the water molecule may 
be expected to be similar in magnitude to the 
one for the same atoms in hydroxyl radical. 

In the case of water.and deuterium oxide, it 
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was found that the ellipses obtained from the 
respective frequencies are nearly alike (Fig. 1, 
full and dotted lines).? Since c, is an interaction 
constant it may be surmised that its value should 
be small in magnitude. From an inspection of 
these ellipses it appears that points near one of 
the ends of the major axes might offer a satis- 
factory choice of one distinctive set of all four 
force constants from the infinite number of sets. 
One such singular point is given by the condition 


or dc,/dc3=0. (1) 


It turns out that this choice also yields a value 
for c: which seems reasonable. The other similar 
point near the other end of the major axes refers 
to a much larger absolute value of c, and a value 
for c; which is very different from the force 
constant c; of the related diatomic molecule. 
Hence this set is disregarded. The program of this 
study is then to calculate all real sets of four force 
constants for triatomic isosceles molecules from the 
three equations of the theory of small vibrations 
using the three known frequencies and then applying 
condition (1) to choose one distinctive set. The 
force constants so determined are then used to 
calculate the three frequencies of isotopic mole- 
cules which can be compared with experiment. 


THEORY OF SMALL VIBRATIONS 


Only a few salient features will be mentioned 
since this topic has been discussed often.'? The 
potential energy expression (V) of the general 
triatomic molecule (XYZ) can be written in 
terms of cartesian displacements: 


where 


2V Ao? Agi? Aa? 
Aa: + 731° Aa: Asi 
Asi. (2) 


The force constants have been defined above and 


Aoi, A31=small bond displacements, 
Aa=small angle deformation, 
731 = equilibrium bond lengths. 


From Eq. (2) and from the kinetic energy (T) of 
the system 

2T=> i 
it is possible to write the Lagrangian equation 


of motion 


where L=T—YV, q; is any coordinate and the 
m; refers to the masses of the atoms. From 
these relations determinantal equations are de- 
rived from which force constants can be cal- 
culated. 

The determinant. of the general triatomic 
molecule is given for comparison with the special 
cases of the isosceles structure, the linear sym- 
metrical molecule and the linear asymmetrical 
case. Some examples of each type have been 
examined. 


GENERAL TRIATOMIC MOLECULE 
The determinantal equation can be written 
A—Ly 
A—L3s 


(3) 


—€1/ port (ro1/731)! sin a— (¢4/my,) COS a, 


—(¢1/my) cos a+ sin a—C4/ ps1, 


—L13= (¢1/my) sin a— (1/1) + (1/2) (731/721) 


+(1/mts) (r21/731) ]+(C4/m1) (131/721)! sin a, 


cos a+(c3’/m1) (ro1/r31)4 sin a—C4/p21, 


(131/721)! sin a= (c4/my) COS 


sin a —cs'[(1/m1) + (1/me2) (131/721) 


+(1/ms)- (ro1/rs1) ]+ (ca/m,) (r21/rs1)4 sin a, 
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—L31= sin a—C3/p21— Cos a, 
— sin a— (¢3/m1) cos a—C3'/ua1, 
—L33= — Col (1/mr) + (1/me) (731/721) + (1/ms) (121/731) J 

+ (¢3/m1) sin a+ (731/721)! sin 


where 1/p21=1/m2+1/mi, 1/u31=1/m3+1/m,; m, is the mass of the central atom. Expanding this 
determinant and comparing it with the general cubic equation yields 


, 


Only three force constants contained in the L;;'s —Lis= —Le3=(¢:/M) sina 
can be determined assuming that the valence , 
angle a and the bond distances and r3; can be M) sin a, 
obtained from other sources such as data on —L3,= —L32=(c2/M) sin a—c;[1/u+cos a/M], 
electron diffraction. The weights of the atoms are wae " ‘ 

labeled as follows: Y(m, or M); X(mz or m); and —2es[1/n—cos M sin a. 


Z(ms3 or m). From the determinantal equation, it is seen that 


ISOSCELES TRIATOMIC MOLECULE 


In this case, c1=c1', C3=C3', M1=M, Are = (5) 
M2=m3=m, and 1/u=1/M+1/m. The quantities 
L;; reduce to 
ij 
With the use of Eqs. (5) and the special condition 
(Eq. (1)), the four force constants of water 
+(c3/M) sin a—(c4/M) cosa, molecule have been determined from its three 
fundamental frequencies. Since the force con- 
stants of deuterium oxide are the same as for 
+(c3/M) sina—cs/u, the water molecule, they can now be used to 


—Ly= 


TABLE I. Calculated frequencies of deuterium oxide TABLE III. Experimental frequencies of isosceles 
vapor molecule. triatomic molecules. 


Experiment (cm~)¢ Calculation (cm~!) Molecule Angle a wi (cm™) w2 (cm~!) w3 (cm™~) 


w3 2883.79 2883.51 =104°31’ 3825.32) 1653.91 3935.39 

1 2758.06 2758.76 H2S 92°20’ 2615 1265 2630° 

we 1210.25 1209.99 >CHz 115° 2970 1444 3000¢ 
SO: 120° 1152 524 1361° 


«2D. M. Dennison, Rev. Mod. Phys. 12, 175 (1940). 
«D. M. Dennison, Rev. Mod. Phys. 12, 175 (1940). 
+ See reference 6. 

TABLE II. Real values of the force constants of the water ¢ W. G. Penney and G. B. B. M. Sutherland, Proc. Roy. Soc. A156, 


vapor molecule (megadynes/cm), (i= imaginary). 654 (1936). 


TABLE IV. Calculated force constants of isosceles triatomic 
molecules. (megadynes/cm). 


Molecules ce C3 


0.0460 


Ss 


ll 


c 

t 

f 

t 

( 

a 

0.8485 i i 0 i i 

0.7485 0.1797 0.1465 0.1000 0.1650 —0.1283 H.0 0.8386 0.0763 0.0234 —0,0099 
0.6485 —0.2000 0.1770 —0.1516 H.S 0.3942 0.0060  —0.0016 
0.5485 0.3689 0.3430 0.1216 —0.1076 123 
05017 0.4011 04011 03468 0.0044 0.0044 >CHe 0.4874 0.0562 0.0170 0.0 
0.4985 i i —0.3500 i i SO2 1.0416 0.1069 0.1346 0.0479 
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Fic. 1. Relation between 
force constants of water and 
deuterium oxide (mega- 
dynes/cm). (Dotted D,O, 
and full line H,O). 


Fic. 2.4Relation between 
force constants of hydrogen 
sulfide molecule (mega- 
dynes/cm). 


calculate the fundamental frequencies of deu- 
terium oxide as shown in Table I. It is seen that 
the calculated values check the experimental 
frequencies most satisfactorily. It appears that 
the smallest absolute value of cs (— 0.0099) repre- 
sents a satisfactory singular point, being near 
one end of the two major axes of both ellipses. 
This special point is represented by condition 
(Eq. (1)). The other extreme of the two major 
axes (c4= —0.3468) would yield a value for c; 
very markedly different from the force constant 


TasBLe V. Comparison of experimental and calculated fre- 
quencies of isosceles triatomic molecules. 


Cl CL 
(em) (cm™) 


2818.60 1403 1407.41 — 3 
1210.25 1209.99 2883.79 
1090 1090.99 
934 930.92 2000 
1072.00 


— - 1214.00 
1438.08 


Mole- 


wl 
Exp. 
cules 


(cm~) 


* See reference 3. 

> See reference 4. 

© See reference 5. 

41D. M. Dennison, Rev. Mod. Phys. 12, 175 (1940). 

¢ A. H. Nielsen and H. H. Nielsen, J. Chem. Phys. 5, 277 (1937). 
1G. M. Murphy and J. E. Vance, J. Chem. Phys. 6, 426 (1938). 
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Fic. 4. Relation between 
force constants of the sulfur 
dioxide molecule (mega- 
dynes/cm). 


Fic. 3. Relation between 
force constants of the 
methylene radical (mega- 
dynes/cm). 


for the diatomic hydroxyl radical (OH). From 
the above relations the following explicit ex- 
pressions of the force constants can be derived: 


], 
C2 ], 
C4=C€1—A;3/C, 


where A =1/yn+cos a/M; B=sina/M; C=1/p 
—cosa/M. A sample calculation is shown in 
Table II. It can readily be seen that the simple 
valency force field without the interaction force 
constant (i.e., cs=0) leads to imaginary values 
of the force constants cz and ¢3. 

In this manner several molecules of this type 
have been studied. The experimental frequencies 
used are collected in Table III. From them the 
force constants of these molecules were calcu- 
lated (Table IV) and their respective ellipses 
indicating the relations between their force con- 
stants, are shown in Figs. 1-4. The force constants 
mentioned in Table IV have been used to calcu- 
late the frequencies of certain isotope molecules 


(6) 
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Fic. 5. Relation between force constants of the carbon 
dioxide molecule (megadynes/cm). 


as shown in Table V. The calculation of the 
deformation frequency w2 of the hybrid water 
molecule (HDO) checks experiment satisfac- 
torily. On the other hand a discrepancy exists 
between calculation and experiment in the case 
of the symmetric frequency w:. The experimental 
value of w: is given by Barker and Sleator*® as 
2720 cm=. On the other hand calculations by 
Libby* and by Shaffer and Schuman® yield 
2818.1 cm! and 2820.3 cm for this frequency. 
These values are not based on the method dis- 
cussed here. The present calculations yield 2818.6 
cm~!. Hence these three calculations agree among 
themselves but they do not check the experi- 
mental value. It must mean that the infra-red 
band at 2720 cm needs to be reinterpreted. 
The frequency w2 of the hybrid hydrogen sulphide 
molecule (HDS) checks with experiment very 
well. For deutero-hydrogen sulphide (D2S) 1 
and we: check reasonably well although the 
difference between experimental and calculated 
values for w; needs further attention. 

The calculated frequencies of the methyl- 
ene radical (>CHg2) and its isotopic relatives 
(>CHD, >CD:) show among themselves the 
same relations as the corresponding isotopic 
water molecules (H2O, HDO, and D.O).and the 


a . a Barker and W. W. Sleator, J. Chem. Phys. 3, 660 
‘W. F. Libby, J. Chem. Phys. 11, 101 (1943). 
5 W. H. Shaffer and R. P. Schuman, J. Chem. Phys. 12, 
504 (1944). 
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similar hydrogen sulfide structures (H2S, HDS, 
and D,S). The introduction of the first deuteron 
into the ordinary proton molecule (H2X, X 
=C, O,S) has a large effect on the frequency 
(Aw; =784, 1006, 733 cm) while the sec. ' 
deuteron changes the frequency very li - 
(Aw; =46, 60, 17 cm~'). The deformation 
quency w2 varies regularly by a few hundreds 
wave numbers (first deuteron: Aw:=230, 2« , 
174 cm; second deuteron: Aws=142, 197, 159 
The asymmetric frequency shows just 
the opposite effect to w;. Introducing the first 
deuteron into the ordinary proton structure has 
only a small effect (Aw3= 16, 53, 8 cm~!) while a 
second deuteron changes w; markedly (Aw3= 753, 
1000, 733 Substituting into the CH:- 
radical does not affect either of the three funda- 
mental frequencies very much as might be 
expected (Table V). 


LINEAR SYMMETRICAL TRIATOMIC MOLECULE 


In this case, ¢1=¢1’, C3=C;’, the valency angle 
a = 180° and 732 =731 +721 = 2731 = 2re1. The quanti- 
ties L;; reduce to 


—Ly= —L2= 
—Ly= 
—Ly3= —L23= 
—L3:= —L32= —c3/m, 

—L33= 


The carbon dioxide and carbon disulphide 
molecules whose experimental frequencies are 


TABLE VI. Experimental frequencies of linear symmetrical 
triatomic molecules. 


Molecules (cm~) w2 (cm~!) (cm~!) 
COz 1321.7 667.9 2362.87 
655.0 397.0 1522.0° 


«T. Y. Wu, Vibrational Spectra and Structure of Polyatomic Molecules 
(Prentice-Hall, Inc., New York, 1941). 
+ See reference 6. 


TABLE VII. Calculated force constants of linear sym- 
metrical triatomic molecules (megadynes/cm). 


Molecules C1 c2 
CO, 1.5339 0.0571 0 0.1052 
CS: 0.7467 0.0234 0 0.0598 


392 
s 
14 n 
t nN 
12 
10 
{ 
b 
t 
b 
t 
| 
Cc 
ti 
Cc 
d 
nN 
| Ct 


FORCE CONSTANTS OF TRIATOMIC MOLECULES 


shown in Table VI have the force constants 
mentioned in Table VII. By their use the isotopic 
molecules C"O, and CS. were calculated as 
shown in Table VIII. Only the asymmetric 
ws of can be compared with 
periment. The check seems satisfactory. The 
‘' tion between the force constants are shown 
in Figs. 5 and 6. It is seen that the angle-bond 
-hteraction constant (c3) is related to the bond- 
bond interaction constant (cs) by an ellipse 
whose major axis is parallel to the cs-axis and 
that cz is equal to zero. This feature is true in 
both instances and may therefore be taken to be 
the case for linear molecules. The angle deforma- 
tion force constant C2 is related linearly to the 


constant c4. In this case the usual ellipse reduced . 


to a straight line. 


LINEAR ASYMMETRIC TRIATOMIC MOLECULES 


Of the six force constants, c3 and c;’ are taken 
to be zero on the assumption that this feature 
can be extended from carbon dioxide and carbon 
disulfide, even though the molecule (mzmms) is 
not symmetrical as for example OCS, although 
linear. The quantities L;; become: 


—L13;=0, 

= 

—Lo3=0, 

—L3=0, 

—L32.=0, 

—L33= (1/1) 

+(1/me) (31/121) + (1/ms) (721/731) J. 


The determinant can now be written as follows: 


0 
0 
0 0 A—Les 


Expanding and comparing with the general 
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TABLE VIII. Calculated frequencies of linear symmetrical 
triatomic molecules. 


we 
Cal. Exp. Cal. 
(em~) 


1321.70 - 649.01 2284.5 
655 383.96 


Mole- 
cules 


Exp. Exp. Cal. 


2295.97¢ 
1472.02 


«A. H, Nielsen, Phys. Rev. 53, 983 (1938). 
quadratic equation yields 
ho=Ls33, 
= — Liha. 


(7) 


An attempt was made to find another relation as 
described above between the force constants 
C1, ¢i’, and c, which would make it possible to 
calculate the three force constants from the two 
frequencies w; and w;. By analogy the following 
conditions were assumed 


de4/dcyx=0 and de,/dc;'=0. (8) 


Fic. 6. Relation between force constants of the carbon 
disulfide molecule (megadynes/cm). 
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Fic. 7. Relation between ratio (R) of atomic distance of 
C=O, and C=S. (R=rcuo/rc=s) and the angle force 


constant ‘‘c”” (megadynes/cm). 


The last two statements in Eq. (7) can be 
written 


(A1+A3) M3161 — (us1/ mor) 
+ 
— (1/mi*) (9) 


Differentiating Eq. (9) gives dcs/dc, and using 
the condition Eq. (8) yields 


(Ai +A3) M31 — +2 0. 


(10) 


Substitution of Eq. (10) into Eq. (9) gives the 
following explicit values of the four force con- 
stants: 


+ 
X (Ai As) } J+ 1") J, 

X As) } J+ [2 01) J, 

(1/91) (121/731) 
+(1/p21) (131/721) 
Co=[{ + 

X (Ai } 1°) J. 


(11) 


The carbonyl sulphide molecule (COS) was in- 
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vestigated by using the following frequencies: 


w1(859), w2(527), w3(2079 ;* 
rowo(1.151), 


The ratio R( =7rc-o0/rc=s = 0.7484) yields the force 
constant c2=0.0364 (megadyne/cm). This value 
is not as close to the average (0.0402) of the 
Ce-values of carbon dioxide (0.0571) and of 
carbon disulphide (0.0234) (megadyne/cm) as 
one might expect for the COS molecule. However 
an investigation of the relationship between R 
and c2 (Eq. (11)) shows that the average value of 
ce (0.0402) is not compatible with any real 


04 12 20 


Fic. 8. Relation between force constants of the carbonyl 
sulfide molecule. c: =CO-bond force constant and c;/=CS 
bond force constant (megadynes/cm). 


6H. Sponer, Molekiilspektren und ihre Anwendung auf 
chemische Probleme (Verlagsbuchhandlung, Julius Springer, 
Berlin, 1936). 

7W. Jevons, Band Spectra of Diatomic Molecules (Uni- 
versity Press, Cambridge, 1932). 

8H. A. Stuart, Molekil Struktur (Verlagsbuchhandlung, 
Julius Springer, Berlin, 1934). 

®L. Pauling, The Nature of the Chemical Bond (Cornell 
University Press, Ithaca, New York, 1940). 

10 G. Herzberg, Molecular Spectra and Molecular Structure 
I. Diatomic Molecules (Prentice-Hall Inc., New York, 
1934). II. Infrared and Raman Spectra of Diatomic Molecules 
(D. Van Nostrand Company, Inc., New York, 1945). 


| 


irbony! 
=CS 


ing auf 
ringer, 


s (Uni- 
ndlung, 
Cornell 
tructure 

York, 


olecules 
5). 
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value of the ratio R. The relation is shown in 
Fig. 7. The ratio R for the largest and smallest 
value of c. were found by differentiation to be 


(12) 
Substituting numerical values: 


(max.) ]= 1.1285; 
c2=0.0383 megadyne/cm. 


R{c2 (min.) ]= —1.1285; 
c2= —0.1774 megadyne/cm. 


Neither of these values of R is satisfactory for 
they would lead to impossible values of either 
the C=O or C=S distance in carbonyl sulphide. 
Hence cz calculated from Eq. (11) has been 
accepted as 0.0364 megadyne/cm. The other 
force constants for the COS molecule are shown 
in Fig. 8 and from Eq. (11) were found to be: 
=0.9461, cy’ =1.2049, and —0.1252. 

The two bond force constants (c; and c;’) show 
however an unexpectedly large change from the 
corresponding constants of = 1.5335) and 
CS, (0.7456). Hence, the conditions Eq. (8) are 
not satisfactory for this triatomic molecule. In 
the present case another procedure was adopted. 
It may be expected that the interaction constant 
cs for COS should not differ greatly from the 
mean of the corresponding constants for CO2 


TABLE IX. Force constants and interatomic distance. 
Badger’s rule.* 


Molecule Remarks 


CO(A 1) x d 
ed 
ed 


—0.1252/ 0.4 


COU 
CO(B 12) 


c4= —0.1252/ 0.4 
c4=0.0823 
¢4 =0.06019 
i 


@ See reference 11. 
See 7. 
we Hibben, The Raman Effect and Its Application (Reinhold 
Publishing Corporation, New York, 1939). 
reference 9. 
* See reference 1. 
! See reference 6. 
2 See reference 8. 
+L. O. Brockway, Rev. Mod. Phys. 8, 231 (1936). 


*T. Y. Wu, Vibrational: oe and Structure of Polyatomic Molecules 


(Prentice-Hall, Inc., New York, 1941). 
reference 10. 


16 


r(A) — 
Fic. 9. Badger’s rule for CO and CS. 


(cs=0.10522) and for (cs=0.05934). This 
mean value for COS would be 0.08228, yielding 
for COS the constant c;= 1.5593 and c;' =0.7243. 
Another method of estimating the interaction | 
constant (cs) for COS is to assume that the bond 
constant ¢; for CO (1.5335) remains the same in 
COS, yielding an interaction cy for COS equal to 
0.06707, and similarly that the force constant c,’ 
for CS» (0.7456) also remains the same in COS. 
This proposal yields cy=0.05333 for COS. The 
average of these values of cy for COS (0.06707 
and 0.05333) or 0.06019 permits the calculation 
of c; for COS (1.52159) and c,’ (0.74015). Either 
of these sets fits into Badger’s rule as can be 
seen from Fig. 9 and Table IX. In the case of the 
CO-bond a reasonably good straight line is 
obtained considering the fact that the inter- 
atomic distances are not known accurately in 
most cases. For the CS-bond only three cases are 
known, insufficient to draw a satisfactory straight 
line. But a line drawn parallel to the CO-line 
fits the CS-points well enough. The force con- 
stants c; (0.9461) and c;’ (1.2049) would not fit 
into this scheme at all. Hence the most reason- 
able values for the force constants of COS are 
C1 (1.5216), cx’ (0.7401), cz (0.0364), cs (0), cs’ (0), 
and c, (0.6019) megadynes/cm. 


1938). M. Badger, J. Chem. Phys. 2, 128 (1934); 3, 710 
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5 
value 100 
f the t oe 
1) as 
en R 
jue of CO 
real 
4 COL") 
12 14 || 
C1 
(mega- 
dynes 
r 
CO(COS) 2079 0.9461 1.0186 1.16 
CO(COS) 2079 1.5593 0.8624 1.16 cs=0.0823 
CO(COs) 2079 1.5216 0.8694 1.16 cs=0.06019 
2363 1.5335 0.8671 1.15 iad 
2168 1.8886 0.8090 1.128 j 
2182 1.9141 0.8054 1.118 b 
Cs(COs) 2079 1.2049 0.9616 1.56 
CS(COs) 2079 0.7248 1.1130 1.56 
CS(COS) 2079 0.7401 1.1055 1.56 
| CS(CSs) 1522 0.7456 1.1028 1.54 
Cs(E+) 1285 0.8416 1.0592 1.536 
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The lower excited energy levels of the molecule diphenyl are calculated by the method of 


antisymmetrized products of molecular orbitals. The wave functions for diphenyl are con- 
structed from linear combinations of benzene molecular orbitals and the perturbation of the 
original energy levels of benzene caused by the presence of the second benzene ring in diphenyl 
is computed. The theory predicts the splitting of the two non-degenerate singlet levels and the 
one degenerate singlet level of benzene into eight levels, four of which give allowed electronic 
transitions to the ground state and four of which do not. Experimentally three strong struc- 
tureless bands are observed, one between \=2200A and \=2800A, one at about \=1900A, 
and one below \=1700A. The calculated spectrum has allowed transitions at about 2400A, 


I. INTRODUCTION 


N a previous paper! the method of antisym- 

metrical molecular orbitals has been applied 
to a calculation of the lower excited electroriic 
energy levels of benzene. The results obtained 
from this calculation have been extensively com- 
pared to experiment; analysis of the observed 
spectrum? indicates good agreement between the 
theory and experiment. Since benzene is one of 
the simplest of aromatic compounds, it was 
decided to extend the usefulness of the results 
obtained in the benzene case to other more 
complicated molecules having benzene as a 
structural element. Sklar* has considered the case 
of the ultraviolet absorption of singly and of 
multiply substituted benzenes. Here, we con- 
sider the case of diphenyl which is a special case 


*A dissertation submitted to the Faculty of the 
Graduate School of Arts and Sciences of the Catholic 
University of America in partial fulfillment of the require- 
ments for the degree of Doctor of Philosophy. 

t National Bureau of Standards, "Washington, 


1M. Goeppert-Mayer and A. L. Sklar, J. Chem. Phys. 
6, 645 (1938). 2 Abbreviation GMS used in future references 
to this paper. 

2H. mang’ G. Nordheim, A. L. Sklar, E. Teller, J. 
Chem. Phys. 7, 207 (1939); F. K. Herzfeld, Proc. Am. 
Phil. Soc. 82, 359 (1940). 
3A. L. Sklar, J. Chem. Phys. 7, 984 (1939). 


1900A, 1500A, and 1400A, and thus is in good agreement with the experimental observations. 


of a substituted benzene in which one hydrogen 
atom is replaced by another benzene ring (minus 
one hydrogen atom). 

For the purposes of this calculation, it is 
assumed that the molecule is flat, i.e., the ben- 
zene rings are coplanar. A number of inves- 
tigators have discussed this point. Dhar has 
concluded from x-ray studies‘ of diphenyl that 
the molecule is flat and that any deviation from 
planarity, if it exists, is less than 0.1A. Dipole 
moment® investigations of para-dinitro and 
para-dichloro derivatives of diphenyl by Smythe 
and others have indicated that dipheny] is flat. 
H. de Laszlo® has obtained electron diffaction 
patterns of diiododiphenyl which are in agree- 
ment with the planar diphenyl structure. 

Pauling’ has discussed the possibility that the 
two benzene rings in diphenyl -are twisted with 
respect to each other about the longitudinal axis 
located in the mean plane of the molecule. 


4 J. Dhar, Indian J. Phys. 7, 43 (1932); J. Hengstenberg 
and H. Mark, Zeits. f. Krist. 70, 283 (1929). 

“<. F. Smythe, Ind. Eng. Chem. 23, 224 (1931); A. 
Weissberger, R. Sangewald, and G. C. Hampson, Trans. 
Faraday Soc. 30, 884 (1934). 

(1938) deLaszlo, Trans. Faraday Soc. 30, 884 and 892 
7L. Pauling, Nature of the Chemical Bond (Cornell 
University Press, Ithaca, New York (1940), p. 217 ff. 


SPECTRUM OF DIPHENYL 


Pauling’s argument may be summarized as 
follows. The conjugated double bond system in 
diphenyl causes the central C—C bond between 
the two benzene rings to have partial double 
bond character. As a consequence the molecule 
tends to remain coplanar. However, because of 
the small double-bond character of approxi- 
mately 20 percent, the force which tends to 
cause coplanarity is not great. The steric repul- 
sion between overlapping hydrogen atoms lo- 
cated in the adjoining benzene rings may cause 
the molecule to twist. Moreover, the molecule 
1,3,5 triphenyl benzene which would be expected 
to have steric interaction of the same order of 
magnitude as diphenyl has been reported to be 
non-planar. Thus, according to Pauling, the pos- 
sibility exists on theoretical grounds that the 
molecule is not planar. The extent to which the 
planarity assumption is valid will be indicated 
by the results of the computation, since con- 
siderable deviation of the diphenyl molecule 
from a flat structure would seriously influence 
the accuracy of the calculation. 


Il. SYMMETRY OF THE DIPHENYL MOLECULE 


The molecule diphenyl belongs to the sym- 
metry® point group Dy». The group Dy», is a 
dihedral group which has three mutually per- 
pendicular twofold axes of symmetry (Fig. 1): 
C27, C2”, and C,*?; and three mutually perpen- 
dicular planes of symmetry: o; the plane of the 
molecule, o,*, and o,”%, normal to the C,? and 
C,” axes, respectively. These 6 symmetry opera- 
tions together with the identity operation E and 
inversion operation I form the complete Dy 
group. The irreducible representations which 


TABLE I. Character table for group Do. 


E Co? I 


T,?(By,”) 
Ty?(Bo,?) 
T;?(Ai,?) 
ou?) 


op Sa 
1 

‘ 


a . — and G. Murphy, Rev. Mod. Phys. 8, 317 


Fic. 1. Symmetry axes and symmetry planes of diphenyl. 


specify the symmetry of the group D., are given 
in Table I. 

In Table I, the irreducible representations are 
represented by the [’’s in which the superscript 
D identifies the representation as belonging to 
diphenyl thus differentiating this set of I'’s from 
the I'’s used by Wilson® in discussing the benzene 
case. The A’s and B’s in parentheses are the 
equivalent symbols in the notation of Mulliken 
and Placzek." 

There are 66 degrees of freedom in the mole- 
cule which are distributed among the symmetry 
groups according to the equation, 
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where I? represents the reducible representation 
of the whole molecule. Of the 66 modes given by 
Eq. (1), only 60 are vibrational modes, since 6 
of the coordinates are required to give the 
rotation and translation of the molecule as a 
whole. The I’s corresponding to the “‘modes”’ 
belonging to rotation and translation of the 
molecule as a whole may be obtained by the 
following considerations. For the translation 
“‘modes”’ find the symmetry of the deformation 
pattern for the molecule when all atoms have a 
displacement in the x direction, corresponding to 
the x coordinate of translation for the molecule, 
and similarly for the y and z coordinates. In the 
case of the rotational coordinates, say rotation 
about the z axis, the deformation pattern is one 
in which all of the atoms have a displacement 
tangential to a circle drawn in the plane of the 
molecule about the center of symmetry of the 


%E. % Wilson, Jr., Phys. Rev. 45, 706 (1934). 


10 R. S. Mulliken, Phys. Rev. 43, 279 (1933). 
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molecule. Thus, Table II gives the irreducible 
representations to which belong the normal coor- 
dinates Q.7, Q,7, Q:7, Q:* where the 
superscript J and w represent translation in the 
direction of, and rotation about the given axes, 
respectively. Hence, the symmetry of the vibra- 
tional modes of diphenyl, ['y?, becomes 


Ty? 


+475? +1006? (2) 


Equation (2) tells us that there are 11 normal 
coordinates having the symmetry I”, 6 normal 
coordinates having the symmetry I”, etc. In 
contradistinction to benzene none of the fre- 
quencies is degenerate, since the characters are 
equal to +1 for all T’s. Whereas in benzene 


there are present 20 different vibrational fre- 


quencies, 10 of which are doubly degenerate, in 
diphenyl, the addition of another benzene ring 
has produced a new molecule of a lower order of 
symmetry such that there are no longer present 
degenerate frequencies. There are 60 different 
frequencies which must be computed. The group 
theory is of assistance here in that from (2) we 
can tell in advance the extent to which the 
secular determinant for the normal modes is 
factored. Thus, the 60th-order determinant 
factors into one 11th-order determinant, three 
different 10th-order determinants; two 6th-order 
determinants, one 4th-order determinant, and 
one 3rd-order determinant. 


Ill. THE DIPHENYL MOLECULAR ORBITALS 


In formulating the wave functions to be used 
to obtain the electronic energy levels of diphenyl 
one may set up functions which are appropriate 
products of the one-electron orbitals of benzene 
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given by Hiickel," as, 


(3) 


o1(v) = (601) exp?*4/6K (yp), 


In the above, K represents the 2pr wave 
function for the (1+)th carbon atom, where K 
takes on the values I, II, ---, VI while & has 
the corresponding values 0, 1, ---, 5 for the 
benzene case. The 6 carbon electrons which give 
benzene its aromatic properties' (GMS) are 
represented by small Greek letters such as », yu, 
and numerically assume the values 1, 2, ---, 6 
in benzene. Again / may assume the values /=0, 
+1, +2, 3 corresponding to different one-electron 
benzene orbitals. The coefficients o; are nor- 
malization coefficients. 

The starting point for setting up the diphenyl 
orbitals is the consideration that diphenyl con- 
sists essentially of two benzene rings. Thus, there 
are now twelve 27 electrons to consider and it 
is evident that if there were no interaction 
between rings the energy levels for diphenyl 
would be the same as those for benzene but 
they would now each be doubly degenerate. 
However, since there does exist an electronic 
interaction between rings the problem becomes 
equivalent to a perturbation calculation in 
which, in order to find the perturbation energy, 
it is convenient to average the interaction energy 
function over the unperturbed state of the mole- 
cule, thus necessitating a determination of the 
unperturbed wave function. This unperturbed 
function can be taken as a linear combination of 
the molecular benzene functions (for the two 
different rings) which belong to the same energy 
level. The linear combinations, in addition, 
should satisfy the condition that they have the 
same symmetry as diphenyl. That is, the trans- 
formation properties of the combination wave 
functions under the action of the symmetry 
operations of the point group D2, must be such 
that each one of the unperturbed functions may 


TABLE II. Symmetry of the translational and rotational 
nornfal coordinates. 


Q.” 


Q." 


Q.° 


T;? 


Q,* 
Ty? 


Q." 


" E. Hiickel, Zeits. f. Physik 70, 204 (1931). 
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be identified with some one of the eight irre- 
ducible representations given in Table I. 

There are a number of possible methods which 
may be used in setting up the unperturbed wave 
function. A straight-forward procedure, anal- 
agous to that used by Goeppert-Mayer and 
Sklar for benzene, may be used in which the 
linear combinations are made up of twelve-fold 
sums and differences of one-electron benzene 
orbitals belonging to the two different benzene 
rings, ring A and ring B, Fig. 2. The orbital for 
the state of the entire molecule consists of 
products of these combined one-electron orbitals. 

Another method consists in first setting up 
linear combinations of the benzene molecular 
functions which have the benzene symmetry 1, 
I's, P10, and T'y2 and which have been found to be 
the correct functions to use in the calculation of 
the benzene energy levels. These benzene func- 
tions are: 


vs, Wa; 


Yo = o-1(6), 
V1 = G0(1) b2(6), 
¥2= o-2(6), 
W3 = o-1(5) o-2(6), 
= 0(1)0(2)b-1(3) b1(5) b2(6) ; 
For diphenyl, we will use the functions 


\o 


which are analagous to W;, and Wz, instead of 
and , 

Taking linear combinations of the benzene 
functions given in Eqs. (4) and (6) as the unper- 
turbed functions of diphenyl makes possible the 
treatment of the diphenyl energy levels as per- 
turbed benzene energy levels, and the calculation 
then makes full use of the previous calculations 
for benzene. Furthermore, the secular deter- 
minant is not as unwieldy as in the first method, 
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for in this case it is automatically factored into 
quadratic terms from which may be obtained the 
splitting of each single benzene energy level into 
the two corresponding dipheny] levels. 

The correct diphenyl functions and the irre- 
ducible representations to which they belong 
are then: 

Wo? = Vos Von 

WP ViaVon} 

=24{ 
W_2? =273{ Wo, Von — Von} 

=274{ 
—VsaVon} 

Wy? Vas Von} 
Wig? = 24 { Van — Via Von} 


In the above, the superscript D identifies the 
diphenyl functions, while the function Wy, 
(J=1, --+, 4) is the benzene orbital defined in 
(4) and (6) for ring A, i.e., it is a function of 
electrons 1 to 6. Similarly Wig is the benzene 
orbital for ring B, i.e., it is a function of electrons 
7 to 12. The function Wo, is that for the ground 
state of benzene so that Wp»? is the molecular 
orbital for the ground state of diphenyl. 

A study of Egs. (7) will indicate the manner 
in which they have been set up and also their 
significance.. Thus, the function for the ground 
state of diphenyl, Yo”, is one in which benzene 
ring A is in the ground state and similarly for 
benzene ring B. A function like V,” represents a 
state of the molecule in which one-half of the 
time ring A is in the ground state while ring B 
is in the VW; state, and the other half of the time 
ring A is in the VW, state while ring B is in the 
ground state. The unperturbed molecular orbitals 
for diphenyl are thus written as products of two 
benzene orbitals, one for ring A and one for ring 
B. Taking linear combinations of the molecular 
orbitals belonging to the same energy, you ob- 
tain the unperturbed wave functions. It will 
be noted that these functions correspond to the 
lower energy levels, since only singly excited 
benzene rings are involved, that is, one of the 
rings always is in the ground state and no more 
than one electron is excited. 
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Fic. 3. Some symmetry elements for benzene. 


IV. SELECTION RULES 


In determining the symmetry of the wave 
functions it is only necessary to consider the 
effect of performing the various symmetry 
operations on the functions. Since benzene 
functions are involved we shall need to use the 
character table’ for the point group Ds, to 
which benzene belongs. Thus, the operation C,’ 
on W,? may be expressed symbolically as 


1 
(8) 


Reference to Fig. 2 shows that the operation C2? 
is equivalent to an interchange of ring A with 
ring B in the normal arrangement of the molecule 
followed by the rotation of each benzene ring by 
the amount z about its own center of symmetry, 
i.e., the operation C; (Fig. 3). Thus, (8) becomes: 


(9) 


Since the function Wo is symmetric to all sym- 
metry operations in the benzene group, we have 
= Voz and C2Vo4 = Vos. Furthermore, since 
Wia and Wig have the symmetry Is, reference 
to the character (x) table given by Wilson® shows 
that x(C2) = —1, so that CoVWi4=—WVi4 or (9) 
becomes 


or 


1 
{ (Won) (— Wis) 


C2? VP (10) 


Proceeding in a similar manner with the other 
symmetry operations one can tabulate the 
results for 


E C2 Cw I oh 
WP 1.,-1 1 -!1 
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where the numbers +1, are the characters for 
the transformations which result when the 
symmetry operations are performed on W,?. 
Comparison of the above with Table I shows that 
WV,” has the symmetry Is?. Proceeding in a 
similar fashion, it is possible to derive the sym- 
metry of the other W? functions as the following: 


VP cv, 
Wy cT,?, 


(11) 


Having established the symmetry of the func- 
tions given in Eq. (7) it is now possible to deter- 
mine the selection rules for the electronic 
transitions. 

The three components of the dipole moment of 
the molecule P,?, P,?, and P,” evidently have 
the symmetry properties—P,? cTs?, 
P,)cT;, since the dipole moments have the 
same symmetry as the translation coordinates 
given in Table II. 

The probability of an electronic transition 
between the ground state and any excited state, 
WV/?, is governed by the integral 


f WP*P dr 


where P¢? is any one of the components of the 
dipole moment P?. Since, if this integral is to be 
different from zero, it must remain invariant to 
any of the group operations, it follows that the 
equivalent group theoretical statement is that 
the integral must be totally symmetric; that is, 
the integral must transform like the irreducible 
representation T)?. As Wo? cI)? it is necessary 
that Po? cI,” to satisfy the condition of 
invariance. Thus, the allowed electronic transi- 


TABLE III. Symmetry of vibrations which make forbidden 
transitions allowed. 


P,P(rs?) P(r?) 
Electronic Symmetry 
3) T;? Allowed T,? 


tions for P,” are those in which the excited state 
has I's? symmetry, while for P,? the excited 
state must have I's? symmetry. P/ transitions 
are not allowed. Transitions involving the func- 
tions Wy”, and W,? result in radiation polarized 
in the x direction, while transitions involving the 
functions ¥,? and ¥;? result inradiation polarized 
in the y direction. 

However, the possibility exists of weak elec- 
tronic transitions to. states involving W_,1,? 
having the symmetry I? and [;? and also to 


states involving I's? and I's? but having different © 


polarization characteristics. All that would be 
required is that the transition occur to a vibra- 
tional level in the excited electronic state such 
that 
D 


where Ty” is the irreducible representation for 


The ground state function must be made anti- 
symmetrical to the interchange of any two 
electrons and is therefore taken as 


= 


where P is any one of the 12! possible permuta- 
tions of the electrons. Derivation of the nor- 
malization factor No» is given in Sections VII 
and VIII below. 

For the excited states, only singlet energy 
levels will be considered. We are mainly inter- 
ested in the perturbation of the benzene singlet 
levels since the triplet states of benzene give 
bands too weak to be observed. The singlet 
excited state spin function for electrons 1 to 6 is 


xsa=[1/v2 Ja(1)8(2)a(3)8(4) 
(14) 


The appropriate spin function for any one of 
the excited states, say V,”, is obtained by con- 
sidering each of the two terms which constitute 
V,” separately. Since in a term such as Wo4Viz 
the molecule is partially in the ground state and 
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x0? = a(1)B(2)ae(3)B(4) a(S) 8(6)a(7) 
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the vibrational level. It will be recalled that the 
To level of 5.0 ev in benzene is a forbidden 
transition of this kind. As indicated in Table III, 
which gives the symmetry of the vibrational 
level which is required to make possible the for- 
bidden transitions to the given electronic level, 
it is evident that these forbidden levels may 
appear weakly in the band spectrum since the 
vibrational modes of diphenyl include all of the 
(Eq. (2)). 


V. ANTISYMMETRIZATION OF THE ORBITALS 


Since the electrons in the various orbitals can 
be in one of two spin states, a or 8, and since no 
two electrons in the same orbital can have the 
same spin, it is necessary toinclude an appro- 
priate spin function as a°factor to the orbitals 
(7). For the ground state the spin function is 


(12) 


partially in the excited state V,, the spin function 
is correspondingly partially ground state and 
partially of the nature of (14), ie., the spin 
factor is xoaxse, Where xoa is analogous to (12) 
with the exception that it is a function of elec- 
trons 1 to 6 only. 

Accordingly, the totally antisymmetric orbital 
for is 


D 


~ (2-121N,)! 
X (15) 


where JN, again is a normalization factor deduced 
in Sections VII and VIII. Equation (15) can be 
shown to be equivalent to a product of coor- 
dinate function ¥,? and the singlet spin function 
for the entire diphenyl molecule. Thus, consider 
the term 


De (—1)?P{ VoaVisxoaxse}. 


It is evident that since the permutations of the 
12 electrons form a group it is possible to perform 
any additional permutation on all members of 
the group without changing the group. Suppose 
the permutation P is performed where P con- 
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sists in interchanging detionins 5 and 6. Then, 
Lp (—1)?P{ VoaVisxoaxsa} 
= —P 
=Lp(- 1)? HPC | 
= (16) 


inasmuch as, only Wo, and xoa can be affected 
by the permutation P of the electrons 5 and 6, 
while Wo, remains unchanged since electrons 
5 and 6 are in the same one-electron orbital. 
In the above 


Thus, 
2 Le 
= Lp (—1)?P{ — Pxoa)xsa} 


= Lp (17) 
Equation (15) then becomes 
D 
(2-121)! (—1)?P Wy xs”, (18) 


where 
xs? =xXsAXsp- 


The antisymmetric orbitals for the other 7 
excited functions are exactly analogous to (15). 


VI. THE SECULAR DETERMINANT 


In setting up the Hamiltonian for diphenyl an 
extension of the method employed in the benzene 
case will be used. In the latter method, the 
interaction energy between the six carbon atoms 
is divided into two parts consisting of (1), the 
mutual electronic interaction energy between the 
six 2pm electrons, and (2) the potential energy 
of the 2pz electrons in the field of the 6 carbon 
atoms stripped of their 2p electrons. The com- 
plete Hamiltonian function for the electronic 
energy of benzene, H, was taken as 


(19) 


where T is the kinetic energy operator for the 6 
electrons, i.e., 


1 
— V7, 


LONDON 


V2 being the Laplacian operator for the ith 
electron; Ho is the potential energy operator 
6 


which is equal to }> Hx, Hx being the potential 


K=1 
energy due to the Ath 2pz-stripped carbon 
atom; and where 


6 
ie 

the electronic interaction. The superscript e will 
be used to indicate electronic repulsion energy 
terms. 

For diphenyl, the kinetic energy term is ex- 
panded to 74+T7s:, since both rings A and B are 
involved ; the electronic interaction term becomes 


12 
Hp? 


corresponding to the mutual interaction between 
the twelve 2pz electrons; while the potential 
energy term becomes where Hoa 
differs from Ho, by an additional term designated 
as Hoa® corresponding to the fact that the 2pz 
electrons in ring B are now in the field due to ring 
A, i.e., Roa =HoatHoa® and 
To all of the above terms must be added a term 
Hp! which is the mutual interaction energy for 
the 12 carbon “‘ions,’’ i.e., the 2pz-less carbon 
atoms. Since the latter term is not a function of 
the 2pz electrons’ coordinates, it does not vary 
with energy level and hence is only of importance 
in determining the energy of the ground state; 
the analogous term was not considered in the 
benzene case for the reason that Goeppert-Mayer 
and Sklar were not directly concerned with the 
ground state energy. Similarly Hp! will not be 
considered here. Accordingly, the Hamiltonian 


for diphenyl, H?, becomes 
(20) 


The energy levels W, then are solutions of the 
secular determinant (23) where 


f i=j or ixj, (21) 


and 


(22) 


SPECTRUM OF DIPHENYL 


For 1=j, Ni:=1 since the ©?’s have been normalized as in (15). 


Hs—W| 


Ha—NaW 
Hu—NuW Hu—W 


Hii-W Hi3s—Ni3sW 
Hss—-W 


Hs2-W 
Hiiu—-W 


Terms not appearing in (23) are excluded since only interactions between states of the same 
symmetry do not vanish. This statement is immediately evident since any energy integral must be 
invariant with respect to any operation in the symmetry group D.; the H;; integrals transform like 
WW? XW/, which are totally symmetric when both ¥,? and ¥;” belong to the same I’. The functions 
,? are not quite mutually orthogonal, hence, the appearance of the orthogonality integrals Ni;;, 
which are evaluated in Sections VII and VIII. However, the @,’s are normalized so that the coef- 
ficient of the W’s appearing along the diagonal of the determinant is unity. The integrals H;; and 
Nj; are real as will be demonstrated later. 

Expansion of (23) leads to four quadratic equations each one of which determines two energy 
levels, given by the expression 


1 [ Hist H 
1—N;,7L 2 

Evidently this relatively simple solution of the eighth order determinant (23) results from our 
choice of the diphenyl functions as linear combinations of the benzene functions. Examination of 
(24) shows that if the interaction energy H;; is zero, and if the ith and jth functions are mutually 
orthogonal, i.e., N;;=0, the two levels given by (24) are Hi; and H;;. 


EVALUATION OF H;; AND Nj; 


It is convenient to separate the H;; integrals into two individual integrals, the electronic inter- 
action energy integral defined as 


f Pdr, (25) 


and the so-called nuclear interaction energy integral defined as 


Evidently the latter integral represents the interaction energy between the 2pz-less carbon atoms 
(or “‘nucleii’’) and the 2pz electrons. The superscript x will be used for ‘‘nuclear’’ * interaction energy 
integrals. 


2 It should be noted here that the “nuclear” interaction energy integrals also include kinetic energy terms; in fact, 
include everything but the electronic interaction energy terms. 
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The normalization integrals N;;, from which are determined the value of the normalization factors, 
Nj, and the orthogonality integrals N;; are similar in form to the H;; and H;; integrals, respectively, 
as may be readily seen by replacing the Hamiltonian operation H” with the operator /. 

The evaluation of all of these integrals may be considerably simplified as a result of symmetry | 
considerations. 

To exemplify the method involved we carry it out for i=j=1. Let Qi:(@) be one of the three 
integrals H;*, H4,:", or Ni;. Which of the three it is depends on whether the operator 6 is 


Hp’, or 1. 


Qn(8) = f (27) 


Then 


» Substituting from Eq. (15) for 6)? we get 


1 
6) (—1)P' [Wor *9 
Qu (8) (—1)?’ Via Vosxsaxon} 


XY (—1)?P{ (28) 


This integral may be simplified" in the usual manner and (28) may be written as 


. 
ON, Lp (-1)? | Via Vonxsaxon} *OP { (29) 
1 


Evidently Q1:(@)=1 for @=J from the definition of N;. Equation (29) may be written as four 
integrals, i.e., 


1 
Qu(8) = yy, (8) J, (30) 


where 


6Qu(0)=2 Le (—1)? 


= 2 (-1)? OP | 


aQu(@)=2 Dr (—1)” 


OP | ViaVonxsaxon} dr. 


It may be readily demonstrated that 
aQu(9) =aQu(A), 


Thus @ is unaltered, irrespective of whether it is Hp*, [74 +7e2+Hos4+Xoz ], or 1, if we interchange 
the symbols A and B. This interchange, which immediately changes .Q11(8) into @Q1:(@) and 4Q11(@) 
into .Q1:(8), is permitted since its only effect is to change the name of the integration variables. 

The general result follows immediately that 


aQii(9) = 
= -Qii(8) 


13 L. Pauling and E. B. Wilson, Introduction to Quantum Mechanics (McGraw-Hill Book Company, Inc., New York, 1935), 
first edition, p. 240. 


(32) 


|; +1, +2, +3, +4, (33) 


| 


SPECTRUM OF DIPHENYL 405 . 


for aQi:(6) and .Q;;(@) are exactly analogous to .Q1:(8) and .Q1:(8), respectively, it being possible to 
obtain the former integrals from the latter by replacing Wia (or Wiz) with Vis (or Vig where 
t=1, 2, 3, 4). The proof of (33) is exactly analogous to that given for (32). From the definition of 
the W_,? functions as given in Eq. (7) it follows that 


aQ_i-i(8) = 
= 
For the next step in the evaluation of Q1:(0) it is advantageous to consider the case where 6= Hp’. 
It is possible to split the integral Qi:(Hp*) or H1i:* into two parts, one part of which is due to the 


benzene electronic interaction in each ring, and the other part of which is due to electronic inter- 
action between rings. Thus let 


i=1, 2, 3,4. (34) 


12 6 12 6 12 
Ap = C/tut C/tutD (35) 


or 
Hp*= Ha +Has*. (36) 


Furthermore, the permutation operator P in Eq. (31) may be subdivided into four different types 
of permutations, i.e., 


(37) 


where P, represents permutations of electrons in the A ring only, and Ps, represents permutations 
of electrons in the B ring only. Pag is a symbol for the following class of permutations: The first 
permutation of the class matches up as many electrons as possible.’ The other permutations of Pas 
are all those permutations which start with the result of the first operation just described and in 
which at least one of the permuted electron pairs contains one electron in ring A and one electron 
in ring B. Under the latter definition, permutations also included in Pag are those for which an 
interchange between electrons in two different rings is concurrent with interchanges confined to 
electrons located in the same ring. P4.g represents permutations in ring A occurring simultaneously 
with permutations in ring B, but without any crossing over between rings. 
Substituting the values of Hp* and P from Eqs. (36) and (37) into .H1:° splits this integral into 12 
separate integrals. Fortunately, however, a considerable number of these are zero. If we consider 
first the integrals arising from H4*, we get for P,4 the integral, apart from the normalization factor, 


A 


The second integral in (38) is equal to unity, while the first integral, has already been evaluated by 
Goeppert-Mayer and Sklar (see GMS Eq. (15) ) and (38) is equal to }Eo* where 


Eo’ = (280+46;). (39) 


For (Ps, Ha*) and (P..s, Ha*) the integral is zero, permutations confined to ring B not being 
allowed as a consequence of the orthogonality of the one-electron orbital (including spin functions). 
This leaves the integral Paz, Ha* to be evaluated. 

Similarly, for Hz, the integral is zero for P4 and P4.z while for Pg we obtain 


3 (—1)? f = = Eis — (2¢0+3e1) ], (40) 


up means that the electron » occurs in this form: or 
Pip’ Vv 
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where E;, is given by GMS, Eq. (16), (see also Eq. (75)). Similarly the integral (Paz, Hg*) remains 
to be evaluated. 

In the third set of integrals involving Hs* only permutations of the type Pag are allowed. This 
follows from the fact that H4z* is the interaction between one electron in the A ring and one in the 
B ring, thus making impossible permutations confined to one ring, since all of the one-electron 
orbitals in the one ring are orthogonal to each other. Accordingly, combining this integral with the 
other two involving Paz which have been previously mentioned, there results the integral 


=4 p> (—1)?42 | { (41) 
AB 


where the lower case letter h* will be used for integrals involving terms containing only P4z and Hp’. 
Thus, we have shown that 


1 1 
ally; = + Eis ali’. (42) 


It is readily demonstrated that 


1 1 
=— =—— DY (—1)?48 { Ho’ (43) 
N, Pas 


Equation (43) follows from the fact that the benzene functions Wo and W; are mutually orthogonal. 
For (Pa, Ha*) and (Pz, H4*) the integral containing the B subscript functions is zero; for (P4, Hx) 
and (Pz, Hz*) the integral containing the A subscript functions is zero; similarly the (P4., H4°) and 

(P4.s, Hp*) integrals are zero. For (Pa, Haz’) the integral is zero since any permutation in the A 
ring will involve at least one A electron which is not functionally related to H4s*, thus making the 
integral vanish because of the orthogonality of the one-electron orbitals (including spin function). 
Similarly for (Ps, Haz*) and (Pa.s, Haz’) the integrals vanish. 

For the nuclear interaction energy aH1,", little advantage is gained in splitting off at this time 
the benzene nuclear terms. Also, for the normalization integrals, it is advantageous to retain the 
more general form. Hence, 

1 
Thus we have to compute the set of integrals and «N11; also phir’, and The 
two sets of three integrals are somewhat similar in form since the same molecular orbital functions 
appear in each integral. Each one of these integrals splits into four integrals if the expression given 
in Eq. (4) is substituted for ¥;. The four integrals which are obtained for .h1:° and which are of 


similar appearance for the other two in the set of integrals represented by the pre-subscript ‘‘a’’ are, 


ally, (1) = 


1 
Nu= (- 1)? 1P dr. (44) 
‘ 


AB 


z (—1)?48 | }dr. 


| 
i 
(45) 
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It is now possible to show that 


= (ashu*)*, (46) 
= 2Re(ayhu’+ ah’), (47) 


where Re signifies the real part of the expression in parenthesis. 

Equation (46) may be proved by taking the complex conjugate of asii°. This changes Poa*Pon* 
to Woaen and similarly Yospo2 becomes Poa*P2n*. Reference to Eq. (5) shows that the integral may 
be written as follows: 


so that 


(aghi*)* (—1)?48 | (48) 


However, since ¢;*=¢_, the only difference between Yosxoa and Woa*xoa is that electrons 3 and 5; 
4 and 6 are mutually interchanged. Since the integral (48) would not be affected by such a change 
in variable designation, Yoa* may be replaced by oa and vice versa, so that Eq. (48) becomes identical 
with A similar demonstration suffices to show the equality of and ashi. 

For o/1,° we similarly get four integrals: 


AB 


bolt’ =} (—1)?48 | dr, 
AB 


AB 


bli? =} p> (—1)P48 f adr. 
AB 


It may be shown that 


/ 


= = (dshu*)*, (50) 
= 2Re( + wh’). (51) 


Equation (50) follows from considerations similar to those given in connection with (46) with the 
exception that in order to convert Yoa to Woa* in the left-hand parenthesis in the integrand of, say, 
bst\,°*, interchange of variable designations y=3 with »=5 and 4 with 6 in the entire integral also 
changes W241. However, it is always possible to regain the original form of Yea by performing the 
inverse double interchange in the parenthesis to the right of P42, since this is an allowed permutation 
in Pag, which, furthermore, has no effect on the value of the integral. Similarly, it is possible to 
convert Wos* to Por. The details of this proof are given in Appendix I. 

We are now in a position to summarize some of the results obtained in the above. Evidently, 
we have shown that 


so that 


1 4 
=—(Eo’+ E1s*) +—Re(arhu’ + + wen’), (S2) 
N. 1 N 1 


and that this integral is real. Furthermore, for Hi:" we have 


Ay" + + mH"), 
1 ‘ 


and for Q1:(1), 4 


e 
n 
1) 
3) 
*) 
id (49) 
1e 
ne 
1e 

4) 
ne 
ns 
en 
of 
e, 
5) 

(53) 
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Proceeding in a similar fashion it is possible to obtain the following set of equations for the elec- 
tronic interaction energy terms: 


1 : 4 
Ay? =—( + +—Re[ (ahr? + + + deh’) |, 
N. 1 N 1 


1 4 
=——( Eo’ + Eis’) +—Re[ + — + beh’) J, 
Nu 


1 4 
Fo9° = —( Eo’ + Exs*) +—Re[ — + — J, 
N. 2 N. 2 


1 4 


1 ‘ (55) 
Hi33° = —( Eo’ + +—Re[ (arhgs° + + (drh33°+ b2hss°) |, 

N3 N3 


e 1 4: 
H_3_3=—(Eo’+ Ess*) + — + J,” 
N_3 N_3 
1 4 
Hae =—( Eo’ + +—Re[ (arhgs° — aghgs°) + — bohss°) |, 
N. 4 N. 4 : 


e 1 4 
H_4_4=—(Eo’ + Ess*) — azhss°) — — |, 
N_4 N_4 


where Ey’, Eis*, E2s*, E3s* are benzene energy terms (see Eq. (75) below). 

Evidently the appearance of the minus sign in H_,_;° and the relationship of the latter to Hi,° 
is in part a consequence of Eq. (34). The minus signs in H22° follow from the relationship between 
WV, and VW, given in Eq. (4). The four equations in which 1= +3, +4 follow in a similar fashion from 
Eqs. (4) and (34). All of these four equations have the common term E35°, due to the degeneracy of 
WV; and W,. In obtaining the term £35, an integral of the form 


B 
is involved. Applying Eq. (6), (56) becomes 


B 


+ f f | (57) 
The first and last integrals are equal to E3s°, while the second and third terms lead to integrals 
and 


Although H is invariant these integrals are not invariant to a rotation of say 27/6 about the sym- 
metry axis of benzene, and hence must be equal to zero. A rotation such as 2x/6=w multiplies the 
first integral by the term e~“e~‘*e‘*e2 =e and, the second integral by e~e?*e?#e—iw = erie, 
Thus, for the electronic interaction energy only eight separate terms need be evaluated for the 
H;;¢ values. Analogous terms containing the same orbital functions in the integrals are required for 


4 
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the H;," and the H;,(/). For each of the latter a set of eight equations similar to (55) may be con- 
structed by analogy from Eqs. (52)—(54). For reference purposes, the eight individual terms required 
to determine the Q;;(6) are written out in Appendix II in terms of the one-electron orbitals and their 
spin functions. 

Next, we consider the evaluation of the H;; terms when 1#j. Four sets of integrals need to be 
determined, i.e., for the combinations, (2, 7)=(1, 3); (2, 4); (—1, —3); (—2, —4). The method is 
analogous to that used for the Q;; terms. For illustrative purposes, we may consider the case of Q13. 

As before, introducing the general operator 6, it is possible to show that 


2 
Q13(4) = + 6Q13(9) J. (58) 


Evidently the reason for this is exactly that given in connection with Eq. (33), i.e., CoWi=— V1, 
and — 3. Thus, the general result that 


2 
is also true. Furthermore, 
aQ-i—;(8) = aQi;(9) 
6Q_:-;(8) = 


Considering specifically ,H13°, it splits into four terms, i.e., 


lia =(1, 3); (2, 4). 


1 


Dp (—1)? | { dr, 


Dip (—1)P { H jar. 


Utilizing Eqs. (46) and (47) in say a:/113° we see, as discussed after Eq. (43), that the (Pa, Ha*) 
integral is zero, since y; is orthogonal to y3. For (Ps, Hz°) the integral is zero since it leads to integrals 
which are not invariant to a rotation of w about the C, axis although Hp° is invariant (compare 
discussion in connection with (57)). For the same reason (P4, Hz°) is zero, while (Ps, H4°) is zero, 
because of the orthogonality of ¥; and 3. The only integrals which remain are those due to 
(Pas, Hp‘), so that becomes Similarly 13° = a3f113°=a3h13°, = agh3°. 
(See text following Eq. (41) for definition of h*.) It is also possible to show (compare derivations of 
(46)) that: 


ayhy3° = (aghy3°)*, aghy3° = (agh3°)*. (63) 


For »H13*, by an exactly analogous procedure, it turns out that 


(64) 


~ 


een 
rom 
y of 
(62) 
(57) 
srals 
the | 
the 
1 for 
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dr, 


=} 


3° = 


4 
V3 


13" + + 13" + (67) 


Also 


Ny 


and 


4 


Using (60) and the relationship between V; and W2, V3; and WY, we get in addition to (66) 


4 
Re[_(ashg° + aghis*) — + J, 


4 
Hoge = — + — |, 


(N24)? 


4 


Similar equations result for 0= Hp" and 6=1 (compare (67) and (68)). 
The four q;;(@) integrals are written out in terms of the one-electron orbitals in Appendix I], 
Eq. (2A). The lower case q is used to represent integrals involving the operator @ and its corresponding 
permutation operator P(6). 


THE NORMALIZATION INTEGRALS AND Ni; 


We next undertake the computation of the integrals N;:;, Ni; or Qi:(1) and Q;;(1), respectively. 
Since N;;=1, the N;; integrals really define the normalization factors N;. The procedure involved 
will be illustrated in the calculation for Nj. 

From Eq. (54) it will be seen that we need to evaluate terms such as a;N;. Equations (1A) in 
Appendix II are helpful here since, if in a191:(8) we substitute for the operator the value @= / and for 
P(@) the value P, we get an explicit formulation of the integral a1 N41. 

For the identity permutation in which the electrons remain in their originally sndaned orbitals, 
the integral a:N3; is equal to 3. The only other permutations which are not zero are those involving 
an interchange between rings. These result in integrals differing from zero, which represent the 
extent of overlapping between an one-electron orbital in one ring with another one-electron orbital 
belonging to the other ring. The diphenyl overlap integrals Az,tg’ are defined as follows: 


Their origin can be seen in the specific example of the allowed permutation which interchanges 
electrons (1) and (7) in a;N,;. This results in a contribution to the total value of a;Ni; in the amount of 


1 . 1 


410 
where 
(65) 
Thus, 
|| 
e 
—1/V_3 
( 
1 
| 
i 


SPECTRUM OF DIPHENYL 411 


A considerable number of A integrals arise as a result of the permutations. Hence, it will be con- 
venient to obtain a general expression for Ai4tg’ in order to obtain the numerical values of the 


various A’s involved. Substituting from Eq. (3) for ¢i4 and ¢1’ there is obtained 


fexp (—2twls) +exp (—4iwls) +exp —3l,)) 


+exp (tw(5lp’ —314)) |sas+Lexp —214)) +exp —414)) 

+ [exp (iw(le’ —414)) +exp (tw(Slp’ —214)) |sss+Lexp (—iwls) (—Siwls) 
+exp (tw(2lp’ —3l4)) +exp (iw(4l9’ —314)) +exp 

+exp (2iw(ls’ +exp +exp —L,)) 


where S47, S4s, °° 
a 2px function in the other. Thus 


+exp (Siw(le’ (72) 


*, Ses, *** are the overlap integrals between a 2pm function located in one ring and 


f IV(v) VII (73) 


The overlap integrals s47, 54s, are functions 
of the internuclear distances between the carbon 
atoms and may be evaluated in closed form.” 
The value of the integrals decreases rapidly as 
the distance increases, so that the overlap 
integrals for 2pm functions located at extreme 
ends of the molecule may be neglected. The 
multiplicity of terms which occur with each s is a 
result of the fact that the same internuclear 
distance is associated with more than one pair of 
atoms. In Table IV the internuclear distances, 
the equivalent atom pairs, and the overlap 
integrals are listed. The atoms or 2pm functions 


The average value, A, is taken as the average 
of the absolute values of the figures given in 
Table V with each possible /4lg’ combination 
given a weight of one, ie., A=0.039 while 
A?=1.52X10-*. From the table it will be ob- 
served that all A’s in which /4= +1 are negative 
in sign. This is a result of the relative phases 
between the $41, and the ¢),° orbitals. Thus 
Asi,ig’ OF Arprti, are to be replaced by —A 
while all other A; ,1,’’s are to be replaced by +A. 

Using Eqs. (1A), (2A) and (3A) in Appendix 
II as outlined in the above, the results given in 
Table VI below are obtained. 


| are numbered as shown in Fig. 2 (normal ; 
in arrangement) ° ( In the derivation of the normalization coef- 

Utilizing the overlap integrals of Table IV we ficient, lg iawn involving powers of A 
\) in may compute the A’s. The results are compiled 8teater than A — neglected, because of the 
1 for in Table V. Since all of the A’s are real by direct Very small numerical magnitude of such terms. 


tals, 
ving 
the 
bital 


(70) 


inges 
int of 


(71) 


calculation, it is readily shown that 
Atalp’= A-l4_lp’= Alg'ta. (74) 


In evaluating the normalization and other 
integrals, a large number of terms of the order 
of magnitude A?, similar to Eq. (71), and also A 
terms appear. As indicated in Table V, the 
absolute value of A varies little, varying between 
0.035 and 0.044. It will therefore be convenient 
and also greatly simplify the complexities of the 
computation, without causing appreciable error, 
if average values for all A’s and A®’s are used. 


*F, Hund, Handbuch der Physik (Verlagsbuchhandlung 
Julius Springer, Berlin, 1933) 24, No. 1, p. 643. 


IX. THE ELECTRONIC INTERACTION ENERGY 


In deriving the electronic interaction energy 
terms it will be instructive to consider in some 
detail the electronic interaction energy for the 
ground state of diphenyl, Eop*. From Eq. (5A), 
we need to calculate 4o(Hp*). Using (3A) we can 
establish the types of integrals involved. Thus, 
considering terms arising from electrons (1) and 
(2) we get a benzene Coulomb integral yoo mul- 
tiplied by A? which results from an interchange 
of an electron in ring B with one of the A elec- 
trons, numbers 3 to 6. The complete term is 
—12yo0A? since one interchange changes the 


65) 1 
(66) 
(67) 
(68) 
(69) 
ding 
. 
= 
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TABLE IV. Internuclear distances and overlap integrals 


for diphenyl. 
Internuclear Equivalent Overlap 
distance A distances integral 
1.48> 0.225¢ 
2.49 d57, ds 12 0.033 
0.014 
dis 3.75 ds 12 0.002 
49 3.76 ds 1, 0.002 
8 4.26 diz, do, ds 10, ds 11, do 12 0.0006 
d3 1 4.90 dz 12, des, ds9 
dis 5.10 112, @3 10, ds 10 
5.65 dou 
do 11 6.15 deg 
dig 6.46 dius, d210, 10 
di 10 7.04 


8 dq =distance between atoms 4 and 7. 

bL. Pauling, Nature of the Chemical Bond (Cornell University Press, 
Ithaca, 1940), p. 218. 

© $47 =0.225, etc. 


sign and since the three a- and §-orbitals in ring 
B interchanging with the two a- and £-orbitals 
in ring A, respectively (excluding electrons 1 
and 2), make possible 12 different permutations. 
If we carry out a similar process for the re- 
mainder of the electron pairs involved in the 
integral, we obtain as the total contribution the 
term —2E*-12A*. This result could readily have 
been foreseen since each one of the rings taken 
by itself gives the energy Eo’ and the 12Paz 
permutations introduce the factor — 12A?. Similar 
results involving terms like Eis* are obtained for 
the excited levels. Hence, it will be convenient 
to list here the value of such terms. Thus: 


Eo’ = 671 — 4601 — 261-1 1T1(A1g), 


E,s° = 6_12 — 612 — (75) 
Egs° = b12 


where 
EM? = + 391+ — 3601 — 502 — 61-1. 


Equations (75) are obtained from Eqs. (15), (16), 
and (18) in GMS by deleting the nuclear inter- 
action energy terms ¢€, €1, and €: from Eo, E™, 
E,s, Ess. 

In the above the integrals yj, are Coulomb 
integrals defined by 


e 


while 6, are exchange integrals, i.e., 


e 


The terms & and 7 are defined thus: 


t= f f 
(78) 


For the diphenyl energy terms we obtain 
similar Coulomb energy integrals involving one 
orbital in the A ring and the other in the B ring, 
i.e., 


e 


also diphenyl exchange integrals, 


e 
|*drdr,. (80) 


Other more complicated diphenyl integrals 
occur, the general form of which may be repre- 
sented as 


e 
(81) 


The /’s in the above may have either subscript 
A or B. 

Returning now to ho(Hp’), for electrons 1 and 
2, we also obtain the terms 


~2[040n, 0404 ]A—2[041p, 0404 
—2[04—13, 0404 JA, (82) 


TABLE V. Value of Ar,7,° for various values of 14 and /,’. 


lalp’ Al alp’ 
0408 0.041 
0.042 
0.038 
+1402 —0.041 
1lals, —la—ip —0.038 
14—1p, —lale —0.044 
1423, —0.035 
14—2p, —1a2e —0.041 
+2408 0.038 
2alz, —2a—1B 0.035 
2a—1p, —2alp 0.041 
2a2p, —2a—2B 0.035 
24—2p, —2a2e 0.042 


‘a 
| 
q 
q 


7) 
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ain 
ng, 


79) 


SPECTRUM OF DIPHENYL 


TaBLeE VI. Value of the normalization coefficients N;, and integrals Nj;. 


Ni, Nig Functional equivalent* Value 
No 1-18A> 
Nu 4[ (a, — (Nut 1-204? 
Ne 4[ (a, Niu — + 1-184? 
N-2 4[ (a, Nu— — 1-18A? 
4[ (a, Ng3-+ + 33+ 0253) 1-164? 
N-s (a; a2N33) — (01. 02. N33) 1-20A? 
Ns 4[ (a; Nas — a2N 33) + (01.33 — b2 1-184? 
Nu 4[ (a; N33— a2N33) — (01. N33 — 02.33) 1-184? 
Nis (a, Ni3+ a2Ni3) + (0: Nist+ is) 


N_1-3 4(N_1N_s)~*L (a1 Nis+ a2.Mi3) — (0:Nis+ 
(a,Ni3— a2Ni3) + Nis— 


4(N_2N_4)4 (a: Nis— a2 is) — Nis— 


Nu 


—2A*(1-20A2)- 
—6A*(1-18A2)-1 
—6A*(1-18A%)-1 


*4aiNi1 =4a1N 33 = 1 —18A2, 
4a2N 11 =4a2N 33 =0, 

=401N 33 =A?, 
4b2N 11 =4b2N 33 =A2, 


4a1N13 = —3A%, 
4a2N 13 =3A2, 
4b1N13 =A’, 

4b2N 13 


which result from interchanging electrons 1 or 2 
with 7 or 8, etc. In addition, it is possible to 
perform another interchange say, in [040z, 0404 |, 
thus obtaining the term 


[010z, 0402 JA? = 60402A?. 


The diphenyl exchange integrals are very small, 
for example, 6040, =0.018 ev. Terms containing 
a diphenyl exchange integral multiplied by A? 
are therefore negligible. 

The hybrid diphenyl integrals given in (82) 
are those for which three of the orbitals are in 
one ring and the fourth is in the second ring. 
However, one of the two electrons is in two 
identical orbitals, i.e., 0404, so that these terms 
might be called mixed Coulomb integrals. For 
the electron pair 1 and 3, in addition to getting 
mixed Coulomb terms such as [0404, 1402 JA we 
obtain mixed exchange terms such as —[040z, 
1404 ]A. The mixed exchange integrals turn out 
to be considerably smaller than the mixed 
Coulomb integrals, i.e., the former range in 
absolute value between 0.05 to 0.10 ev, while the 
latter vary little in absolute value ranging from 


ho(Hp*) =4(yoson +4 (1—12.5A2) 


With these approximations we obtain for ho(Hp°), 


(260408 + 860412 —360MA—24E)°A?. 


Similarly, for the excited electronic energy terms, a;/;;° (Eq. (1A)) we get 


0.26 to 0.30 with an average of 0.27 ev (average 
for six different integrals). All of the mixed 
Coulomb terms contribute negative energy terms 
to ho(Hp*) while all of the mixed exchange terms 
are positive contributors. The number of dif- 
ferent mixed energy integrals which appear is 
considerable yet they add little to the value of 
the energy levels. Hence, in order to simplify 
the calculation, we replace all of the mixed 
Coulomb integrals by the average value MM, 
without causing any appreciable error, while 
we neglect the mixed exchange integrals. It is 
estimated that the error in the total value of 
these mixed energy terms caused by neglecting 
the exchange integrals is less than 10 percent of 
this total value. Since the latter only affects the 
location of the final energy levels by possibly 
0.2 or 0.3 ev the error caused by dropping the 
mixed exchange integrals is negligible. In car- 
rying out this simplification it is well to keep 
in mind that all mixed energy integrals are 
positive with the exception of those in which one 
or three of the four orbitals belong to ring A, 
in which case they are negative. 


(83) 


= (1 —12.5A?) —12(Eo’+ s*) A? — 12(£—2n) A? 


= (260402+ 7 60418 + 60422 + 351418 +3614 56-1422) —360MA, 


413. 
J 
N_o-4 | 
| 
80) 
rals 
pre- 
(81) 
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and 
(82) 
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+2[1414, A’, 
=2[—1424, ](1—12.5A?) —[—14—12, 2524 ] 
—5{4[14—24, 0208 ]+6[14—2a4, 1212 A’, 
=2[14—24, ](1—12.5A2) —[14—12, 22—2a | (84) 
—5{4[14—24, 0202 1+6[14—24, ]} A’, 
Aashss° = (4yos0n + (1 —12.5A?) — 12(Eo°+ A? 
50422 + 351418 +3614 —360MA, 
—2222 ]+2[1414, ]} A’, 
= 2[1424, ](1—12.5A*) —[1 412, 2824 ]—5{4[1424, 0205 ]+6[1424, ]} A’, 
—2e—12 ](1—12.5A?) —[14—12, —2824]—5{4[1424, 0202 ]+6[1424, 1512 |} A’. 
In the above, reference to Eqs. (1A) shows that aih33° and a2h33* may be obtained from ai/1,° and 
respectively, by replacing J=2g by 1=—2z, while and may be obtained from 


and by replacing —24 by /=+2,4 and /=2, by /=—2z. 
In deriving Eqs. (84) use was made of the relationship 


(85) 


: This may readily be demonstrated by a method analagous to that used in deriving Eq. (32). 
We also determine the a,/;;*, from Eqs. (2A). These follow: 


4 4ayhys° = 2((0404, 28—22 |+2[14la, ])(1—12.5A*) —[04—22, 2204 ]—[14—2e, | 
2p—14]—60MA+ 21 yu +17 61-1 
— 100402 — A’, 


1214]+[14—12, 12—14]+60MA+ (Sy00+ 2871+ 15 y+ 6702+ 6712 
— 14601 — 3602 — 761-1 +10 3070412 +5 0428+ 20y1418+ 1071428) A’, (86) 
20 —14—24, 18 ](1—12.54%) —[-14—12, 23-24] +10MA 


= 2[ 1424, ](1—12.5A*) —[14—12, 2824 ]+10MA 

| +5{3[1ala, 22—12]+2[0404, 22—12]—3[1424, 1212 ]—2[1424, 0208 ]} A®. 
may be obtained from by replacing —14 by J/=14 and /=—2,4 by 1/=2,. 


X. NUCLEAR INTERACTION ENERGY TERMS 


The nuclear interaction energy terms a,/;;" (Eqs. 1A) may be expressed in terms of the following 
integrals: 


a= + Hoa, | ola, subscript A or B. (87) 


This is a benzene integral which gives the energy of the /th orbital. In diphenyl, the energy of an 
electron in the /th orbital of, say, ring A, is affected because it is in the potential field of the ‘‘nucleii”’ 
of ring B. This results in the integral 
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The total energy of an electron in the /th orbital in diphenyl is therefore given by 
(89) 


In addition, integrals occur in which an electron is partially in both an A and B orbital. Thus, 
the energy of an electron partially in both rings due to the nucleii of both rings is given by 


f $14*(v)[Ta(v) + Hon (v) +Hoa®(v) (90) 

where in a manner similar to (89) 
ctain'= f + Hone) (91) 
Mata’ Hoa (92) 


Also an integral of the form 
Alla’ = f (93) 


is obtained. This integral evidently represents the potential energy of an electron belonging to two 
different orbitals of the same ring due to the action of the nucleii in the second ring. In all of the 
above integrals, the Hamiltonian operators operate on the function which follows the operator. 

From symmetry considerations we can reduce the number of integrals which need to be evaluated. 
The Hamiltonian operators are totally symmetric so that if we perform the rotation C2* on say (92) 
we get 


Hox") or Alalp’=AColz, Cola’. (94) 


A similar symmetry equation holds for both (91) and (93). In the ultimate evaluation of the integrals 
appearing in the excited state terms, only the real parts need be considered; hence, wherever con- 
sidered desirable, either the integral or its complex conjugate was used. 

The type of interchanges allowed are similar to those obtained in the electronic interaction energy 
terms, i.e., terms involving both A and A? appear. Thus, for the ground state the nuclear energy 
is given by 


NoE op" = (4w9+ 8) (1 = 15?) + ( — 4w040z3 — -ip)A — 480414 (95) 


The term (4w9+8w:) arises from the fact that for the ground state each one of the benzene rings 
has two electrons in the /=0 orbital, two in the /=1 orbital, and two in the /= —1 orbital. In addition, 
for any one of the energy terms w; which is obtained when electron » is in orbital /, 15 permutations 
are possible between the remaining 11 electrons. Thus, suppose » is in orbital /4 and has the spin a. 
Then there will remain 2 electrons in ring A having spin a which can be interchanged with 3 electrons 
in ring B having spin a, thus giving rise to 6 interchanges; also there will remain 3 electrons in ring 
A having spin 6 which can be interchanged with three B electrons of spin 8, thus causing 9 more 
interchanges, or a total of 15 interchanges. Each one of the interchanges gives rise to the factor — A’. 
Since No>=1—18A? (Table VI), we may write: 


(4wo+ (1 — 15.A?) = (4w+ 801) No+3(4w9+ 8a) A’. (96) 
For the excited states, the terms are: 
= (4wo+ + we) (1 — 15.A?) + ( — — 1400413 — 20422 + Alp + 
+ A+ (—42d0414 + 60424 — 6A14 -14 — 6A1424 + 614-24) A’, 
= 3(A14-14 +214 -24-+A24 -24) 
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= (w242p — wialp) A— 10A14-24A", 
= (w24-22 — wia—1p) A— 1014-24 A", 
= (4wo+ 71 + we) (1 — 15?) + ( — 1400413 — 200422 + sip 
+2w142p+ 2w-1428) A+( ~4200414+ 600424 — 6A14-14+6A1424 — 6A14 -24) A’, 
= 3(A14—-14+2A1424 +A24-24) A’, 
= (w242p — wiaip) A— 10A1424A?, 


(97) 


= (w24 —wis-ip)A 10A14244”. 


For the excited states, one of the electrons is in orbital /= +2, thus giving rise to the term 4wo+7w1 
From Table VI, we have 


(49+ (1 — 15.A*) = + 7001 + Ni + (400+ + w2) A?. (98) 
Similarly we may obtain equations for the ah;;". Thus: 
= 24 -24(1 —15.A?) — 2(wo42p — wis — wi A— (119+ — 6d0424) A’, 
Aaghy3" = —d14-14(1—15.A*) 
+ (1109+ 160 + +3A1424+3A14-24) (99) 
Adjhy3" = A+ +401 — 51424 — 514-24) 
Abghy3" = (—wia-1p+ A+ (2wo+4e1 — 5A1424 — 5A14 -24) A”. 


XI. EVALUATION OF THE ENERGY LEVELS 


The numerical values of the energy levels are given by W;;*—Eop where W;;* is given by Eq. 
(24). From the latter it is seen that we need expressions for }(H;:+H;;)—Eop, Hii, Hij, and 
(Hii— 

Using Table VI for the normalization coefficients, Eqs. (53) and (55) for determining the manner 
in which the individual «,q;;(6) are compounded for both the nuclear and electronic terms, and using 
Eqs. (83), (84), (95), and (97), Table VII may be constructed for the values of H;;, Hi:i—Hj;, and 
4(H;:+H;;) —E,?. In this table, the following abbreviations are used: 


(88) Ao=4A0+8A1, 
(93) A! = —42)0414+6d0424 
(91) = — 1400415 — 200428 + + 6014-18 + + 20-1422, 
Qo’ = —4w0408 — 16w041p+ 
(79) G=4y040p+ 14y041p +2 +12y1418 +4 71422, (100) 


G =4y040p+ 
G—Go=2(yoa2z — +4(y142B— ; 
(80) D=260408+7 60412 + 60423 +361418 +3614 61428 + 56-1428, 
Do = 260403 + 860413 +46141g +4614 -12, 


D-Dy= 60423 — 6142p — + 614-1. 


Numbers in parentheses on the left hand sides of the above equations are references to the equa- 


tions in the text which define the terms appearing in Eqs. (100). 
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All of the integrals listed in the table with the 
exception of some of the bracket integrals such as 
[0z—12, 120. ], [142a, —2z12], and [142a, 
22—1s8] when evaluated are found to be real. 
For the complex integrals only the real parts 
need be taken in accordance with Eq. (55), etc. 

Similarly Table VIII may be constructed for 
the H;; terms from Eqs. (66), (67), (69), (86), 
and (99). 

The numerical value of the various integrals 
involved may be obtained from the tables given 
by Sklar and Lyddane'*® for the Coulomb, 
exchange, and hybrid Coulomb exchange repul- 
sion energies resulting from the repulsion of two 
2pm electrons. In addition the Coulomb and 
exchange attraction energies of a 2pm electron 
for a carbon ion are given. In order to use these 
tables, a screening constant of z=3.18'" for the 
carbon atom was used. It was then possible to 
plot the value of the above integrals as a function 
of the internuclear distance. Using the distances 
given in Table IV the various energy integral 
terms were evaluated. The numerical value of 
the Coulomb energy integrals given by Sklar 
and Lyddane extend to internuclear distances 
of 3.33A only. However, at this value of distance 
the Coulomb energy has practically become e?/r, 
so that this law was used to extrapolate the 
Coulomb integrals to distances as large as 7A. 

In this connection the Coulomb energy terms 
given in GMS require some revision as a con- 
sequence of an error!® in sign in the original 
tables of Bartlett’? from which the Coulomb 
integrals were determined. The new values as 
determined from Table 1 of Sklar and Lyddane’s 


paper are: 


A\= f f —I*(v) 1 1?(u)dz,d7, =9.05 volts, 
Tons 
2 
A,= ff —]?(v) = 5.68 volts, 
Tos 


e 
f IV?(u)d7,dr, = 4.98 volts. 


'® A. L. Sklar and R. H. Lyddane, J. Chem. Phys. 7, 
374 (1939). 

17 C, Zener, Phys. Rev. 36, 51 (1930). 

18W. H. Furry and J. H. Bartlett, Jr., Phys. Rev. 44, 
210 (1938), footnote 11, p. 2 

18 J. H. Bartlett, Jr., P 


11. 
ys. Rev. 37, 507 (1931). 


The integrals evaluated by Sklar and Lyddane 
are the only ones which will occur if interactions 
between electrons which are not on neighboring 
atoms are neglected. This assumption was used 
in the GMS calculation with the exception that 
Coulomb interaction integrals arising from non- 
neighbors were retained. In the present work it 
was found necessary to include the energy terms 
arising from non-neighbors, since their neglect 
caused appreciable error. As a consequence ex- 
pansion of the energy integral terms given in 
Tables VII and VIII resulted in general in- 
tegrals of the form 


(0b; cd) = f f a(v)b(v) 
X 


where a(v) is the real lpr function for electron ». 
The integral (ab;cd) was approximated by 
(reference 3, p. 990) 


(101) 


(ab; cd)~S arSea f K 
X 


where Sq, Sea are overlap integrals, Eq. (73), and 
the function K.?(v) is a 2pm function located 
midway between the original a and 6 functions 
and similarly for K.a. The integral (ab; cd) then 
becomes equal to a product of two overlap 
integrals by a Coulomb integral which is evalu- 
ated for an internuclear distance corresponding 
to the distance from the midpoint of a and b to 
the midpoint of c and d. Similarly, an integral of 
the form (ab;cc) was approximated by the 


(102) 


‘ product of the overlap s,s, with a Coulomb in- 


tegral between an electron on the 27 function 
c and a 2pm function located midway between 
a and b, This same type of approximation was 
used in evaluating ionic energy terms which had 
not been tabulated. 

Reference to Table VII shows that the center 
of gravity of the diphenyl levels depends on 
benzene energy level terms and also on diphenyl 
terms arising from the interaction between 
benzene rings. The numerical value (Table IX) 
of the benzene terms }3(£:,+£3,)—Eo and 
3(Eo,+ 3s) — Eo are obtained from the paper by 
Goeppert-Mayer and Sklar. 
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TABLE VII. Values of Hii, Hii —Hj;, }(Hii+Hj;) — Eon. 


Numbers in table are 


xA 


=2Eo + Ao +Go —Do 


+A +G@ —D 
=Eo +E3s +A —D 
=E\s —E3s 
=}(Eis +E3s) —Eo +A —Ao 
+4 —Go —(D —Do) 


+A +6 —D 
+E3s +A —D 
Hii-H3-3 =his 
+H_3-3) —Eop=4 (Eis +E3s) —Eo +A —Ao 
+G —Go —(D —Do) 


+A +G —D 
+E35 +A —D 
—E3g 
=}(E2s +E3s) —Eo +A —Ao 
—(D —Dpo) 


=Eo +A —D 
=Eo+E3s +A +@ —D 
Ho-2—-H4-4 —E3s 
+H4-4) —Eop=4(E2s +E 3s) +A —Ao 
+G@ —Go —(D —Do) 


[0404, 23 
(lala, 1p 


to 


—120 —2 —2 


-1-7 -7 1 
—5 13 13 5 
“3-3 -3 3 
3-3 -—3 3 


-1 -4-l1 3 6 6 
—1 —22 —35 —3 -—6 6 

1 —22 —35 —3 -—6 6 


In the course of this investigation the benzene 
energy levels were recomputed including all 
energy terms arising from non-neighbors. It was 
found that the levels computed in GMS were too 
low by approximately 3 ev as given in Table X. 

Goeppert-Mayer, in a private communication 
has indicated that similar results were obtained 
in a recalculation of the benzene levels including 
the effect of neighbors and neighbors once- 
removed. The new calculations are in poorer 


agreement with experiment in benzene and if 
used to determine the center of gravity of the 
diphenyl levels will cause poor agreement with 
experiment. Consequently we use the original 
values of 8.0, 5.8 and 5.0 for 'Tys, 'I'9, 1I'i0 in the 
diphenyl! calculation as the best estimate of the 
experimental values for the benzene energy 
levels. 

The term A—Ap+G—Go+D—Dy, is common 
to all energy levels. 


418 
= 
- 
Ay 2 2 -1 -1 1 -1 1 
Hyg 2 2 -i -1 -1 1 -1 -1 1 1 
Hy —Hs3 3-2 -1 -1 1 
3(Hi +H33) —Eop 
-2 1 1-1 1 1 1 -1 
8s 8 4 = 
-1 -1 -1 -1 4 & 4 1-1 1 
Hee 2 -1 1 -1 1 -1 1 1-1 
Hu 2 -2 2 1 
Hop $ -8 -8 1 #1 =1 
3(H22 +H44) —Eop 
2 3 1-8 1 
—1 1-2 12 1-1 2-1 1 
q TABLE VIII. Values of H;;. Numbers in table are coefficients t) 
q 7 a | 
q 
q -1 1-2 42 8 1-1 1-8-1-1-1 @ 3-8-6-832 3 1 
Base 1-3 4-8-8 1-1-1 1 2 8-8-1-1 2 3 -1 
Hes= 11 2 4-8 2-1 -1-1-1 1 1-1-1 1 3 -1 
4 
st 
q 
4 


—-24 


SPECTRUM OF DIPHENYL 


coefficients of terms appearing at top—Multiply by A or A?, as indicated. 


18 —18 
18 -—18 —18 18 


= 


-~6 6-6 6 3-3 7 16-16 16 16 —16 —16 —16 —16 
—6 6 -6 6 3-3 -15 20-16 2 2 2 20 —16 —16 — 
—62 —6 —36 —6 31 3 14 —11 —18 —18 —18 —18 —18 —18 —18 —18 
—62 —6 —36 -—6 31 3 14 11 —18 —18 —18 —18 18 18 —18 —18 


XII. DISCUSSION OF TABLE IX 


The Coulomb integrals yi4!g’ are approximately independent of the values of / and /’. To under- 
stand this, proceed as follows: 


+2 cos IV+IV V+V VIFVIT) 


+2000 IV+I1 V+IV VIFV I+VI 
42(—1)4{1 IV+ITV-+IIL VI, (103) 
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are 
=a 
3 12 5.5 -18 —48 
1 1 44 3.5 —16 1 3-6-8 3 7 7 —16 —16 5 5 10 10 40 60 
1 1 4 4 3.5 —16 1 3-8 -6 3 7 7 -16-16 5 5 10 10 40 60 
4 —4 2-2 7 7-7 -7 16 = 16 —40 —60 40 60 
1 1 -3 -8 2 23.5 —5.5 -16 18 1 488 3 -7-7 3 35 35 35 35 -—-8 -—8 -8 -8 5 5 10 10 20 30 20 30 
-1 5 8 8 7.5 —20 1 3 -6 32 3-15 —15 20 «20 5 5 10 10 40 60 
-1 5 8 8 7.5 —20 1 3 32-6 3 -15 -—15 20 20 5 5 10 10 —40 —60 
8 —8 —38 38 —15 —15 1 15 2. 20 —20 —20 40 60 40 60 
-1 5 -3 -—4 4 475 —5.5 —20 18 1 48 3 13 18 3 -75 -75 -75 -75 10 10 10 10 5 5 10 10 —20 —30 20 30 
-1 3 6 6 5.5 —18 1 -3 -6 0-3 11 —Iil -18 18 —5 —5 —10 -—10 
-1 3 6 6 5.5 —18 1 0-6 -3 11 
6 -6 6 11 -11 -11 11 —18 18 18 —18 
-1 3 -3 3 35.5 -5.5 -18 18 1 48 -3 -3-3-3 55 -55 55-55 —9 9 -9 9 —5 —10 —10 
1 3 6 6 5.5 —18 1 -3 -6 0-3 -11 18 —18 —5 —10 —10 
1 3 6 6 5.5 —18 1 —3 0 -3 -ll il —5 —10 —10 
6 - —6 6 il ll —1l 
1 3 -3 -—6 3 355 -18 18 1 48 -3 -—3 -3 —5.5 5.5 —-5.5 5.5 9 -9 9 -—5 —5 —10 —10 
ents of terms appearing at top—Multiply by A or A? as indicated. 
es =e ssTS SSIS Sx 
30 620 —30 —20 
30 —20 30 20 
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TABLE IX. Numerical values of various integrals and energy level terms. 


Ess) — Eo 6.9 volts — 1.30 volts 
4( Fo, +E3s) — Eo 6.5. QA? —1.32 
7040B 3.49 0.180 
yOAlB 3.51 0.181 
‘ yOA2B 3.54 0.179 
y1aAlB 3.49. E35°A? 0.186 
3.52 GA? 0.192 
—0.09 0.191 
A—Aop+G—Go +0.09 A’A? — 0.234 
6040B 0.022 0404, 0.124 
6041B 0.023 0404, 0.148 
6042B 0.018 lala, 0.179 
614lB 0.021 —0.059 
—0.045 Re — 0.044 


(1424, —28—128] 
Q! — 


w0 A0B 


@W0A1B — 3.57 Yo1 9.02 
w042B —3.38 Yo2 8.78 
wlalp 3.44 vu 8.95 
wi4—1B 3.86 v12 8.43 
W1A2B 3.18 3.12 
w-142B 3.71 502 1.67 
w2A2B — 3.00 
w24—2B — 3.68 
3.67 
042A —1.13 
—1.18 
AlA—2A 0.65 
AlA2A 3.75 
—0.99 


with a similar expression holding for | ¢iz’|*, we get 


4nl, 
{| cos (=) Si24(v) +2 cos (=) Si34(v) +(— 


Lp’ 
) +2 cos ( 


x| Seta) +2 cos ( 


where 


VI XII 
S2= Ke; Sz’ = K? 
K;=I1 K,;=VII 


VI XII 
S124 = K Kins; K Kin ( 
K;=1 K;=VII ( 
(105) 
VI XII 
S134 = KK S133 = Zz. K 
K;=I é 
VI- XII 
K;j=I K;=VII 
Kos = Ksy2= Kays= I, Key 


and 


Kiy2=Kuys=VII, Kiys=IX. 


420 
1424, — 0.090 —0.112 
—0.359 A? —0.109 
—1.20 —0.0947 
—3.47 9.11 
3 | 
152 
11 
N_i-3 — 0.00295 
Nu = N_o4 — 0.00937 
in 
i 
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The integral (104) can be separated into a number of terms, a typical one of which is 


(106) 


In accordance with our previous discussion an integral such as 


f f 11 (u) VIL (u)d7,d7, 
is replaced by ; 
f f ete. 


However, if (VII, VIII) were replaced by VII? instead of K;s?, (VII, 1X) by VIP’, etc. (106) would 


become 
| cos ( ) +2 cos +(—1)!*’s Af. fro (107) 


where s; is the overlap integral between neighbors in the benzene ring, S2 between neighbors once- 
removed and s3; between neighbors twice-removed. The expression in brackets in (107) is exactly 
oiz’ (Eq. (19) GMS). In addition to the term (106) there are similar terms in which I?(v) is replaced by 
and (—1)!4I(v)1V(v) successively. These terms combined with (107) result in the term 


Evidently, under these assumptions (104) may be written as (see 105) 


e 
360 f (110) 
so that 


1 
f [seo serwarar, = 3.52 volts. (111) 


Comparison of (111) with the values given in Table IX, which were obtained by using the fictitious 
centrally-located 2px function approximation, shows the extent of agreement between the two types 
of approximations. 

The diphenyl exchange integrals 6747’ are all very small, approximately 0.02 volt. This results 
from the fact that the double exchange between rings gives rise to integrals whose numerical values 
are proportional to the product of two overlap integrals (Table V). As a consequence, D—D, is 
completely negligible. 

The \, integrals (Eq. (88)) are practically independent of the value of J. As in the benzene case, 
we set the potential term 


=X Hxa(v)— | (112) 
Kp Kp 


i.e., the potential acting on an electron in ring A as a result of the carbon ions in ring B is the sum- 


422 ALBERT LONDON 


mation of the potential Hxg due to each individual carbon ion and the same as that due to the 
neutral carbons atoms in ring B; }>kg Hxg(v), minus the potential caused by the 27 electrons in 
ring B. From (88), (103), and (105) we have 


1 Qala, 
cos Si24(v)+2 cos (~~) Si34(v)+(— Hx3(v)dz, 
60, 6 6 Kp 


For reasons similar to those used in obtaining the approximate value of yisig’, (113) becomes ap- 
proximately 

1 1 e 

f seo > ff [seo — 21.23 volts. (114) 

6 Kp 6 Top 


Note that in Table IX, \2—); was calculated directly by an appropriate grouping of differences of 
small terms. 
- The location of the energy levels also depends on the square bracket terms such as [— 1424, 28—12 ] 
which are relatively small in magnitude. The bracket terms such as [—14—12, 224] are really 
exchange integrals since each electron is partially in ring A and in ring B. Their approximate value 
is 0.02 volts; hence, in combination with the factor of } they may be neglected. 

The terms which contain A as a factor are all of the form wisig’ (Eq. (90)). These are computed 
in the following manner: From (92) and (112) we have 


Ka vp 


also from (91) and (3) 


1 
€lalp’= (or)! four] +e"VI11)+ (vy) ]dz,. (116) 
: Kp 
Equation (116) may also be written as 


1 
four { [Tp(v) + HAyn(v) 2(r) + Ayin(v) Je@’ VIII + 


Kp 


1 
ol 
Since VII(v), VIII(v), --- are atomic wave functions we may write (GMS) 
[T(v) + Avn(r) VII = W2,VI1(r), (118) 


etc., where W2, is the energy of a 2p electron, in the unexcited carbon atom. Thus, the first integral 
of (117) becomes equal to Ai4ig’W2,. Expanding }>xsHkz in the second integral of (117) and adding 


| 
| 
1 


the 
s in 


113) 
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114) 
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the resulting value of €741g’ to Xi41g’, we obtain the final result 


Wart f > Hx(v) 


K=!I 


e[K=XII 


where >~’ is a summation in which all terms are zero other than those for which Kz in ¢ig’ has the 
same value as Kz in Hxzg, as in the third integral of (117). 

For the value of W2,, the experimental value”® of — 11.24 volts obtained from ionization potential 
data was used. Since W2, occurs with the factor Ai,ig’, the value of wiaiz’ is not strongly dependent 


on Wp. 

The evaluation of terms containing A? as a 
factor involves the determination of the value 
of integrals not previously considered such as 
Nala’ (Eq. (93)). An expansion similar in nature 
to (115) for Az47g’ holds for these integrals. 

Other A? terms which occur such as Q, Q, 
Fo’, Exs’, Eos’, have rather large numerical 
values since they depend on the energy of the 
orbitals, i.e., w; (Eq. (89)), or on benzene in- 
tegrals (Eq. (75)). However, in combination 
with the factor A’, their contribution is moderate. 

The effect of the bracket integrals, such as 
[0404, 1g—1s] when multiplied by the factor 
A’ is completely negligible. 

Terms in H;; such as woA*, w:A?, and wd? 
strongly influence the value of H;; because of 
their relatively large values and the large coef- 
ficients with which they appear. However, their 
contribution to the value of H;; is generally 
partially compensated for by the benzene y and 
6 integrals and by diphenyl y integrals. 

Reference to Eq. (24) shows that the energy 
levels also depend on the value of N;;. From 
Table IX, evidently NV;?<1, so that in expressing 
W:;*—Eop as given in Table VII, the coefficient 
1/(1—N;;;*) is simply taken as 1. 


TABLE X. Benzene energy levels—ev. 


8.0 
IT, 5.8 8.8 
5.0 7.3 
Tho 3.0 7.2 

2.2 5.8 

15 4.5 


* A. E. Ruark and H. C. Urey, Atoms, Molecules, and 
Quanta (McGraw-Hill Book Company, Inc., New York, 
1930), p. 280. 


Using the numerical values listed in Table IX, 
and Eq. (24), the energy levels depicted in Fig. 4 
are obtained. From considerations following Eq. 
(24) it is clear that if both H;; and N;; were 
zero, the energy levels would simply be W;; and 
W,;. Thus the 5-volt level in benzene would 
result in levels and W_2». However, 
because of the interaction between ®,? and #4”, 
the value of We: shifts slightly to Was, while 
W_s_2 becomes W_2_4-. Thus the 5-volt level in 
benzene has been split into the two levels, 4.9 
volts and 5.2 volts. It will be recalled that the 
5-volt level in benzene gave a forbidden transi- 
tion to the ground state. In diphenyl, however, 
the 5.2-voit level gives an allowed electronic 
transition, so that the weak 5.0-volt transition 
in benzene gives rise to a strong transition in 
diphenyl. 

Similarly, the forbidden 5.8-volt transition in 
benzene gives rise to an allowed 6.6-volt and 
forbidden 5.4-volt transition. The 8.0-volt 'T12 
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level which gives an allowed transition in 
benzene gives rise to four levels two of which 
have allowed transitions. 

Experimentally, the ultraviolet absorption 
spectrum of diphenyl in the vapor phase”! shows 
a minimum absorption at 5.9 volts and also at 
6.9 volts with a maximum in the absorption 
band at about 6.4 volts. Also a broad maximum 
occurs from about 7.3 volts to 7.6 volts which 
was the shortest wave-length used in these 
experiments. The ultraviolet absorption spec- 
trum of diphenyl in the vapor state has been 
investigated by Beck” at the Catholic University 
of America. These measurements show the 
existence of a very broad absorption band from 
about 2200A to 2800A, i.e., from about 5.6 to 
4.4 volts. Thus, the levels shown in Fig. 4 are in 
good agreement with the experimental observa- 
tions. 

In Carr and Stiicklen’s paper the absorption 
bands are very broad and show very little vibra- 
tion structure. The band observed by Beck is 
completely diffuse and shows no vibration 
structure. It is believed that the diffuseness in 
the observed bands may be partly due to the 
overlapping of the several electronic levels 
involved. 


XIII. SUMMARY AND CONCLUSIONS 


The diphenyl wave functions are set up as 
linear combinations of molecular benzene func- 
tions, and in such a way as to have the same 
symmetry as diphenyl. These wave functions are 
next antisymmetrized with respect to the inter- 
change of any two electrons by including the 
appropriate spin factor in the functions. As a 
consequence of this procedure the secular de- 
terminant automatically factors into four quad- 
ratic equations from which the diphenyl energy 
levels may be determined. The center of gravity 
of the dipheny] (singlet) levels is shown to depend 
on the benzene singlet levels previously com- 
puted by Goeppert-Mayer and Sklar. In the 
calculation of the position of the levels, a large 
number of integrals occur as a consequence of 
the interaction between the two benzene rings 
of diphenyl. All integrals are considered in the 


21 E. Carr and J. Stiicklen, J. Chem. Phys. 4, 760 (1936). 
22 Private communication. 
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evaluation of the energy levels up to terms con- 
taining A’, the square of the overlap integral, as 
a factor. The integrals which must be evaluated 
are of two types, electronic, involving the inter- 
action between two 27 electrons, and ionic or 
nuclear involving the interaction between a 2p7 
electron and a carbon ion or 2pz-less carbon 
atom. The numerical value of integrals of the 
required type for the case where only interactions 
between neighbors are considered, with the ex- 
ception of Coulomb integrals, are available in the 
literature. In this paper the restriction that the 
calculation be confined to neighbors, was 
dropped. The effect of all non-neighbors was 
determined by suitable approximation formulae. 
The resulting energy levels are in fairly good 
agreement with the experimental observations. 
All of the numerical values of the energy 
integrals are based on internuclear distances 
computed on the basis of a flat plane structure for 
diphenyl. The good agreement between theory 
and experiment therefore supports the hypothesis 
of the coplanarity of the two benzene rings. 
The method carried out in this paper for 
progressing from a known solution for a rela- 
tively simple molecule, such as benzene, to the 
problem of obtaining a solution for a more com- 
plicated molecule, such as diphenyl, is believed 
illustrative of the procedure which may be 
carried out in similar cases or in building up to 
more complicated molecules. Furthermore, it has 
also been demonstrated here that the method of 
antisymmetric orbitals is suitable for use in 
determining the splitting of known energy levels 
which may be caused by: the presence of a per- 
turbing influence or interaction of some kind. 
In carrying out the calculations of the energy 
levels, it was found necessary to consider the 
electronic and nuclear interactions separately. 
A method which would appear to have inherently 
greater accuracy would be to consider the elec- 
tronic and nuclear interactions jointly. In this 
case, the energy integrals which occur are essen- 
tially due to the interaction between a 2p7 
electron and a neutral atom. Since these integrals 
have small numerical values, the levels would be 
obtained as a summation of small terms. In the 
present method the levels are partially dependent 
on differences between electronic and nuclear 
integrals which are larger in magnitude than the 
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electronic-neutral atom interaction. However, it 
was found difficult to carry out the proposed 
method, since a one-to-one matching between 
electronic and nuclear terms was not possible in 
all cases, and since the values of the resulting 
penetration integrals had not previously been 


evaluated, with the exception of neighbor inter- 
actions. 

The author is indebted to Professors K. F. 
Herzfeld and A. L. Sklar for many helpful dis- 
cussions and much assistance and guidance in 
connection with this work. 


APPENDIX I 


In order to bring out the meaning of the permutation operator Pag given in Eq. (37) we carry 
out the details of the proof sketched in the above for Eq. (50). 

Consider the one-electron orbitals which make up the integrand of 6:/11° [see 0:g1:(0) in Eq. (1A) 
in Appendix II]. We do not write the orbitals for electrons (1), (2), (7), and (8) since they will remain 
unchanged in future operations. Thus, writing the terms of the integrand only, we have 


in 


a B a .B a B 
X (1’) 
Lap — Ba] 


a B a B a B B-B 


a [eB—Ba] ao B a B 
X 


which is obtained by taking the complex conjugate of everything inside the square brackets. 
Now let us change the name of the integration variables throughout the entire integrand in the 
following manner: 305, 46, and (2’) becomes 


a B a B a B [aB— Ba] 


B a B a B 
The integrand for b4h1:° is 


a B a B a B [aB— Ba] 


[aB—fa] a 


a B B a B 
X  (4’) 


Terms in 41,’ will be different from zero only if they contain no more than two unmatched orbitals. 
In accordance with the definition of Pag, in Eq. (4’) we perform the matching permutations 35, 
9911, 10<+12, thus leaving the unmatched pair 


J. (5’) 
Evidently the above permutation transforms Eq. (4’) into (3’) thus showing that 
= 


The unmatched pair (5’) represents an identity permutation in the class P,g since it involves one 
electron in ring A and one in ring B. Furthermore it gives rise to an integral which results from the 
interaction between an electron in ring A and one in ring B. Consequently, it will be seen that the 
class of Pag permutations leads to integrals which represent the interaction between rings. 
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APPENDIX II 


Value of the integrals a,g;;(0) in terms of the one-electron orbitals 


1 a 


X 
a B a B a B 
X 
a B a B LaB—Ba | 
X ]} dr, 


1 a a a 


P(6) 


Ba] 


a B a B 
X 
a B a B a B 
X P(0) [G04 (1) 
a B a B Lap — Ba] 
X dz, 


i 
8 
a B a B LaB—Ba] 
X 
a B a B LaB— Ba] 
X P(0) [G04 (1) p04 (2) 614 


a B a B a B 
X ]} dr, 


1 a a 


P(6) 


[ap — Ba] 


a B B 
a B a B Lap — Ba] 
X P (6) [04 (1) G04 (2)-14(3)-14 (4) 
a B a B a B 
X 1)@-1(12) ]}dr, 


1 a a 


a B a B [ap — Ba] 
X 
a B a B a B 
X P(8)[ 


a B a B Ba } 
X ]}dr, (1A) 
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8 Pe) 


a B a B [aB—Ba] 
X 


a B a B a B 
Lap — Ba | 


8 


.«.4 [ap —Ba] 
[aB— Ba] 


B a B 
1)@-12(12) }}dr, 


[ap — Ba] 


B 
0)@-12(11)¢-22(12) 


Lap — Ba | 


14(3)- 14(4)14(5) (6) 


B B a B 
1)@-15(12) ]}dr. 


In the above set of integrals, for 06=Hp*, the axqii(@) become a;h;;° and P(Hp*) becomes Pz; 
for Hp", the a:gii(@) become a,f7;;" and P(Hp") becomes P; while for 6=/, the a:q;:(@) become 
a,N;; and P(1) becomes P. The spin function which corresponds to each one-electron orbital is 
written on top of the orbital, while the bracket [a8—a] over say electrons 11 and 12 is read as 
[a(11)8(12) —6(11)a(12)]. The coefficient 1/V2 which appears in front of each of the [a8—Ba ] 
bracket is already accounted for in the coefficient } in front of the integral. 

In a similar manner we may write the integrals a,q; ;(4) 


8 


a B a B [aB— Ba | 
X 


[ap — Ba'] 


x |}dr. 


8 Pe) 


a B a B [aB— Ba] 
X 


a B a B a B 
14(S)- 
LaB— Ba] 
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1 a a a 


[ap — 


B a B 
« 
X P(6) [G04 (1) G04 | 


1 a a a 


P() 
a B a B [aB— Ba] 
X 
a B a B [aB— Ba] 
(0) 


In addition to Eqs. (1A) and (2A) we need an expression for determining the ground state energy 
and normalization integral for diphenyl. The following integral serves this purpose. 


a B a B B 
X (3A) 


For 6=Hp*, P(@)=P, proceeding in the usual fashion we can split off permutations confined to 
the A ring and also those confined to the B ring to give 


Io(Hp*) = | (4A) 


where /o(Hp*) is the same integral as (3A).with P(@)=P.az and @=Hp’. 
The electronic energy for the ground state of diphenyl, Eop’, is given by 


(5A) 


Similarly the nuclear energy is obtained from Jo(@)/No by replacing @ by Hp" and P(@) by P, while 
No=Jo(1) where P(1) =P. 
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The effect of the presence of 10-* M thorium chloride on the zeta-potential at the interface 


between vitreous silica and three solutions of potassium chloride has been determined by the 
streaming potential method. This trace of thorium changed the zeta-potential, which was 
strongly negative with the same solutions without added thorium, to positive values of smaller 
magnitude. The theory of Langmuir by which can be calculated the thickness of the wetting 
film owing to the zeta-potential has been extended to apply to these 4-1, 1-1 electrolyte mix- 
tures, and the effect of this wetting film on the capillary rise of these same solutions has been 
studied. With only one solution (10-° M KCl) containing added thorium was the zeta-potential 


low enough so that the wetting film was negligible in the capillarimeter of Jones and Ray. 


INTRODUCTION 


ONES and Ray’ have accurately measured the 
relative surface tension of water and of solu- 
tions of electrolytes with their differential capil- 
larimeter, and they observed that the first trace 
of added salt apparently decreased the surface 
tension below that of pure water. On the surface 
tension-concentration curve for potassium chlo- 
ride solutions, the minimum surface tension was 
found at about 0.001. concentration, where it 
is 99.982 percent of the surface tension of pure 
water. Langmuir? has suggested that the mini- 
mum is not real, but is caused by variation in the 
thickness of the wetting film on the capillary wall, 
dependent on the concentration of the solution. 
He advanced a quantitative theory, developed 
more rigorously by Jones and Frizzell,’ which 
gives values for the wetting film thickness large 
enough to decrease appreciably the effective 
diameter of the capillary. With a particular 
solution, the thickness of the film depends on 
the magnitude of the zeta-potential, which is the 
electric potential at the interface between the 
silica capillary and the liquid. With water and 
with dilute solutions the film may be several 
hundred angstroms thicker, if the zeta-potential 


* Contribution from the Mallinckrodt Chemical Labo- 
ratory of Harvard University. 

‘Grinnell Jones and Wendell A. Ray, J. Am. Chem. Soc. 
(ise. (1937); 63, 288 (1941); 63, 3262 (1941); 64, 2744 

*Irving Langmuir, Science 88, 430 (1938); J. Chem. 
Phys, 6, 873 (1938). 

* Grinnell Jones and Laurence D. Frizzell, J. Chem. Phys. 
8, 986 (1940). 
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is large enough, than it is with more concen- 
trated solutions. A film of this thickness is 
sufficient to account for the apparent surface 
tension minimum observed without »admitting 
the existence of a true minimum in the surface 
tension. In order to test the Langmuir theory, 
the zeta-potential of silica in contact with solu- 
tions of potassium chloride was measured and 
reported by Jones and Wood.‘ It was found that 
correction of the apparent relative surface tension 
for the thickness of the wetting film as computed 
by the Langmuir equation eliminated the mini- 
mum in the surface tension-concentration curve. 

An earlier test of the Langmuir theory was 
made by Jones and Frizzell,? but their results 
were inconclusive. Langmuir? had suggested that 
if the zeta-potential were to be lowered to zero, 
or nearly so, by the addition of thorium ions, the 
film thickness should be decreased, and the 
minimum eliminated. The best evidence, accord- 
ing to Abramson,® indicated that the zeta- 
potential of glass was near zero in a solution of 
10-* M thorium ion. No similar information 
existed regarding silica. Jones and Frizzell made 
numerous measurements of the relative surface 
tensions of thorium chloride solutions in the 
concentration range 10~* to 5X 10-5 M, but were 
unable to deduce from these measurements the 
particular concentration at which the zeta-poten- 
tial of silica was zero. They stated that direct 


4 Grinnell Jones and Lloyd A. Wood, J. Chem. Phys. 13, 
106 (1945). 

5H. A. Abramson, Electrokinetic Phenomena (Chemical 
Catalog Company, New York, 1934), p. 186. 
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TABLE I. rE cme (millivolts) with 
1.0X10-§ M ThCl,—1.00X 10-5 M KCl. 


Average 


Number 
Time of meas- «(X10*) 
deviation 


(days) urements (av) 


E/P 
(av) 


Capillary I 
2.035 
2.252 
2.208 
2.211 


3 
< 


—6.8 
-7.1 
—44 
— 20.9 
—19.1 
+14.1 
+7.0 
+10. 


00 00 CO COCO COCO > 


2 
3 
4 
8 
9 
10 
16* 
17 
18 


+37.1 
+39.9 
+34.4 
+33.0 
+35.7 
+37.5 
+30.5 
+18.2 
+23.6 


~100 00 


2 
3 
4 
8 
9 
10 
16* 
17 
18 


* Preceding this set of measurements, a new sample of solution was 
et in the apparatus, as it was supposed that thorium may have been 
ost from the solution in apparatus I, due perhaps to adsorption. 


measurement of the zeta-potential was needed to 
settle this question. However, they redetermined 
the apparent relative surface tension of potassium 
chloride solutions, using 10~§ M thorium chloride 
solution as the solvent component rather than 
pure water, on the presumption that if the 
minimum were eliminated, it would be evident 
that the addition of thorium had decreased the 
zeta-potential enough so that the wetting film 
was negligible, and that Langmuir’s theory was 
correct. They found that the minimum was not 
eliminated, although it was not as pronounced 
as without the thorium. The question then arose 
as to whether the zeta-potential was brought 
near enough to zero by the added thorium, and 
the work herein reported was undertaken to 
provide an answer. The zeta-potential of silica 
in contact with solutions of 10-*, 10-‘, and 
10-* M potassium chloride, each 10-§ M in 
thorium chloride, has been determined by the 
measurement of streaming potentials. 


EXPERIMENTAL 


The streaming potential apparatus has been 
fully described in an earlier publication,‘ and in 
obtaining the measurements reported below, no 


WOOD 


changes were made either in the apparatus or in 
the experimental procedure. _ 

Thorium chloride was prepared following the 
method described by Jones and Frizzell.* A solu- 
tion of 0.020 M concentration was prepared with 
conductivity water and stored in a non-sol glass 
bottle. No turbidity appeared on long standing. 
The 10-® M solutions were prepared by dilution 
of this stock solution. The potassium chloride 
was recrystallized from analytical reagent KCl, 
and the solid was weighed ‘out for addition to 
the proper quantity of the 10-* M thorium 
chloride solution. Conductivity water from a still 
with block tin condensers was used in preparing 
the solutions, and when freed from carbon di- 
oxide, had a specific conductance of less than 
0.8X10-* mho per cm. The solutions used were 
made up in contact with laboratory air, but after 
the solutions were put in the streaming potential 
apparatus, carbon dioxide-free air was bubbled 
through them until the dissolved carbon dioxide 
was removed. 

In referring to these solutions in the discussions 
following, it will be considered sufficient for the 
sake of brevity to refer only to the KCI concen- 
trations, with the understanding that each con- 
tains 1.0X10-* VW thorium chloride in addition. 
The temperature in all experiments was 25°C 
+0.03°. 

The first solution measured was 1.00X 107° 
KCl and the results obtained are given in Table I. 
The “Time” is the number of days that the 
solution had been in contact with the capillary 
when the given measurement was made. The 
“‘Number of measurements” gives the number of 
individual streaming potential measurements 
which were averaged to give the reported values 
of « and E/P. x is the effective specific con- 
ductance of the liquid in the capillary, and 
includes any contribution of surface conductance. 
E/P is the ratio of the streaming potential £ to 
the pressure P which produced it. ¢ is obtained 
from the equation of Helmholtz® 


in which 7 is the viscosity of the solution and D is 
the average dielectric constant of the liquid in 
the ionic double layer. Both these quantities are 


6H. Helmholtz, Wied. Ann. 7, 337 (1879). 


430 
ol 
al 
d 
th 
pi 
tl 
la 
Capillary II ( 
: 1 al 
4 
di 
a 
of 
T 
pr 
di 
pr 
n 
gi 
sc 
Ir 
al 
tl 
W 
ac 
of 
of 
pt 
of 
m 
K 
th 
Ww 
1g 


or in 


ig the 
solu- 
with 
glass 
iding. 
lution 
loride 
KCl, 
on to 
orium 
a still 
aring 
mn di- 
than 
were 
after 
ential 
bbled 


oxide 


ssions 
yr the 
ncen- 
con- 
ition. 

25°C 


<10-° 
ble I. 
t the 
illary 

The 
der of 
nents 
ralues 

con- 
and 
ance. 

E to 
ained 


1 Dis 
id in 
are 


assumed to be the same as those for pure water 
(n= 0.00895 and D = 78.49, according to Dorsey).? 
It should be mentioned that the values given for 
¢ were obtained by taking the average of the ¢’s 
obtained from each individual measurement of « 
and E/P, so that in general the ¢ given is slightly 
different from what would be obtained by using 
the average values of « and E/P for com- 
puting it. 

When the first measurement was made, after 
the solution had been in contact with the capil- 
lary for two days, the zeta-potential of capillary | 
(diameter: 0.032 cm) was found to be —11 mv 
and that of capillary II] (diameter : 0.017 cm) was 
+37 mv. The zeta-potential varied somewhat 
during the succeeding eight days, but it was 
apparent that the potentials of the two capil- 
laries were not reaching agreement. These results 
seemed to indicate that the concentration of the 
thorium might have changed in one or both sets 
of apparatus, perhaps owing to unequal amounts 
of adsorption of the thorium on the walls. 
Therefore, both were refilled with a second 
portion of the solution. The capillaries were not 
disturbed. It was then observed that the zeta- 
potential of both capillaries was positive although 
not yet in very good agreement. These values are 
given in Table I beginning 16 days after this 
solution first came in contact with the capillary. 
Inasmuch as there was nothing to indicate that 
any one of the zeta-potentials obtained on the 
16th, 17th, and 18th days was more probable 
than another, the average of these six values 
was accepted as the most probable value. This 
accepted value is +17 mv and the probable error 
of a single measurement is +6 mv. The addition 
of 1.0X10-* M thorium chloride to the 10-5 M 
potassium chloride has thus reversed the polarity 
of the zeta-potential and greatly decreased the 
magnitude, changing it from —148 mv, reported 
earlier,4 to about +17 mv. 

The next solution measured was 1.000 X 10-4 M@ 
KCI, also with 1.0X10-§ M ThCly. The capil- 
laries were treated with boiling aqua regia, and 
thoroughly rinsed and steamed with conductivity 
water. The other parts of the apparatus were 
rinsed thoroughly with the solution, but not with 


™N. E. Dorsey, ———— of Ordinary Water-Substance 
orporation, New York, 1940), pp. 
an 
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acid or water, to avoid removal of the thorium 
ions which may have been adsorbed on the walls 
while the previous solution was in them. The 
presumption was that while the relatively small 
surface of the capillary walls would not adsorb 
enough thorium to change the concentration 
appreciably, the much larger surfaces of the 
other parts of the apparatus might do so, as 
they evidently did with the 10-> M KCI solution. 
Whether this was true or not, the difficulties 
encountered with the first solution did not arise 
again. A summary of the zeta-potential measure- 
ments with the 10~* M KCI solution is given in 
Table I], and it will be seen that the potentials 
observed were in better agreement than with the 
10-5 M solution. After five days, the potentials 
of the two capillaries were 4 mv apart, and the 
accepted value was +57 mv, the average of 
these. The addition of the thorium has thus 
changed the zeta-potential from —135 to about 
+57 mv for silica in contact with 10~* WM KCl. 

Following this series, the apparatus was thor- 
oughly cleaned to remove the thorium. After 
repeated treatment of the capillaries with aqua 
regia, concentrated nitric acid, and hot saturated 
KCI solution, zeta-potential values of —145 mv 
for capillary I and — 140 mv for capillary II were 
obtained with 1.000 X 10-4 M KCI without added 
thorium. This. demonstrated that the thorium 
ions could be and were removed from the 
capillary surface, and that the surface was 
substantially the same as it was before it came 
in contact with the thorium, when —135 mv 
was obtained with the same solution. 


TABLE II. Zeta-potentials (millivolts) with 
1.0 10-* M ThCl,—1.000 x 10-4 M KCI. 


Number 
Time of meas- K E/P i? Average 
(days) urements (X10) (av) (mv) deviation 
Capillary I 
1 4 16.24 34.06 +53.4 0.1 mv 
2 8 16.04 33.23 +51.5 0.2 
3 8 16.06 31.73 +49.2 0.2 
4 8 16.14 32.54 +50.8 0.7 
5 8 16.15 35.29 +55.1 0.2 
Capillary II 
1 4 1642 45.28 +718 04 
2 8 16.23 4064 +63.7 0.5 
| 8 16.27 39.54 +62.1 0.2 
4 8 1635 37.81 +59.7 
5 8 1636 3749 +592 0.5 
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Then the third solution with added thorium, 
1.00 10-* M KCI and 1.0X10-§ M ThCly, was 
placed in the apparatus. The various parts were 
of course thoroughly rinsed with the solution 
before assembly. Table III gives the average 
zeta-potentials, and the accepted value was +67 
mv, the average of the last pair of measurements, 
which were closest together. The addition of the 
thorium has therefore changed the zeta-poten- 
tial of the 10-* M KCI solution from —112 to 
about +67 mv. 

With none of these solutions was the study of 
the change of the zeta-potential with time carried 
on as long as with the KCI solutions reported 
earlier because the purpose here was to obtain an 
idea of the value of the zeta-potential as it most 
probably was in the surface tension capillarimeter. 
Since it is well established that a considerable 
period of time is usually necessary for an equi- 
librium potential to be obtained, and inasmuch 
as the solutions were not continuously in contact 
with the surface tension capillarimeter as long as 
these were with the streaming potential appa- 
ratus, the zeta-potentials reported here are prob- 
ably only an approximation to what they were 
in the experiments of Jones, Ray, and Frizzell. 

The chief experimental difficulty seemed to be 
in maintaining the thorium ion concentration at 
10-* M, observed in the experiments with the 
10-5 M KCI. Although there was no difficulty in 
preparing the solution, it is inevitable that some 
adsorption of thorium must occur at all the 
surfaces of the containing vessels, since thorium 
must be adsorbed rather strongly to be capable 
of so great an effect on the zeta-potential. Even 
though there may be actually only a small 
amount of thorium removed from the solution 
by adsorption, this small quantity may be an 
appreciable fraction of the thorium in the 10-* M 
solution. This probably explains the lower order 
of accuracy of the measurements with the solu- 
tions containing 10-§ M thorium chloride than 
was obtained in the measurements reported 
earlier. 

However, it is reasonable to assume that the 
zeta-potentials measured in the streaming poten- 
tial apparatus are essentially the same as they 
were in the surface tension capillarimeter with 
the same solutions, and the Langmuir theory 
will be applied using these data. 


TABLE III. Zeta-potentials (millivolts) with 
-1.0X10-§ MThCI,—1.000 x M KCI. 


Number 
Time of meas- K E/P c Average 
(days) urements (X<10®) (av) (mv) deviation 
Capillary I 
1 8 148.36 5.117 +73.3 0.3 mv 
2 8 148.31 5.088 +72.9 0.7 
3 8 148.30 5.064 +72.5 0.1 
4 8 148.25 4.865 +69.7 0.8 
Capillary II 
1 8 147.68 4610 +4658 0.8 
2 8 147.69 4.554 +650 0.8 
3 8 147.75 4.527 +646 0.4 
4 8 147.72 4.510 +4642 0.7 


It is evident that the 10-* M thorium chloride 
has not caused the zeta-potential of the silica 
with these solutions to become zero, but it has 
gone further and made the potential positive. A 
lower concentration of thorium would probably 
suffice to bring the potential to zero. Further ex- 
perimental work is needed to determine exactly 
what this concentration would be, but it is 
apparent that only a slightly lower thorium 
concentration would be needed for the 10-5 KCI 
solution. Since the zero potential was overshot by 
a much larger margin for the 10-* and 10-* M KCl 
solutions, it seems probable that a concentration 
of thorium much lower than 10-* would be 
enough to reduce the potential just to zero with 
these solutions. It seems reasonable to conclude 
that a different thorium concentration will be 
needed for each of these solutions to bring about 
a zero potential, because it is evident that less 
thorium would be required to reduce the potential 
of —112 mv of the 10-* M KCI to zero than the 
greater potential of —148 mv for the 10-5 M@ KCI. 
Therefore, a series of solutions with different 
thorium concentrations with each of the KCl 
solutions will have to be measured in the 
streaming potential apparatus, and when the 
proper amount of added thorium to produce a 
zeta-potential of zero is ascertained, the surface 
tension measurements will need to be repeated on 


the proper solutions. The extreme dilution of the. 


thorium with the attendant dangers of adsorption 
on the walls discussed above gives rise to the idea 
that a trivalent positive ion might be used more 
conveniently instead of thorium. 

However, according to the Langmuir theory, it 
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is not necessary for the zeta-potential to be 
reduced exactly to zero for the wetting film to be 
very thin. Langmuir considers the film to be 
attributed to the fact that the potential on the 
wall of the capillary causes a distribution of the 
ions in the solution such that the osmotic pres- 
sure of the water in the film is less than it is in the 
bulk of the solution. If the potential is high 
enough, the difference in osmotic pressure be- 
comes sufficient to create a thick film. As will be 
made clear below, there is a certain value of the 
zeta-potential for each solution which must be 
attained or exceeded before the difference in 
osmotic pressure becomes large enough to begin 
the formation of a thick wetting film on the 
capillary wall. If the zeta-potential is below this 
critical value, the wetting film will be attributed 
only to van der Waals or other short range forces. 
It will be shown below that the zeta-potential of 
the 10-§ M KCI solution was converted by the 
thorium toa low enough value so that the wetting 
film vanished, as far as the zeta-potential con- 
tribution was concerned, but that with the other 
two solutions, the potential was great enough so 
that the wetting film was appreciable. Therefore, 
it is evident that the thorium has not reduced the 
wetting film to zero, except for the 10-' M KCl 
solution. 

However, Jones and Wood‘ found that applica- 
tion of the Langmuir correction to the surface 
tension measurements, using the zeta-potential 
data, eliminated the minimum in the surface 
tension curve for the KCI solutions. Since the 
zeta-potentials for these same solutions with 
added thorium are now determined with suffi- 
cient accuracy, the Langmuir theory can again be 
tested by application of the correction to the 
surface tension data on these solutions. It is true 
that it would have been perhaps a neater test had 
the thorium concentration been such as to cause 
the wetting film to vanish, but since the thickness 
of the wetting film can be calculated according to 
Langmuir’s quantitative theory, a second critical 
test can be made with the zeta-potentials re- 
ported above. 


APPLICATION OF THE LANGMUIR THEORY 


Jones and Frizzell developed Langmuir’s 
theory only for a 1-1 electrolyte, so it is necessary 


to generalize their treatment somewhat to apply 
it to a mixture of electrolytes. Therefore, we 
begin with the general Poisson-Boltzmann equa- 
tion for the electric potential in an electrolyte 
solution in the region of a charged plane instead 
of the equation for a 1-1 electrolyte: 


dy 
—=—— ne exp (— 
ic Dp = 


where y is the positive electric potential, D is the 
dielectric constant of the solution, ; is: the 
number of ions of species 7 per cc, 2; is the charge 
on an ion of species 7, and ¢ is the magnitude of 
the electronic charge. d*f/dx* is used instead of 
Vy since the potential is assumed to vary only 
in the direction perpendicular to the planar 
interface. This equation is conveniently simplified 
by the following substitutions: 


n= ep/kT, 
0= kx, 


xk=[(43e?/DkT) n2;° |}, 


i=1 


(2) 
1/A=>. 
i=1 
dy 
Equation (1) then becomes 
dp 
p—=—A exp (3) 
dn i=1 


which on integration gives 
p?= (2A exp (—zm))+C. (4) 


To evaluate the constant of integration, C, we 
follow Langmuir and Jones and Frizzell, who say 
that at the air interface of the film when @=0, 
the potential is at a minimum, so that 7=m, and 
p=0. The value of C is then 


=-—2A n; exp (—zm). (5) 
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Combining Eqs. (4) and (5) and rearranging, we 
have 


1 


(6) 


[2A (exp (—2) —exp (—2m)) }! 


Since there is no known integral of this expres- 
sion, the general case can be carried no further. 
We now make the substitutions 


K=e, 
y= 


(7) 


and Eq. (6) becomes 


dé 
(8) 


dy 8 
It is now evident that if z; has integral values from 
—2 to +2, Eq. (8) can be expressed as an elliptic 
integral.* However, if any 3; is an integer greater 
in absolute value than 2, as with thorium ion, 
which has a charge of +4, Eq. (8) is than a 
hyperelliptic integral and cannot be integrated 
explicitly. It will be shown that a close approxi- 
mation for Ar, the thickness of the wetting film, 
can be obtained even so. 
Let 


Th*‘ be species 1: +4; 
K+ and H* be species 2: 22= +1; 
Cl- and OH~ be species 3: z3= —1. 


Substituting these values in Eq. (8), and also 
recalling that 1;=c;N where c; is the ionic 
concentration and N is Avogadro’s number, we 
then have for the solution of thorium chloride 
and potassium chloride 


dé —(10¢;+¢2)! 
dy (€2)§[ + (C3/c2)y 
+¢3K—)y?/ce 


It is evident that 


f de 
vo 


8 The reduction of such expressions to the Legendrian 
form is discussed in J. Edwards, The Integral Calculus 
(Macmillan and Company, Ltd., London, 1922), Vol. II, 
p. 578 and following. 


(9) 
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where Ar is the thickness of the wetting film, 
measured from the air interface in the direction 
of the wall; yo=K, the value of y at the air 
interface; and y; is the value of y when y=¢. 

Before proceeding further with Eq. (9), a 
digression will be made in order to find out how 
values of K are determined. Langmuir postulates 
that there is an equilibrium between three forces 
at the air-solution interface, provided that the 
film thickness is sufficiently large so that van der 
Waals or other short range forces do not play an 
important role in determining the thickness. The 
three forces in equilibrium at the level of the 
capillary meniscus are (1) the negative pressure, 
—(d—8)hg, tending to make the film thinner, 
owing to the force of gravity g, h being the 
capillary rise, d the density of the solution, and 8 
the density of air; (2) the pressure (d—£6)hg/2, 
owing to the surface tension of the cylindrical 
surface, tending to make the film thicker; and 
(3) the difference in osmotic pressure between the 
bulk solution and the air-solution interface of the 
wetting film. The osmotic pressure is determined 
from van’t Hoff’s simple law: P=nkT, where n 
is the total number of ions per cc. At the air inter- 
face of the film, 2 is obtained from the Maxwell- 
Boltzmann equation: im; exp (—z:ef/kT) 
or n=)>_n;K* at the air interface. The osmotic 
pressure in the bulk solution is P= > ,kT. The 
difference in osmotic pressure, tending to make 
the film thicker, is then 

di nikT(K*i—1). 

Combining these pressures at equilibrium, when 
their sum is zero, 


For thorium chloride-potassium chloride mix- 
tures, Eq. (10) becomes after subtituting n;=c,N 
for 7=1, 2, 3: 


Citcetcs S00(d—B)hg 
C2 coNRT 


Cy C3 
=—Ki+K+—K=2). (11) 


Co 


K is then determined by solving this equation for 
its root between 0 and 1, since 0<K <1 if the 
zeta-potential is positive. The value computed for 


do 


film, 
ection 
he air 


)hg/2, 
drical 
and 
en the 
of the 
mined 
1ere 
inter- 
xwell- 
b/kT) 
motic 
The 
make 


when 


(10) 


mix- 


(11) 


on for 
if the 


K corresponds to a certain y,, at the air interface, 
which is the minimum y in the film. This means 
that if the potential y due to the zeta-potential is 
sufficiently high so that enough ions are attracted 
into the film to make the osmotic pressure differ- 
ence large enough to overcome the pressure 
—(d—B)hg/2, then the value of y at the air 
interface is ¥m, which can then be computed 
according to Eq. (10). Obviously ¢ must be 
larger than the computed y,, for this to be possi- 
ble. If ¢ is found experimentally to be smaller 
than Wm, then there can be no wetting film 
attributed to the zeta-potential, because the 
osmotic pressure is insufficient to form it. 

In proceeding with Eq. (9), substitute 2A as 
expressed in Eq. (11), and make the following 
further simplifications of notation: 


B=C3/C2; 
b=(c1/c2)y®; 
Thus Eq. (9) becomes 
(d6/dy) = —B/[a+b}}. (13) 


In the region of convergence, this equation can 
now be expanded into a Taylor series: 


do —B b 3-175? 
dy Va 2a 2!2?X\a 


(2n—1)(2n—3)--+1 
n!2" 


This series, represented by S(y), is Eoetints 
convergent in the region where |b| <|a|, and in 
this region is therefore integrable termwise. To 
determine whether the limits between which we 
wish to integrate lie within this region, we ex- 
amine first the limit yo=K, determined by 
Eq. (11), which can be rewritten as 


(15) 
Ce 


Compare this with the requirement for con- 
vergence of the series, and we see immediately 
that |b| =|a| when y=K, and |b| <|a| when 
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y<XK. Therefore the series is not convergent at 
y=K, although it is convergent for any y less 
than K. Let yp = K —e and let be taken so small 
that any error introduced by using this yo instead 
of K becomes infinitesimal. If the values of y in 
the region of integration are restricted to 
0<yi<y<~yo, which is also a requirement of the 
physical assumption |~m|<|¢| of the theory, 
they are all in the region of convergence, the 
series expansion is valid, and Eq. (14) is integ- 
rable termwise. 

This integration would be extremely laborious 
to perform, and fortunately is unnecessary, since 
an idea of the value of the integral can be had by 
a relatively simple method. The series S(y) in 
Eq. (14) will have its minimum value in its 
convergent region at yo when b/a has its greatest 
value, slightly less than unity. Its maximum 
value is at yi, where b/a is very small, and S(y:) 
is slightly less than unity if y; is appreciably 
smaller than yo. The value of the series can be 
said to be S(yo) <S(y) <1 for all values of y in 
yiSy<yo. The sum of the first 15 terms of S(yo) 
is 0.78053, and the 16th term is —0.13825. In a 
uniformly convergent, alternating series, the 
value of the remainder of the series after summa- 
tion of the first m terms will be of the same sign 
as, and less in absolute value than the (7+1)st 


term. Thus we can say® 

vo dy 40 S(y)dy vo dy 

f —> >0.64 f 

However, the spread between the right and left 
members of (16) can be narrowed considerably by 
considering the integration to be performed in 
two steps: first between y; and an intermediate 
value y;, and then between y; and yo, where y; is 
defined by 


jal 100° Wi <Vi<Vo. (17) 


9 Should it occur to the reader that the limits of the inte- 
grals appear to be reversed for a correct statement of the 


inequality, it should be mentioned that Js " dy//a is 


- negative when y:<yo. This is demonstrated below when it 


is resolved into the Legendrian form of an elliptic integral 
(Eq. (22)), and it is also required by the ph gy nature of 


—B 
the problem since Ar =x, —x9=—— 


must be positive. Actually the et of sign is immaterial, 
since the proper sign of the square root may be taken to 
give a positive Ar. 
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Now the minimum value of the series S(y) in the 
interval y; to y;isat y= y,, where |b|/|a| =1/100, 
and evaluation of the first three terms tells us 


vis vi d 
vi Va a 


Then from (17) and (19) we have 


vo d vo S(y) “id 
| 


(18) 


> 0.995 


/a 
vo dy 


+0.64 
vi V/a 


(19) 
It is perhaps worthwhile to point out that the 
reasoning between Eqs. (13) and (19) merely 
constitutes a rigorous proof, when the numerical 
values for the solutions of present interest are 
applied, as will be shown, that the y® term in 
Eq. (9) can be neglected without introducing a 
large error. It also makes possible the calculation 
of two numbers between which Ar must lie, and it 
will be shown that these numbers are relatively 
close together.. 

It remains to express fdy/\/a in a form which 
can be evaluated. 


(1-4 J/(1+4), 


The integral then reduces to the Legendrian form 
of an elliptic integral, and we have 


dy dt 
—B 


TABLE IV. Calculation of parameters, and y;, K and ym. 


(20) 


KCI concentration 


X 10°) 


d— 
(1/x)(X 108) 


(mv) 
¢ (experimental) 
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ZETA- POTENTIAL IN MILLIVOLTS 


300 
OT iN ANGSTROMS 
Fic. 1. 


2By/u 


To convert this into the form in which Legendre’s 
tables of elliptic functions can be used, change 
the variable again by the substitution 


q°t? =q* cos? ¢+sin? ¢, 


and we have 


udy —H 


=——_[F (a, 
q 
in which sin?a=(g—1)/qg. Thus we have 
evaluated Ar as being confined between two 
readily evaluated integrals: 


do 


[1—sin? sin? ¢ |} 


F(a, go) ], (22) 


H 
(a, — F(a, $1) ]>Ar 


H 
{ F(a, — F(a, 


+0.64{ F(a, ¢o) — F(a, 
COMPUTATION OF Ar 


(23) 


In computing Ar it is necessary to determine 
the various constants and the limits of integration 
defined in the above equations. These depend on 


80 __ 
| 
| | 
where 
|| 
i= 
Let 
|| 
= 
10-3 10-4 10-5 
1 1 1 
€2( X 10°) 1002 102 12 
€3( X 10°) 1006 106 16 
P| 0.995,95 0.995,91 0.995,90 
95.62 287.5 649.1 
1.004,00 1.039,22 1.333,33 
1.005,0 1.048,9 1.354,0 
1.109,28 2.073,50 10.125,0 
4.339,25 1.444,81 0.639,34 | 
4.436,11 1.712,96 1.121,36 
76.972° 54.283° 26.903° 
0.623,00 0.267,60 0.066,05 
0.633,46 0.267,92 0.066,06 
0.456,56 1.317,08 2.717,19 
+11.73 +33.83 +69.79 
+67 +57 +17 


nge 


23) 


ine 
ion 
on 


the concentrations of the different ions (¢1, C2, €3), 
the capillary rise (z), and the density in air of the 
solution (d—8), all of which depend on the 
particular solution. The temperature was 25°C in 
all cases. The fundamental physical constants 
required are taken from the tabulation of Birge.'® 
The concentration of the thorium ion (c¢;) is the 
same in all solutions: 1.0 10-* M. The concen- 
tration cz is the sum of the concentrations of all 
ions having a charge of +1, in these instances, of 
potassium and hydrogen, the latter caused by the 
ionization of water. It is assumed that all ions 
derived from water are univalent, and that since 
the specific conductance of the water used was 
about 0.8X10-*, the concentrations of positive 
and negative ions from the water are each 
2X10-* M (considering A, the equivalent con- 
ductance to be 400). For example, cz for the solu- 
tion of 10-4 M KCl is 10210-*. The concen- 
tration c3 is the sum of the concentrations of ions 
of —1 charge, namely, the chloride from both 
KCl and ThCly, and hydroxide from water. As an 
illustration, cz for the solution with 10-4 KCI is 
106 10-°. The capillary rise h is taken as 10.862 
cm in each case, and the dielectric constant D 
is 78.49. 

In Table IV are given numerical values of the 
various quantities required for the computation 
of Ar, according to the definitions given above. 
Solution of Eq. (17) gives y; (actually the y; 
values give || /|a| somewhat smaller than 0.01) ; 
solution of Eq. (11) gives K; and pm, the po- 
tential at the air-solution interface, is given by 
m=eblm/kT. The experimental values for ¢ are 
given for comparison with y,,, and it should be 
noted that for the 10-5 7 KCI solution the ex- 
perimental value of ¢ is less than ~m, which is 
assumed to be the minimum value in the film. 
This means that the theory is not applicable to 
this solution, since the thorium ion has evidently 
reduced the zeta-potential to such a low value 
that it is no longer an important factor in 
determining the film thickness. For this solution, 
the film would presumably be negligible and 
attributed to van der Waals or other adsorptive 
forces operating at short range of a few angstroms. 
The inapplicability of the Langmuir theory to 
this solution should not be regarded asa weakness 


*R. T. Birge, Rev. Mod. Phys. 13, 233 (1941). 
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of the theory, but rather as the achievement of 
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the result intended by Langmuir when he pro- 
posed the addition of thorium to reduce the zeta- 
potential to such a low value that it would not 
create a thick film. 

Using the values of the parameters in Table IV, 
the thickness of the film has been calculated for 


‘various values of the zeta-potential, and the 


curves described have been plotted in Fig. 1. The 
curves A, B, and C are for the solutions con- 
taining 10-*, 10-*, and 10-5 M KCI, respectively. 
The values plotted are the “upper bound” of Ar, 
the left member of (23). The “lower bound,” or 
right member of (23) is at most 9 angstroms less 
than the given Ar values. This demonstrates that 
the approximate value obtained for the hyper- 
elliptic integral in Eq. (9) is satisfactory for 
application of the Langmuir correction. The ex- 
perimental values obtained for ¢ are shown on the 
curves A and B by circles, and the half-circle 
marked C, on the ordinate axis, represents the 
experimental ¢ obtained for the solution with 
10-5 M KCL. This point is really not off the curve 
C, because these curves should all be considered 
as starting at the origin. Thus as ¢ increases from 
0 to Ym, Ar is always 0, according to the Langmuir 
theory, and when ¢ increases beyond pm, Ar 


TABLE V. Wetting film thicknesses. 


KCI concentration: 10-*M 


+67 + 

n 2.608,645 

0.073,635 0 

0.623,00 

y(=K) 0.633,46 

t 0.225,343 0.233,231 0.863,086 1 

sin 1.0000 0.998,11 0.518,401 
90.00° 86.475° 31.225° 0° 

F(a, o) 2.900,60 2.630,64 0.572,46 0 

Fla, — Fla, $1) 0.269,96 2.328,14 2.900,60 

Upper bound, Ar 22.5 194.5 242.3 

(angstroms) 
Lower bound, Ar 14.4 185.4 233.0 
KCI concentration: 1074 

4 +57 + 

n 2.219,295 

0.108,684 0 

0.26760 


vi 
yo(=K) 0.26792 
t 0.583,762 0.584,156 0.807,316 1 
sin 1.0000 0.999,670 0.726,831 0 
90.00° 88.530° 46.622° 0° 
Fla,¢) 2.019,34 1.975,28 0.876,40 0 
Fla, $0) — F(a, $1) 


0.044,06 1.142,94 2.019,34 
Upper bound, Ar 10.7 277.2 9.7 
(angstroms) 
Lower bound, Ar 6.8 271.9 483.4 


| | 

22) 
ave 
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increases according to the abscissas of the curves 
shown. 

Examples of the computations are given in 
Table V, for the “‘upper’’ and “lower bounds”’ of 
Ar (i.e., the left and right members of (23)). These 
are given for the experimental values of ¢ and for 
infinite ¢ for the 10-* and 10-4 M solutions. The 
upper and lower bounds for Ar for infinite ¢ for 
the 10-5 M solution are 615 and 610 angstroms, 
respectively, determined in the same manner asin 
Table V. 


HYDROLYSIS OF THORIUM CHLORIDE 


The application of the Langmuir theory given 
above is based on the simplifying assumption 
that thorium chloride is not hydrolyzed. This is 
unquestionably an over-simplification, as there is 
no doubt that it does hydrolyze appreciably in 
water, because its solutions are acidic. It is 
possible that there may exist in solution the ions 
Th+4, Th(OH)+*, Th(OH).+, and Th(OH)3s*, and 
there may be other ions of thorium as well. 
There is no quantitative information in the 
literature on this hydrolysis. Experiments are in 
progress in this laboratory to determine the 
extent of hydrolysis, and until they are com- 
pleted, there is no way of making a precise 
computation of the thickness of the wetting film 
taking this into account. However, it is evident 
that Ar would be little affected even by relatively 
complete hydrolysis, because the concentration 
of potassium chloride is greater than that of the 
thorium chloride. Without attempting to discuss 
in detail the necessary alterations in the mathe- 
matical treatment to take hydrolysis into ac- 
count, it can be stated that the effect is to shift 
the curves A, B, and C in Fig. 1 a short distance 
in the direction of greater concentration. Thus C 
would be shifted slightly in the direction of B, 
and B slightly in the direction of A, while A’ 
would be scarcely changed since the potassium 
chloride concentration is so very much larger 
than the thorium chloride concentration. With 
complete hydrolysis to Th(OH)3* or Th(OH),, 
the computation would revert to the simple 1-1 
electrolyte case in which c=[KCI]+4 xX 10-°, (or 
c=[KCI]+6X10-, taking into account the ions 
from the water). It is therefore evident that even 
if the information were available to take the 
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hydrolysis into account, the curves A, B, and C 
would be essentially unchanged. 


APPLICATION OF THE WETTING FILM CORREC- 
TION TO THE SURFACE TENSIONS 


The wetting film correction factor to be ap- 
plied to the apparent relative surface tension is 
(r—Ar.)/(r—Aro) where r is 0.0136 cm, the 
radius of the capillary. The surface tension of 
10-* M KCI with 10-§ M ThCl, is taken as the 
average of the two values given by Jones and 
Frizzell, and Dr. W. A. Ray kindly measured the 
surface tension of 10-4 and 10-> M KCI with 
added thorium for the writer. His previously un- 
published results are, given in Table VI, along 
with the Jones-Frizzell measurement, and are 
referred to 10-§ M ThCl, instead of pure water. 
Since the zeta-potential of the reference solution 
has not been measured, the Langmuir correction 
is here applied with reference to the 10-° A/ 
solution, and the corrected surface tensions are 
relative to this solution. If the minimum is 
eliminated, the surface tensions of the 10-* and 
10-4 M solutions should be greater than (or 
equal to) that of the 10-* 1 solution, but it is 
seen in Table VI that on the contrary the 
application of the correction serves to lower the 
minimum. These results then do not support 
Langmuir’s suggestion that there is no minimum 
in the curve. 

It is worthwhile to inquire what the film 
thicknesses would have to be for the Langmuir 
correction to eliminate the minimum, in order to 
see whether experimental error in the zeta- 
potential measurements might account for the 
failure of the theory. Rearranging the equation 

r—Ar, 
r—Aro 


where op is the corrected relative surface tension, 
and gg is the observed apparent relative surface 
tension, we have 


OR OR 
Ar,.——Aro= (: 
OR OE 


Since ¢r/cz is very nearly unity, we can write 


oR 
Ar.—Ary= (: 
3 
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In order for the relative surface tensions to be 
exactly 1, i.e., to give a horizontal line for the 
surface tension-concentration curve, it is found 
that Ar for the 10-5 M KCI solution must be 
13.6 angstroms greater than Ar for the 10-4 
KCI solution, and the latter must in turn be 122 
angstroms greater than Ar for the 10-* MM KCl 
solution. If ¢ is +67 for the 10-* solution, then 
it would have to be +66 for the 10-4 M solution, 
and +91 for the 10-> M solution. It does not 
seem probable that an experimental error of this 
magnitude could have been made in measuring 
the zeta-potential of the 10-° VM solution. 

Since the failure of the Langmuir theory in 
this instance cannot be attributed to experi- 
mental error in the zeta-potential measurements, 
consideration should be given to the possibility 
that the charge at the air-solution interface of 
the wetting film may be large enough to play a 
more important part in determining the thickness 
of the wetting film than does the zeta-potential at 
the silica interface. If instead of the electric 
potential having its minimum value at the air 
interface, it has a comparatively high value, there 


TABLE VI. Corrected surface tensions. * 


Apparent 
relative 
surface 
tension 


Langmuir 
correction Corrected relative 
factor surface tension 


KCl 
conc. 


Ar 
(angstrom) 


(1.00000) 
0.99979 
0.99976 


10-5 d 0 
10-4 277 
10-8 194 


0.99980 
0.99986 


* Referred to the solution of 10-5M KCl—10-6M ThCl«. This refer- 
ence solution has the same apparent relative surface tension to that of 
water as a solution, namely 0.99995. 


might be a thick film due to this potential even 
though the zeta-potential has a very low value. 
As was pointed out by Jones and Frizzell,’ the 
Langmuir theory entirely neglects the effect of 
electric charges and the resulting potential at the 
air interface on the wetting film thickness. It 
may be that if the theory is extended to take this 
potential into account and when information 
about its magnitude is obtained, then the mini- 
mum in the surface tension curve will be satis- 
factorily explained as owing to the Langmuir- 
type wetting film with these solutions, just as it 
was with the same solutions without added 
thorium. 
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Acetone has been photolyzed in a flow system at several temperatures in the pressure range 


0.5 to 11 mm with exciting radiation of wave-length 2650-2900A and 2537A. The free methyl 
radicals formed in the decomposition have been determined as lead tetramethyl. Approxi- 
mately one methyl radical is formed for every four molecules of carbon monoxide, the actual 
proportion depending on temperature and pressure. Relatively more free radicals are formed 
at 100°C than at lower temperatures and more are formed at lower pressure than at high. 
A process involving the existence of a long lived metastable state for photoactivated acetone 
has been suggested to explain the results. The proposed scheme is found consistent with most 


of the published work on acetone photolysis. 


INTRODUCTION 


UBSEQUENT to the absorption of ultra- 
violet light, acetone molecules decompose 
either by rearrangement into the two principal 
products, ethane and carbon monoxide,'? or by 
fission into methyl and acetyl radicals which in 
turn give ethane, biacetyl, and carbon mon- 
oxide.*-> According to some authors it is unneces- 
sary to consider the rearrangement into ultimate 
molecules under any circumstances since the 
main features of the reaction are accounted for in 
terms of free radicals. The absence of any definite 
information on the relative importance of the two 
possible paths of decomposition has given rise to 
numerous and sometimes contradictory sugges- 
tions for the detailed mechanism of the total 
reaction. To explain different sets of data, for 
example, various authors have had to suggest 
several possibilities for the disappearance of the 
acetyl radicals: formation of biacetyl homo- 
geneously, through bi-molecular and tri-molecular 
reactions, and heterogeneously on the walls; 
decomposition of the acetyl radicals both hetero- 
geneously and homogeneously, and both uni- and 
bi-molecularly ; and several combinations of these 


* Present address: National Distillers Products Corpora- 
tion, Cincinnati, Ohio. 

** Present address: Monsanto Chemical Company, 
Clinton Laboratories, Oak Ridge, Tennessee. 

1 R. Spence and W. Wild, Nature 138, 206 (1936); J. 
Chem. Soc. 352 (1937), 590 (1941). 

20. D. Saltmarsh and R. G. W. Norrish, J. Chem. Soc. 
455 (1935). 

3 J. J. Howland and W. A. Noyes, (a) J. Am. Chem. Soc. 
66, 974 (1944); (b) 63, 3404 (1941). . 
a oa} S. Taylor and C. Rosenblum, J. Chem. Phys. 6, 119 

5 E. Gorin, J. Chem. Phys. 7, 256 (1939). 
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possibilities." * *-* It has also been suggested that 
the acetyl radicals disappear by reaction with 
methyl radicals to produce ethane.® 

A clear-cut decision concerning the behavior of 
the radicals seems to be impossible until at least 
the relative importance of the ultimate molecule 
reaction has been settled. Since the free alkyl 
radicals can be quantitatively converted into 
volatile lead compounds by allowing them to pass 
over lead mirrors,’ it should be possible in the 
absence of chain reactions, as is the case in 
acetone photolysis," to determine exactly what 
fraction of the products of reaction arises from 
free radicals and what fraction, if any, comes 
from rearrangement of the photoactivated mole- 
cules. We have undertaken such an investigation. 


EXPERIMENTAL 
Procedure 


Carefully degassed acetone was introduced as a 
vapor from a storage bulb of the substance and 
was allowed to flow through a section of quartz 
tubing surrounded by a Hanovia S 2537 lamp, the 
radiated energy from which is at least 90 percent 


6D. S. Herr and W. A. Noyes, J. Am. Chem. Soc. 62, 
2052 (1940). 

7A. O. Allen, J. Am. Chem. Soc. 63, 708 (1941). 

8H. W. Anderson and G. K. Rollefson, J. Am. Chem. 
Soc. 63, 816 (1941). 

9S. W. Benson and G. S. Forbes, J. Am. Chem. Soc. 65, 
1399 (1943). 

10(a) F. A. Paneth and W. Hofeditz, Ber. 62, 1335 
(1929); (b) F. A. Paneth, W. Hofeditz, and A. Wunsch, 
J. Chem. Soc. 372 (1935). , 

1 E, I. Akeroyd and R. G. W. Norrish, J. Chem. Soc. 890 
(fost: J. A. Leermakers, J. Am. Chem. Soc. 56, 1899 
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in the 2537A line. The illuminated acetone was 
allowed to pass over a mirror of lead deposited 
immediately adjacent to the lighted section of 
tubing (except in a few cases where the mirror 
was further down stream). The mirror was of 
such length as to give full lead transport” at a 
pressure of 4.4 mm. The gas stream, after it had 
passed over the mirror, was sent through a trap 
kept at liquid nitrogen températures. The gase- 
ous products not condensed in this trap were 
removed by a Toepler pump and were analyzed 


for carbon monoxide, and occasionally also for 


ethane, in a semi-micro gas analysis apparatus." 


Lead Analysis 


The determination of the lead alkyl in the cold 
trap, by dithizone titration, was carried out as 
previously described." In later experiments 
(serial numbers above 200), a modified end-point 
procedure was used which was more convenient 
and probably more dependable. In this method, 
an excess of dithizone is added to the sample 
whose lead content is to be measured and the 
excess determined by titration with a standard 
lead nitrate solution. This procedure eliminated 
the blank corrections necessitated by loss of 
dithizone in the direct titration. 


RESULTS 


The experiments were conducted with inlet 
pressures of 0.5, 2.1, 4.4, and 11.5 mm of acetone 
and with the illuminated region at 0°, 45°, 60°, 
80°, and 100°C. In Table I are recorded the 
results from the experiments in which the exciting 
radiation was 2537A. 

The volume of gas recorded in Table I is in 
mm corrected to S.T.P. The amounts of lead are 
in micrograms. It is assumed, in computing 
methyl production, that each lead atom trans- 
ported from the mirror is combined with four 
methyl groups. One microgram of lead is there- 
fore equivalent to 0.432 mm of CH; or 0.216 mm? 
of potential ethane, at S.T.P. 

Some of the radicals, however, are not re- 
covered because they disappear on the walls of 
the tube before they reach the mirror; but the 


 M. H. Feldman, J. E. Ricci, and M. Burton, J. Chem. 
Phys. 10, 618 (1942). 

8K. W. Saunders and H. A. Taylor, J. Chem. Phys. 9, 
616 (1941). 
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number that escapes detection in this way can be 
computed through the following considerations. 
The methyl radicals disappear in a first-order 
process.'* 24 We assume that this process is the 
same in both illuminated and dark regions. 
During illumination of acetone, however, radicals 
are released at a practically constant rate because 
the fraction of molecules that decomposes is 
small, actually something in the neighborhood of 
10-*. The rate at which methyl] radicals accumu- 
late as acetone passes through the illuminated 
section is 


(1) 


TABLE I. Photolysis of acetone with 2537A radiation. 


per 
min- 
ute, 
in 
mm’? CO/(}CHs+C2He) CHsX1.12/CO 
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144, May, H. A. Taylor, and M. Burton, J. Am. Chem. 
Soc. 63, 249 (1941). 
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where (CH3;) is the concentration of methyl 
radicals at any point, ¢ is the time for the gas to 
travel from the inlet to that point and ko and k, 
are the coefficients of the zero-order formation of 
radicals for a fixed set of conditions and of the 
* first-order disappearance, respectively. Assuming 
the pressure gradients to be constant and the 
expansion that of a perfect gas, ¢ is given by the 


formula,!5 
t=x/v—ax?/2Pv. (2) 


In this equation, x is the distance the gas has 
travelled, v is the streaming velocity (approxi- 
‘mately 3X 10* cm sec.~! in our experiments), P is 
the entrance pressure, and a is the pressure 
gradient along the tube. The value of ki was 
calculated from half-life experiments performed 
in the same apparatus,” being 63 sec.~! when the 
radicals are formed at 45° and 87 sec.-! for 
photolysis at 100°C. The coefficient ko was then 
evaluated from the integrated form of Eq. (1), or 


(CHs) =[Ro/ki JL1—exp (—hit) J, (3) 


by substituting the amount of methyl recovered 
when the mirror was immediately adjacent to the 
illuminated section of quartz tubing. The total 
amount of methyl which would be recovered at 
this point if no disappearance had occurred in the 
lighted region is then simply kot. The calculation 


based on the experimental figures shows that . 


approximately 12 percent of the radicals disap- 
peared before they were carried away from the 
lamp. The fraction of radicals thus disappearing 
in an illuminated zone for which ¢ is not too large 
is moreover nearly independent of both the 
temperature and the pressure. The final column 
in Table I gives ‘‘corrected”’ ratios of CH;/CO; 
these are simply the experimental figures in- 
creased by 12 percent to allow for the radicals 
lost before they‘treached the mirrors. 

The column in Table I headed CO/($CHs; 
+C.H¢) constitutes a kind of materials-balance 
for our experiments. If the sole products’ are 
carbon monoxide, ethane, and lead tetramethyl, 
the ratio should be unity. That the average of all 
the ratios is close to one, with relatively few 
ratios deviating by more than 20 percent from 
this figure, is evidence that the principal products 
are indeed the three substances named. The 


16 F, Paneth and W. Lautsch, Ber. 64, 2708 (1931). 
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TABLE II. Photolysis of acetone at 100°C and 2.1 mm, 
using 2650-2900A. 


per 


CO/(}CHs+C2He) (CHs/CO)X1.12 


@ No lead mirror. 
+ Mirror used, but lead not determined. 


deviations probably reflect the difficulties of 
working with micro-quantities of products rather - 
than the formation of important amounts of other 
substances such as biacetyl. 

Table II contains the experimental results of 
runs where the activating radiation lies in the 
region 2650-2900A. A General Electric H-6 lamp 
was the source of light and the light was filtered 
by a chlorine-mercuric chloride cell.!* The runs 
were made at 100°C with an entering acetone 
pressure of 2.1 mm. 


Completeness of Radical Recovery 


Before interpreting these data, it is necessary 
to consider whether the lead mirrors actually 
convert all the methyl radicals in the photolyzed 
vapor into lead tetramethyl. Paneth, Hofeditz, 
and Wunsch'™ have shown that over 80 percent 
of the methyls formed by decomposition of lead 
tetramethyl are converted by contact with a lead 
mirror into lead-alkyl. It seems probable that 
every free radical striking the metal reacts 
with it. In their experiments, however, Paneth, 
Hofeditz, and Wunsch tested the behavior of the 
radicals when all of the reactant was decomposed 
and only inert carrier gas and radicals came into 
contact with the mirrors. In our experiments, 
the radicals accompany a very large excess of 
undecomposed ketone, and a question arises 
concerning the completeness of radical recovery 
under these circumstances. 

All the methyl radicals actually present are 


16 E. J. Bowen, J. Chem. Soc. 76 (1935). 
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removed from the gas by a 2 cm mirror” and they 
must be converted either into lead tetramethyl 
or into ethane, for no other methyl derivatives 
are formed insofar as we can judge from the 
materials-balance. It is not easy to picture a 
specific catalytic effect of acetone in causing the 
methyl radicals adsorbed on the surface of lead 
to combine giving ethane rather than lead tetra- 
methyl. Were such an effect in operation, how- 
ever, it would be pressure dependent; yet at 100° 
a fivefold increase in pressure leaves the CH3;/CO 
ratio nearly unchanged. One further possibility 
must be considered, namely, a conversion of 
methyls adsorbed on the lead surface into some 
product other than lead alkyl by reaction with 
acetone molecules. The observation, however, 
that a lead mirror has no effect on carbon 
monoxide production or on the ratio, CO/($CH; 
+C.H¢), means that the methyls cannot dis- 
appear by reaction with acetone. We must 
conclude that all the free radicals reaching the 
mirror are changed to lead alkyl. 

In connection with the calculation of methyl 
production, we have assumed a maximum alkyla- 
tion of lead. The formation of less highly 
alkylated compounds of lead would mean that 
we have over-calculated the methyl concentra- 
tions and that our estimates of the fractions of 
molecules decomposing by way of free radicals 
are maxima. 


DISCUSSION 


The results summarized in Table I, contrary 
to all previous opinion, indicate that an im- 
portant part of the acetone decomposition at 
2537A occurs as a rearrangement. If the reaction 
proceeded wholly through free radical formation, 
the CH;/CO ratios, under conditions such that 
all the acetyl radicals decompose, i.e., at high 
temperature, would have a value of 2. The ratios 
in all cases are much below this figure, ranging 
from 0.2 to 0.6. At high temperatures, where it 
may be assumed that two mirror-active particles 
are produced per molecule splitting into free 
radicals, the fraction, a, of decomposing acetone 
molecules producing free radicals, would be 
numerically one-half of the mean “corrected” 
CH;/CO ratios of Table I. At lower temperatures 
this fraction, a, would be somewhat greater than 
half of the tabulated ratio, but it cannot be 


greater than the ratio CH;/CO even if the 
acetyl radicals formed are assumed not to decom- 
pose at all. It is evident, therefore, that much of 
the carbon monoxide is formed unaccompanied 


‘by free radical formation, and that the fraction 


of the total decomposition proceeding through 
free radical split is low, and of the order of 0.3. 

These data make it necessary to explain, in 
terms of a reaction that proceeds partially by 
way of a rearrangement, the numerous reported 
experimental results connecting both quantum 
yields and products with such variables as 
pressure, temperature, and intensity. It is not 
enough to propose two alternative modes of 
decomposition of activated acetone, one to give 
radicals, one to give molecules. Such a simple 
proposal (without introducing the assumption of 
very effective recombination of the primary 
radicals, a question which will be discussed 
below) cannot explain, for example, quantum 
yields below (1—a), the fraction of molecules 
undergoing the rearrangement reaction. A some- 
what more complex scheme is required. Such a 
scheme has been suggested by Spence and Wild! 
for the photo-decomposition at 3130A. The Spence 
and Wild mechanism is consistent with most of 
the established facts at 3130A. A modification 
and amplification is here suggested for the reac- 
tion at 2537A. 


The Spence and Wild Reaction Scheme 
The reactions may be summarized as follows: 
(1) 
CH;COCH;*—CH;CO+CHs, (2) 
CH;COCH;*+M-—CH;COCH;**+ M, (3) 
(4) 
CH,COCH +0, (5) 
CH;COCH;**—>CH;COCHs. (6) 


Briefly this scheme suggests that photo-activated 
acetone molecules either decompose in a rela- 
tively short time to give radicals or are converted 
(by collision) into a metastable state from which 
they may escape by collisional deactivation or 
by rearrangement into molecules.'!? The lifetime 


17 In photolytic reactions, the principle of least motion as 
proposed byRice and Teller (J. Chem. Phys. 6, 489 (1938)), 
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Fic. 1. Relationships in acetone photolysis at 100° 
(schematic). A: free radical split; B: ultimate molecular 
rearrangement; C: total decomposition. 


of the metastable molecules produced by radia- 
tion of wave-length 3130A, is, according to the 
estimates of Spence and Wild, in the neighbor- 
hood of 10-* sec. The life of the initially activated 
molecules, however, must be very short as shown 
by the complete diffuseness of the absorption 
spectrum.!8 

The efficiency of various types of molecules in 
bringing about the transition of excited molecules 
‘to the long lived state cannot be predicted. Nor 
can their efficiency in inducing rearrangement or 
deactivation be accurately estimated prior to 
actual experiments. It is only possible to suppose 
that the effectiveness will depend on the type of 
molecule and on the temperature. 

If one rejects consideration of collisional deac- 
tivation, the low quantum yields under ordinary 
conditions must be ascribed to recombination of 
the initially formed radicals.!* The reaction-com- 
plex of methyl and acetyl according to this view 
must undergo a transition to acetone with much 
greater efficiency than a similar association of 
two methyl radicals to give ethane or two acetyl 
radicals to give biacetyl. Under some conditions, 
the quantum yields are as low as 0.2, which would 
mean recombination of 80 percent of the primary 
radicals. The specificity of the reaction, (i.e., the 
efficiency compared with the total of all other 


according to which the lighter ketones would not be ex- 

ted to rearrange directly into stable products, is not 
strictly applicable as, indeed, Rice and Teller themselves 
point out (page 496). Consequently we are basing our 


discussion on the experimental findings. 
-18R, G. W. Norrish, H. G. Crone, and O. D. Saltmarsh, 
J. Chem. Soc. 1456 (1934). 
19H. S. Taylor and J. C. Jungers, Trans. Faraday Soc. 33, 
1353 (1937); suggest that part of the effect is due to the 
reaction 
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possible reactions of these radicals) would mark 
this photolysis as unusual. It would be a striking 
circumstance not merely that the photo-activated 
molecules should be so unstable that all of them 
decompose in ~ 10-8 sec., but that the radicals 
would endure long enough to enter this process 
almost exclusively. Although the higher energy 
content of the initially excited molecules as 
compared with the united radicals, which have 
only thermal energy, would lead to greater 
instability of the former, the yield needed to 
explain the result would have to be so great that 
some alternative way of accounting for them 
seems desirable. In our own experiments, the 
carbon monoxide yield is not affected by the 
mirrors beyond a matter of 5 to 10 percent, and 
that result is possibly attributable to other 
variations. The removal of methyl radicals from 
the gas mixture with consequent elimination of 
union of methyl and acetyl radicals does not 
change the basic course of the reaction. In other 
words the recombination of methyl and acetyl 
radicals is unimportant in our experiments, al- 
though it may occur under other conditions.” 


EFFECTS OF EXPERIMENTAL VARIABLES 
ON THE REACTION 


At a given wave-length, the proportion of 
excited acetone which undergoes a free radical 
split and the proportion that rearranges, will be 
a function of the acetone pressure and of the 
temperature. The probable relations are shown 


1.0 


Fic. 2. Relationships in acetone photolysis at room 
temperature (schematic). A : free radical split; B: ultimate 
molecular rearrangement; C: total decomposition. 


2 Cf. K. Faltings, Ber. 72, 1207 (1937), who found 
acetone to be formed with a quantum efficiency of ~0.1 
when ethane is photolyzed in the presence of carbon 
monoxide. 
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schematically in Fig. 1 for high temperatures 
(100°C) and in Fig. 2 for room temperature. 
At very low pressures the decomposition will be 
completely effective in producing free radicals, 
for in the absence of collisions, the unstable 
molecules can do nothing else but split in this 
way. As the pressure increases, more and more 
of the activated molecules will be sent, as a 
result of collision with other gaseous molecules, 
into a long lived state. At high temperature, it is 
assumed that practically all the molecules in this 
photo-activated metastable state decompose be- 
cause of their great internal energy, so that the 
total quantum yield (or the efficiency of the total 
decomposition) should remain substantially at 
unity at all pressures. At lower temperatures 
however, the rearrangement passes ultimately 
through a maximum because of the increasing 
importance of collisional deactivation. At those 
temperatures, therefore, where quantum yields 
are appreciably less than unity, the quantum 
yields would probably continue to fall with 
increasing pressure, and would approach zero as 
a limiting value. 

While the experimental quantum yields do not 
reach unity at ordinary temperature, it is true 
that several investigators have found that 
quantum yields increase rapidly below 50 mm 
pressure both at 2537A and at 3130A. Yields at 
the longer wave-length are smaller than at 
shorter under comparable conditions because the 
less energetic molecules are more apt to be 
deactivated by collision before decomposition. 
Despite the drop in quantum yield, experimental 
values of the ratio CH;/CO are about the same 
at the two different wave-lengths (Tables I and 
11); apparently, a molecule which has absorbed 
a quantum at A=3130A must have both a 
smaller chance of a split into radicals in a given 
time and a correspondingly greater chance of 
entry into a metastable state. . 

The effects of pressure and temperature on the 
experimental ratio, CH;/CO, may be discussed 
on the basis of these curves. That this ratio 
should fall with increasing pressure follows at 
once from the effect of pressure on the two 
competing paths of decomposition as shown on 
Figs. 1 and 2. In Fig. 3 the average of the observed 
ratios, CH3;/CO, (Table 1) is plotted as function 
of pressure for the three principal temperatures 


CH; 
CO 


P(mm) 


Fic. 3. Observed ratio, CH;/CO, as function of pres- 
a i 2537A. Squares: 0°C; circles: 45°C; triangles: 


studied. As expected according to the proposed 
scheme, the ratio rises with decreasing pressure, 
at each temperature, and may be expected to 
reach the theoretical value of 2 as a limit. The 
actual shape of the curves drawn, however, does 
not have great significance, since the uncertainty 
of each point is rather large. It would be difficult 
to predict how the effect of pressure should 
change with the temperature, since a rise in 
temperature would presumably increase the effi- 
ciency of both paths relative to deactivation. 
The figure shows that the ratio increases with 
rising temperature, especially at the higher pres- 
sures, in spite of the greater efficiency of the 
rearrangement reaction at higher temperature. 
This state of affairs may result partly from a 
greater proportion of the initial reaction giving 
free radicals (even though the collision frequency 
at a given pressure is increased) but more 
especially from an increase in the average number 
of methyls arising from a molecule undergoing 
the free radical split. There seems to be con- 
siderable evidence that acetyl radicals become 
almost completely unstable at 100°C, so that two 
methyls would be formed from each acetone 
decomposing in this way, whereas an average of 
less than two per molecule would be formed at 
lower temperatures. The fraction, a, of the total 
decomposition proceeding through the free 
radical path would therefore be half of the ob- 
served CH;/CO ratio at 100°C, and somewhat 
greater than half this ratio at lower temperatures, 
depending on the probability of reactions of 
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CH;CO other than decomposition yielding CH; 
as a product. 

The effect of narrower tubes (a decrease in size 
of tube results in increase of the ratio CoHg/CO)** 
appears to be to leave the free radical process 
substantially unchanged but to increase the total 
quantum yield by inducing rearrangement of the 
metastable molecules. At 100°C, where the total 
quantum yield is about unity anyway, the effect 
of changes in the geometry of the apparatus 
should be a minimum. 

There is some question of how the quantum 
yields are affected by a diminished light intensity. 
If there is any effect at all, it is to increase the 
yield at low intensities. While Spence and Wild 
did not discuss this effect, it can be explained on 
the basis of some reaction of activated molecules 
with products, such as radicals, other activated 
molecules and final products of the reaction. It 
seems very likely that biacetyl in particular 
would account for the decrease in yield with 
rising intensity. The transfer of the excitation 
energy from acetone to biacetyl would be ac- 
companied by a corresponding increase in light 
emission from the latter. Since fluorescence is 
observed” even above 100°C, there should be 
some intensity effect on quantum yields at 100° 
although this would be difficult to observe. 

There have been numerous reports of the 
effects of foreign gases on the reaction. The effect 
of carbon dioxide,” oxygen,” hydro- 
gen,**5 propane,’ deutero-biacetyl,™ and_ nitric 
oxide® have been investigated. Some of these ex- 
periments have been cited in-support of a com- 
plete radical mechanism ; but it is in these cases 
especially that the course of the reaction is very 
seriously changed, as, for example, when nitric 
oxide is introduced. In each case, however, it 
seems necessary to consider specific reactions 
between acetone molecules in the metastable 
state and the admixed gas. Thus the extremely 
efficient production of iodides suggests two im- 
portant reactions, namely, 


CH;COCH;**+1.—-2CH3I+CO, (7) 


21 W. A. Noyes and F. C. Henriques, J. Chem. Phys. 7, 
767 (1939). 

22 P, Fugassi, J. Am. Chem. Soc. 59, 2092 (1937); F. O. 
(i . E. Schildknecht, J. Am. Chem. Soc. 60, 3044 

1 

23W. West, J. Am. Chem. Soc. 57, 1931 (1935). 

24. —D.S. Herr, M.S. Matheson, and W. D. Walters, J. Am. 
Chem. Soc. 63, 1464 (1941). 
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(8) 


which might be considered together with the 
reaction 


CH;CO+1.2CH;COI +I, (9) 


suggested by Gorin to explain the dependence of 
iodide production on the concentration of iodine.® 


Other Radical Reactions 


Since all of our measurements of free radical 
concentrations are based on the idea that the 
disappearance reactions are the same in the 
illuminated region as beyond, it seems desirable 
to consider certain reactions which might lead to 
an extra fast clean-up in the illuminated region. 
One of the possibilities arises from Gorin’s con- 
clusion® that acetyl radicals are stable at 100°, 
disappearing at that temperature not so much 
by decomposition as by reaction with methyl 
radicals, 


This reaction would be not only an additional 
source of ethane but also an outlet for methyl 
radicals. It may be argued that some, perhaps a 
large fraction, of these radicals disappear in this 
way at 100° and so are not picked up by the 
mirrors. Were this bimolecular reaction con- 
tributing significantly to the over-all changes, 
however, it would be much more important at 
high radical concentrations than at low. A five- 
fold change in concentration of the radicals, as 
measured by the rate of carbon monoxide pro- 
duction, causes nothing greater than a 30 percent 
change in CH;/CO ratio at any temperature and 
only a 20 percent change at 100°. We must con- 
clude that reaction (10) is not important in our 
experiments. For similar reasons, the reaction, 


(11) 


cannot have an important role in our experiments. 
The union of methyl and acetyl radicals to 
reform acetone, moreover, can have no effect on 
the interpretation of results in Tables I and II 
because such a union can affect only the quantum 
efficiency of the total process and not the distri- 
bution of products. The formation of biacetyl 


(12) 


cannot account for small values of CH3/CO but 
on the contrary should lead to ratios greater than 


| 
_ 
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unity. While biacetyl seems not to make up more 
than a minor part-of the products, each biacetyl 
molecule removes two potential molecules of 
carbon monoxide from the products and its effect 
would be to increase the net CH;/CO ratio. 

Another type of reaction which might be of 
importance within the illuminated zone but not 
outside involves hot radicals, that is, energy-rich 
radicals emerging from decomposing acetone 
molecules. The behavior of such radicals from 
propionaldehyde has been considered by Garrison 
and Burton.** The chief effect of such hot radicals, 
may be to induce rearrangement of other acetone 
molecules to give stable products. If the initial 
reaction of excited acetone molecules, however, 
produced in every case two such hot radicals, 
quantum yields should be three rather than <1 
as observed. Furthermore, even if all the initial 
splitting gave rise to such specially active radicals, 
there would still be more CO produced than 
corresponds to recovered radicals. 

We have omitted from the reaction list those 
reactions which follow the initial free-radical 
split. There are numerous possibilities for forma- 
tion of ethane, biacetyl and carbon monoxide. 


These possibilities have already been extensively 
discussed in the literature and the present work 
affords little additional information concerning 
them. 


CONCLUSION 


Most of the photochemical studies of acetone 
can be explained satisfactorily (in the absence of 
free radical determinations) either by a com- 
pletely free radical mechanism or by a mecha- 
nism such as that of Spence and Wild. The 
particular interest in the present work is that the 
free radicals have been counted and the actual 
proportion of reaction proceeding through radical 
intermediates has been estimated quantitatively. 
Decision as to the secondary free radical reactions 
must as before depend on kinetic studies. Like- 
wise, opinions concerning the behavior of the 
metastable state must be based on experimental 
determination of quantum yields, fluorescence 
and nature of the products. Nevertheless, the 
Spence and Wild scheme apparently is consistent 
with published data, for 2537A. 


25 W. M. Garrison and M. Burton, J. Chem. Phys. 10, 
730 (1942). 
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COMPARISON OF PROPIONALDEHYDE AND 
ACETONE PHOTOLYSIS 

The experiments of May, Taylor, and Burton" 
indicated that equal numbers of free radicals are 
released when propionaldehyde or acetone vapors 
are photolyzed under comparable conditions. It 
was argued that this identity, assuming the same 
quantum yields, means that if one photolysis is 
an exclusively free radical process the other must 
also be such. It was later found, however, that 
the only mirror-active particles formed in the 
photolysis of propionaldehyde are ethyl radicals,”* 
and there can be only one such radical coming 
from one molecule of the aldehyde. But an 
acetone molecule may give either one or two 
methyl radicals and on the average will give 
something between. The experiments of May, 
Taylor, and Burton were conducted at a tempera- 
ture (about 55°C)? where the acetyl radicals 
produced from acetone may be considered to be 
at least partially if not completely decomposed. 
The observed approximate equality in mirror 
activity then cannot be consistent with a 100 
percent free radical mechanism for acetone, if the 
quantum yield of the two photolytic decomposi- 
tions is approximately the same. On the basis of 
such an assumption, the observation involved 
would be further interesting evidence of the 
reality of the results here reported, which indi- 
cate so strongly that the photolysis of acetone 
does not proceed exclusively through free radical 
formation. The observation of May, Taylor, and 
Burton would mean either that both processes 
involve partial free radical formation or that at 
least the acetone photolysis does. A determina- 
tion of the quantum yields, then, of the two 
decompositions under the conditions of the ex- 
periments of May, Taylor, and Burton, would be 
of interest; if the quantum yield of the pro- 
pionaldehyde photolysis is not greater than twice 
that of the acetone it follows that some of the 
acetone must decompose without production of 
free radicals. The data of Garrison and Burton 
for the ratio C:H;/CO in propionaldehyde 
photolysis*® show that about 45 percent of the 
propionaldehyde molecules producing carbon 
monoxide also produce ethyl radicals (allowance 
being made for about 25 percent disappearance of 
radicals in the illuminated zone). 


26. May, Thesis, New York University, March, 1941. 
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few halftone cuts. The Editorial Board will not hold itself 
responsible for opinions expressed by the correspondents. 
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later than the 15th of the month preceding that of the issue in 
which the letter is to appear. No proof will be sent to the 
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not be made and no reprints will be furnished free. 


Some Observations on the Theory 
of Contact Angles 


HAROLD M. SCHOLBERG AND W. W. WETZEL 
Minnesota Mining & Manufacturing Company, Saint Paul, Minnesota 
August 27, 1945 
F a plane solid surface is dipped into a large liquid 
surface, a meniscus is formed. The shape of this menis- 
cus is governed by the law: 


1/Ri+1/R2=pgy/o1,. (1) 


where the R’s are the principal radii of curvature of the 
liquid surface, p is the density, g the acceleration of gravity, 
y is the height of the meniscus at any point, and a; is the 
free surface energy of the liquid surface. 

If the work necessary to stretch the liquid surface and to 
raise the water is calculated, the following expression is 
obtained : 


(2 —ah?)3/2}+ oih cot r(1 -*), (2) 


where h is the height of the meniscus at the point of con- 
tact with the solid, I is the angle of tilt of the solid surface 
with the horizontal, and: 


a?=pg/2e1. . 


The free energy available to do the work is given by the 
expression : 


(3) 


where oq is the free surface energy the solid-air interface, 
and o,; is the free surface energy of the solid-liquid inter- 
face. These two expressions are not equal. However, if 
Eq. (3) is differentiated with respect to I and the resulting 
equation set equal to zero, the derivative of h with respect 
to I is obtained. If Eq. (2) plus an unknown factor, f(h, T), 
is also differentiated with respect to I and set equal to zero 
and the expression for dh/dI substituted in this second 
differentiation, the necessary terms to be inserted in (2) 
to make it equal to (3) are obtained. Since this is an ex- 
treme condition it is the necessary and sufficient condition 
for equilibrium. The additional terms inserted in expres- 
sion (2) are exactly the terms obtained by calculating the 
reversible work done on the slide to rotate it from the 
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angle Io, where h=0, to the angle I. If we write 


cos agen =a constant, 
1 


the final equation takes the form 


sin 3 
+. pein a 
where a is the angle of contact and is a function of h and I. 

Experimentally Eq. (5) is obeyed well within the limits 
of error for measurements on paraffin against water. This 
is all that is necessary to show, for plane solid surfaces 
intersecting a liquid in such a manner that one of the 
principal radii of curvature of the meniscus is infinite, that 
differences in the advancing and receding contact angles 
are impossible providing aq is constant. 

If the two principal radii of curvature of the liquid 
surface are both finite, an additional variable but not an 
additional equation of condition is introduced; the system 
then becomes univariant instead of non-variant and a 
difference between the advancing and receding contact 
angles becomes possible. 

A more complete paper on this subject will appear later. 


The Constitutional Difference Between a- and 
6-Derivatives of Naphthalene 


J. A. A. KETELAAR AND G. W. VAN OOSTERHOUT 


Laboratory for General and Inorganic Chemistry, 


University of Amsterdam, The Netherlands 


August 21, 1945 

HE difference with respect, e.g. to the bond distance, 
electric dipole moment, etc., for chlorine, bond to a 
saturated carbon atom, or to a carbon atom carrying 
multiple bonds as in the substituted ethylene and benzene 
derivatives was explained as due to a partial double bond 
character in the latter case.! In order to see whether the 
rather crude method of approach? used before would be 
also valuable in explaining finer constitutional details, an 
analogous calculation was set up for a- and 6-chloronaph- 

thalene, respectively. 

The complete results will be published elsewhere, here 
only some conclusions will be given. For the energy differ- 
ence between the ground state and the first excited state, 
corresponding to the long wave limit of the U.V. absorp- 
tion we calculated 1.15028 and 1.16298 for the a- and 
B-derivative, respectively, as compared with 30845 cm™ 
and 31130 cm~!? found experimentally. The correspond- 
ence is quite remarkable. The double bond character for 
the C—Cl bond was calculated as 14.17 percent and 12.98 
percent, respectively (12.6 percent for monochlorobenzene). 

For the electric dipole moment of both molecules we 
may calculate (using an effective electronic charge‘ that 
gives the right experimental value for the a-compound) 
Experimentally: Parts’ 1.59 resp. 1.72; 
Nakata® 1.50 resp. 1.57 for a- and 8-chloronaphthalene. 
Both the sign and the magnitude are correctly given by the 
theory. . 
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For negative atoms with p, electrons other than chlorine, 
e.g., oxygen and nitrogen atoms, qualitatively the same 
result will be obtained indicating a larger double bond 
bond character for the a-derivative. According to current 
conceptions’ this would lead to a lesser basic strength for 
a-naphthylamine as compared with the corresponding 
B-derivative. This conclusion is corroborated by the experi- 
mental evidence. 


a-naphthylamine 
B-naphthylamine 


K=0,836 10-8 
K=1,29 10-”. 


The higher amount of double bond character as calcu- 
lated for the a-derivative as compared with #-chloro- 
naphthalene, etc., is thus found to be in complete accord- 
ance with various experimental evidence. We feel justified 
to conclude that the molecular orbital treatment in the 
crude form used here is able to account even for minor 
details of molecular constitution. 


1J. A. A. Ketelaar, Rec. trav. chim. 58, 266 (1939). J. Sherman and 
J. A. A. Ketelaar, Physica 6, 572 (1939). 

2G. W. Wheland and L. Pauling, J. Am. Chem. Soc. 57, 2086 (1935). 

3H. G. de Laszlo, Proc. Roy. Soc. A111, 355 (1926). 

4L. E. Sutton, Trans. Faraday Soc. 30, 789 (1934). J. A. A. Ketelaar, 
Rec. trav. chim. 58, 311 (1939). 

5 A. Pacts, Zeits. f. physik. hemie B10, 264 (1930). 

6 N. Nakata, Ber. 64, 2059 (1931). 

*L. Pauling, The Nature of the ae Bond (Cornell University 
Press, Ithaca, New York, 1939), p. 192 

8 N.F. Hall and M. R.S. Sprinkle, J. Am. Chem. Soc. 54, 3469 (1932). 


An Absolute Method for the Determination of 
the Area of a Fine Crystalline Powder 


D. HARKINS 
University of Chicago, Chicago, Illinois 
AND 


GEORGE JURA 
Universal Oil Products Company, Chicago, Illinois 
July 12, 1945 


N a letter in this Journal! with the above title the writers 
presented a very accurate method for the determination 
of the area of a very fine solid by a calorimetric method. 
The results obtained by this method were found to be in 
excellent agreement with those given by the method of 
Emmett and the application of the theory of Brunauer, 
Emmett, and Teller.? In this work a crystalline powder is 
first outgassed and then exposed to the saturated vapor of 
any liquid which gives a duplex film on the surface of the 
solid. This film has the same surface energy as the liquid 
itself and with water is about 25A, or nearly 8 molecules 
thick. The area of the film-coated powder is found from the 
total heat in ergs evolved when the powder is dropped into 
and dispersed in the liquid in a sensitive calorimeter. 

In a letter to The Journal of Chemical Physics Hans M. 
Cassel* expresses agreement with the theory of this method. 
He fears, however, that when water is the liquid, the ad- 
sorption of the duplex film will give adhesion between the 
crystals which ‘‘may cause inseparable clusters to form.” 
To the writers it appears that Cassel has two objections to 
the experimental application of the method: (1) the forma- 
tion of clusters, which do not disperse in the liquid, and 
(2) the effect of capillary condensation, which reduces the 
area of the duplex film. 


THE EDITOR 449 

As to (1), Cassel’s references to Tomlinson‘ and Stone 
give evidence against his own idea, in the more complete 
work of Stone. In considering (2), Cassel states “this point 
has been overlooked” by Harkins and Jura. Had he read 
our papers carefully he would have known that this state- 
ment is not true, since it is discussed in J. Am. Chem Soc. 
66, 926-7 (1944). 

From a careful consideration of this particular problem 
we are able to state that with powders of the degree of 
fineness thus far used this method is far more accurate 
than any other; its only disadvantage is its difficulty. The 
error due to capillary condensation is so small that it is 
unimportant. 

1. While the adsorption is taking place there is mostly 
point contact between the crystals, since only about one- 
tenth of the space in the mass of crystals is occupied by 
them, while nine-tenths is empty space. This condition 
makes capillary condensation a very slow process. 

2. To adsorb a duplex film the mass of crystals is ex- 
posed to saturated vapor for a number of hours. The energy 
of immersion when the film is adsorbed at 96 percent of 
saturation, where only half as much water is adsorbed, 
is only 3 percent higher than at the saturation pressure. A 
decrease as small as 3 percent with an increase of vapor 
pressure from p/po=0.96 to p/po=1 is possible only if the 
area remains essentially constant in this region, since a 
decrease of about this magnitude is to be expected to occur 
from the considerable thickening of the film itself by ad- 
sorption without any effect of capillary condensation. 

3. The fact that the effect of capillary condensation is 
shown by Fig. 1, which gives the heat of emersion of a 
sample of anatase (covered with various thicknesses of 
water) as determined by (a) direct calorimetric measure- 
ments and by (b) calculation of the heat from the adsorp- 
tion isotherm. The agreement is excellent, although the effect 
of capillary condensation would be to make the calorimetric 
value low and the adsorption value high. In the adsorption 
method the energy due to capillary condensation is in- 
cluded, making the value high if capillary condensation is 
appreciable, while it lowers the area and thus the calorimet- 
ric heat. Thus, the agreement between the two methods 
shows that capillary condensation is so small as not to 
affect the area obtained by the calorimetric method. 
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4. Concrete experimental evidence that capillary con- 
densation does not play an important role in adsorption 
on non-porous surfaces is shown by the work of Bangham 
and Mosallam, who showed that the adsorption of benzene 
on cleaved mica sheets was the same even at the highest 
vapor pressures, whether the sheets were separated by fine 
wires or bound tightly together. 

Other evidence which shows that the effect of capillary 
condensation is so small that it may be entirely neglected 
will be presented in later papers, since it is too lengthy to 
be given here. 

IW. D. Harkins and G. Jura, J. Chem. Phys. 11, 430 (1943). 

2S. Braunauer, P. H. Emmett, and E. Teller, J. Am. Chem. Soc. 60, 
att a Cassel, J. Chem. Phys. 13, 249 (1945). 


4G. A. Tomlinson, Phil. Mag. 6, 695 (1928). 
5 Stone, Phil. Mag. 9, 610 (1930). 


Birefringence of Polyisobutylene in Solution 
and the Frenkel Theory of 
Macromolecular Rupture 


Joun REHNER, JR. 


Esso Laboratories, The Standard Oil Development Company, 
Elizabeth, New Jersey 


August 27, 1945 

HE problem of the shape of dissolved macromolecules 
has long held the attention of many investigators 
because of its important bearing on the questions of macro- 
molecular structure, branching, deviations from the Ein- 
stein viscosity formula, and the shear stability of the solute 
under certain conditions of viscous flow. In a recent paper! 
Zvetkov and Frisman have published experimental data 
showing how the birefringence An of various samples of 
polyisobutylene, dissolved in a hydrocarbon solvent, de- 
pends on the velocity gradient g. Their curves (Fig. 11 in 
their paper) show that for a polyisobutylene sample of 
sufficiently high molecular weight, the dependence An =f(g) 
is already curvilinear at the smallest observable velocity 
gradients (g~100 sec.~') and that there exists a definite 
value of g at which the slope of the curve increases more or 
less abruptly. The magnitude of this value of the gradient 
is higher the lower the viscosity of the solution, and there- 
fore the lower the hydrodynamic force acting on the 

molecule in flow. 

Recognizing that linear macromolecules are ruptured 
under the influence of sufficiently high shear forces, Zvetkov 
and Frisman have attempted to account for the shape of 
their curves on the basis of a theory, recently published by 
Frenkel,? in which it has been proposed that a critical 
velocity gradient g; exists at which the coiled chain becomes 
rod-like in its central portion, and above which rupture will 
occur. The theoretical value of g; is related to the viscosity 
n, the effective radius r of a monomeric chain unit, the 
length / of a monomeric link, and the degree of polymeriza- 
tion N, according to 


(1) 


For a polyisobutylene sample of average molecular weight 
132,000 as determined by the Staudinger formula, Zvetkov 
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and Frisman applied Eq. (1) and found g,2500 sec.~', 
in good agreement with the point of steep ascent in their 
corresponding curve. 

As this agreement might be regarded as supporting the 
theory of rupture, it is felt advisable to call attention to 
two possible sources of error. (1) Zvetkov and Frisman 
were apparently unaware that polyisobutylene solutions 
do not obey the Staudinger formula’ over a wide range in 
molecular weight. The true molecular weight of their 
particular sample, if determined by osmotic pressure meas- 
urements, would be of the order of 10 times as large as the 
value they employed, and substitution of this in Eq. (1) 
would reduce g; by roughly 100-fold, in disagreement with 
their experimental observations. (2) In deriving Eq. (1) 
Frenkel has assumed: that, in the process of rupture, the 
central chain link in the macromolecule is extended in 
accordance with Hooke’s law. This is not only physically 
unlikely and quantitatively erroneous but, what is more 
important, his final relationships leading to Eq. (1) 
appear to depend for their “catastrophic” form on the 
analytical linearity of Hooke’s law as well as the assumed 
validity of the latter for bond extensions as large as Al~/,‘ 
even though Frenkel admits that under conditions of 
rupture Hooke’s law can no longer be applied. These facts 
suggest that the agreement between experiment and theory 
indicated by Zvetkov and Frisman is fortuitous, and that 
in order to account for the point of a steep ascent in the 
birefringence curve it is probably necessary either to find 
a different cause or to develop a theory of rupture which 
places greater emphasis on the force law governing the 
extension of the valence linkages. Results obtained in 
these laboratories® with ordinary and fractionated samples 
indicate that the shear degradation characteristics of 
polyisobutylene in solution display a more complex be- 
havior, in their dependence on temperature, solvent, and 
molecular weight, than that suggested by Eq. (1). 


(1945) Zvetkov and E. Frisman, Acta Physicochim. U.S.S.R. 20, 61 
2J. Frenkel, Acta Physicochim. U.S.S.R. 19, 51 (1944). 
3P, J. Flory, J. Am. Chem. Soc. 65, 372 (1943). . 
4 The writer is indebted to Dr. R. F. Tuckett for an interesting dis- 
cussion of this point. 
5 J. I. Wasson, unpublished data. 


Inversion of Low Angle Scattering Data for 
Particle Size Distribution 
S. H. BAvER 
Department of Chemistry, Cornell University, Ithaca, New York 
August 21, 1945 
HE purpose of this note is to indicate how the expres- 
sion for the observed intensity as a function of the 
scattering angle, which involves an integral over the de- 
sired particle size distribution, may be inverted ina straight- . 
forward manner. 
The equation generaly used for the ratio of scattered. 
to primary intensity is of the form: 


p°N(R) exp (1) 


where a?=.221 (Biscoe and Warren), k= (4x sin 0/2)/d, R 
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is the radium of gyration of the particle, N(R) is the num- 
ber distribution in radii, p is the number of electrons per 
particle, and K contains fundamental and apparatus 
constants, as well as the polarization factor. Since 


= M(R), 
the mass distribution, and 
where dp is the mass density and d is the electron density, 
Eq. (1) reduces to: 
I/Ip=K']} R'M(R) exp (—a*?Rk?)dR. (2) 

0 
To obtain M(R), substitute x =a?k?, and y= R?; the equa- 
tion can then be inverted according to well-known methods 
of Laplace transforms. In practice, instead of so deducing 
distributions, ranges in particle radii are determined by 
measuring slopes of the curve obtained by plotting log 
(I/Io) against k®. However, Eq. (1) follows from an ap- 
proximation introduced in order to simplify. the process of 
reducing the data.! 

For spherical particles, ‘‘widely and randomly separated 
to produce incoherency in the radiation from individual 
particles,”’ the ratio of secondary to primary intensity is:? 


where 
(x) = (3/x?)((sin x/x)—cos x). 


The replacement of ¢? by the exponential which appears in 
Eq. (1) is unnecessary since it is possible to invert Eq. (3). 
After correcting the original data for polarization, finite 
slit length, etc. the expression for the intensity may be 
rewritten, 


F(2k)/K" =k*(I/Io corrected 


=K" (4) 


where 
K" = 
since 
J 3/2") = (2/xx)*((sin x/x)—cos x). 


Employing the generalized Neumann integral,* the above 
equation reduces to: 


"Ta(2kR cos q)dqdR, (5) 


F(2k) = (2/n) J 
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Eq. (5) may be rewritten in the form of the Schlémilch 
integral equation after changing the order of integration: 


F2k)=(2/n) "G(2k cos )dq, (6) 


where 
G(z)= J (s=2kcos p). (7) 


The Hankel transformation’ permits one to express the 
number distribution in terms of an integral over G(s), 


R°N(R)= (8) 


whereas G(z) may be expressed as an integral involving 
F(2k):5 


m2 9 
G(2k) = F(0)+2k J sin (9) 


Admittedly, the procedure in reducing the data is some- 

what involved. The program consists of the following steps: 

(a) Convert the photographic density to (1/J0)corrected- 

(b) — F (2k); plot the derivative of F with respect to 2k, vs. 2k. 

(c) R — the argument in the latter by (2k sin g), and thus get a 
mily of curves, one for each value of 2k, as functions of g. 

(d) Integration according to Eq. (9) gives G(2k). 

(e) Replot G(2k) as a function of 2; integration according to Eq. (8) 

gives N(R) 

That the original Eq. (3) may be applied only to a 
dilute suspension of spherical particles should be stressed. 
It is not valid where the separation between the particles 
approaches the order of magnitude of the mean radius. 
Therefore, there is some question as to the significance of 
particle size distributions based on this method which 
have been presented in the literature. In the case where 
the holes are dimensionally comparable to the particles, a 
straightforward analysis and inversion procedure have 
been formulated by Professor P. Debye.* What is significant 
then is the fluctuation of the product of electron density at 
two points separated by R, as these points scan the avail- 
able volume of diffraction. 

I wish to thank Dr. J. M. Richardson for discussions 
pertaining to the inversions used. 

1 For a simple discussion refer to Louis A. Pipes, J. App. Phys. 10, 
172, 258, 301 (1939). 

2 J. Biscoe and B. E. Warren, J. App. Phys. 13, 371 (1942). 

3G. N. Watson, Treatise on the Theory of Bessel Functions (The Mac- 
millan, New York, 1944), p. 150. 

4M. A. Stone, Linear Transformations ay Mathematical 
Society Colloquium Publications, 1932), Vol. XV, +p i 

5 E. T. Whittaker and G. N. Watson, Modern aie (Cambridge 


University Press, New York, 1940), p. 229. 
6 Private communication. 
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